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Preface 


The research in the field of advanced materials is significantly driven by the devel- 
opment of *smart" materials that are designed to react with environmental condi- 
tions and other stimuli. Smart materials that have attracted considerable attention 
in the last years have now become an intrinsic field in dealing with both societal 
challenges and significant engineering problems in various sectors such as energy 
and water consumption, recycling and engineering. 

The field of advanced materials benefits from multidisciplinary basic and ap- 
plied research in physics, chemistry, mathematical modeling and materials engi- 
neering. Technical challenges related to developing suitable smart materials and 
intelligent structures are applicable in a wide range of fields, such as sensors, spin- 
tronics, supercapacitors, biomedical applications (drug delivery and hyperthermia 
therapy), self-healing, flexible memories, construction technologies and other emerg- 
ing technologies. 

Carbon-based materials have gained major attention for replacing other con- 
ventional materials as a promising basis for developing novel smart materials due 
to their abundance and low cost, high-temperature stability and corrosion resis- 
tance, high electrical and thermal conductivity and enhanced mechanical proper- 
ties, particularly in the microscopic and nanoscopic scale. 

This book consists of chapters written by expert scientists and engineers from 
the international carbon-based materials community, who contribute important re- 
search work with emphasis to advanced applications. Contributions focus on novel 
applications of carbon-based smart materials, as well as on functionalization pro- 
cesses and multiscale modeling. 

In the first chapter *New class of graphene-based devices for the next generation 
of nonvolatile memories," the main scientific works concerning the properties of mem- 
ories based on graphene layers, graphene oxide, and reduced graphene oxide are re- 
viewed by showing the promising results in terms of writing time and endurance/ 
cyclability compared to existing technologies. Moreover, these materials can poten- 
tially allow developing a highly scalable complementary metal-oxide-semiconductor- 
compatible technology. This chapter suggests that graphene-based memory technol- 
ogy can be exploited for space applications, as it is being insensitive to radiation and 
easily reprogrammable. 

Second chapter is devoted to the research undertaken at the University of 
Birmingham on the functionalization of carbon-based materials using the active 
screen plasma (ASP) technology. The chapter entitled *Plasma surface activation 
and functionalization of carbon-based Materials" exploits the typical physical and 
chemical phenomena associated with the plasma technology, further optimized by 
the additional control and degrees of freedom offered by the ASP technology. In 
this way, the highly reactive species in the plasma were used to reduce graphene 
oxide films, activate the surface of carbon paper membranes or functionalize the 
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surface of carbon fibers. All the above seem to be great opportunities to create pat- 
terned surfaces using the ASP technology, in combination with other processes, 
with potential application in multiple fields such as functional materials, sensors 
and biosensors, electrical, electronic and optoelectronic devices. Without question, 
the ASP technology remains a valuable enabling technology for the surface func- 
tionalization of engineering materials, showing broader prospects of practical ap- 
plications than ever before. 

Modeling of carbon-based smart materials is the subject of the third chapter. A 
better understanding of the thermophysical behavior of carbon-based polymer 
nanocomposites at different length- and time-scales could pave the way of discover- 
ing novel smart materials and could be facilitated by proper materials modeling. 
Advanced modeling approaches for the estimation of nanocomposite properties are 
based on linking and/or coupling various model types that are applicable at differ- 
ent length- and time-scales. This chapter review could be particularly useful for a 
broader audience of readers, whose research interest is focused on modeling issues 
in carbon-based materials. Atomistic, mesoscopic and macroscopic simulation tools 
that allow exploring the main properties of nanocomposites are extensively pre- 
sented, and future perspectives of advanced modeling of carbon-based smart mate- 
rials are discussed. 

The fourth chapter entitled “Carbon nanotube-based materials for space appli- 
cations” summarizes the activities performed in three sequential European Space 
Agency projects, highlighting the introduction of carbon nanotubes (CNTs) into 
composite materials for space applications. The first activity aimed at proving the 
feasibility of using CNT skeletons as an alternative to conventionally filled resins, 
overcoming high viscosity issues, which limits the processability of materials. For 
the second activity, the aim was to scale up the nanotube structures and integrate 
them in carbon fiber-reinforced composites. More recently, activities were focused 
on selecting two applications with highest potential for performance improvements, 
taking into account all knowledge gained, followed by the development and testing 
of two demonstrators. This chapter presents the development steps and demonstra- 
tor manufacturing and testing for one of these applications, a carbon fiber- 
reinforced polymer (CFRP) optical mirror. For this mirror demonstrator, a complete 
assembly consisting of two CFRP skins with outer CNT layers bonded to a carbon- 
based honeycomb center was manufactured. The principle of producing a highly 
accurate optical mirror in CFRP-CNT technology is rated as very promising, with im- 
portant mass savings at required thermoelastic distortion values. 

The next chapter deals with molecular dynamics (MD) simulations of graphene- 
based polymer nanocomposites. The unique, advanced properties of graphene justify 
why it is widely considered as the material of the future. Within this chapter, graphe- 
ne’s great potential for new applications is highlighted, and the recent progress in the 
design of graphene-based polymer nanocomposites from both an experimental and a 
computational point of view is discussed. The basic concepts of molecular simulations 
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focusing particularly on the MD technique and simulation results for the structural, 
conformational and mechanical properties of a test nanocomposite system based on 
poly(methyl methacrylate) (PMMA) filled with two different types of graphene sheets 
are also presented: simple (i.e., unfunctionalized) graphene sheets and functionalized 
(i.e., graphene oxide). This study was driven by recent experimental reports according 
to which the addition of a small fraction of graphene sheets in a polar polymer matrix 
such as PMMA can lead to a remarkable enhancement of its elastic constants. 

The last chapter *Carbon from waste source for Li-ion battery" suggests that the 
accomplishment of *12 principles of green chemistry" is a fundamental goal to be 
pursued by means of designing/optimizing environmentally *conscious" processes, 
materials and devices to reduce their environmental impact throughout all the 
phases of their life cycle, thus benefiting the economy, protecting people and 
achieving the real goal of a truly sustainable world. This goal can be achieved by 
exploiting carbonized nanocellulosic binders and water-based papermaking cou- 
pled with high-performing green composite polymer electrolytes. 

The book is useful for everyone who is interested in applications of carbon- 
based smart materials and is especially oriented to the graduate students and 
young researchers to facilitate their introduction in this new area. The editors thank 
all the authors who contributed their work to this book and also acknowledge the 
superb assistance that the staff of De Gruyter has provided for the publication of 
this work. 
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Paolo Bondavalli 

1 New class of graphene-based devices 
for the next generation of nonvolatile 
memories 


1.1 Introduction 


Graphene is a one-atom thick layer of carbon atoms arranged in a hexagonal lattice. 
Graphene potentialities are attracting a lot of researchers to probe opportunities in a 
number of directions in the “more Moore" or “beyond CMOS” optics in order to iden- 
tify the new future technologies [1-3]. Another promising utilization of graphene and 
related nanomaterials is to fabricate nonvolatile memories (NVM) exploiting their 
*memresistive" behavior storing a value of electrical resistance in a permanent way. 
This happens when a current passing through the materials changes the level of re- 
sistance. Therefore, resistive memory exploits the change in the resistance of a mate- 
rial under the effect of an electric field as an information write/erase principle for 
nonvolatile data storage. The reading of resistance states is nondestructive, and the 
memory devices can be operated without transistors in every cell [4, 5], as for flash- 
type memories [6-10] (see Section 1.2.2), thus achieving a classic cross-bar structure. 
This kind of memories is called resistive random-access memory (RRAM or ReRAM) 
and is only one of the possible types of nonvolatile ways to store information in a 
permanent way. One of the most important advantages of these new classes of 2D 
materials is that these materials can be implemented in flexible electronics [11-15], in 
the form of one-thick atom layers as for graphene or in the form of layers of flakes 
of graphene oxide (GO) or reduced GO (R-GO), thereby reducing the final cost of the 
final device exploiting roll-to-roll fabrication [16]. Another great advantage of 
ReRAMs is their potential to implement them by exploiting only two terminals to 
work (two contacts and not three as a common transistor, which has drain, source 
and gate, e.g., flash-type memories), which could dramatically reduce the circuitry 
and allow to implement easily in 3D architectures by using the roll-to-roll fabrication 
technique. Potentially, this is applied in various fields such as health monitoring 
[17—21], intelligent packaging [22, 23], cards, labels, badges, value paper and medical 
disposables. Another potential field is the radiofrequency identification (RFID). This 
is a technology to electronically record the presence of an object using radio signals. 
Indeed, an innovative alternative pathway to reduce RFID costs and integrate a mem- 
ory chip to store data is to eliminate the silicon substrate completely, and produce 
RFID and memory on the same flexible plastic substrate as the antenna [24-26]. 


Paolo Bondavalli, Thales Research and Technology, Palaiseau Cedex, France 
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Thanks to graphitic layers, the antenna and chip can be built on the same low-cost 
substrate, and attachment costs can be removed. 


1.2 Graphene-based NVM 


As outlined in the International Technology Roadmap for Semiconductors 2011 sec- 
tion concerning Emerging Research Devices, “ultrathin graphite layers are interest- 
ing materials for macromolecular memories thanks to the potential fabrication 
costs that are considered as the primary driver for this type of memory, while ex- 
treme scaling is de-emphasized." The main drawback is related to the fact that 
memory operation mechanisms and the physics are still unclear and that a deeper 
research in this field is necessary to improve the comprehension of the phenome- 
non and the efficiency of the devices. These are not the same physical mechanisms 
exploited in memories based on graphene-related materials, which will be dis- 
cussed in the next paragraphs. 


1.2.1 Graphene and graphitic layers 


The first paper highlighting the utilization of graphene for resistive NVM exploiting 
two terminal structures was issued in 2008 by Stadley et al. at Caltech [27]. They 
reported the development of an NVM element based on graphene break junctions. 
These junctions were obtained by achieving two-terminal devices, transferring gra- 
phene sheets on SiO,/Si substrate and depositing metal electrodes on them using 
e-beam lithography process. After applying a specific voltage under ultravacuum con- 
dition (107 torr), they were able to create a break in the graphene layer (see Figure 1.1). 
In these papers, all the measurements were performed under ultravacuum to avoid 
environmental gas interferences, which strongly reduces the potential utilization of 
the device. 

After cycling the memories, they observed that the resistance values clearly 
changed as a function of the voltage (see Figure 1.1b): they moved from a high resis- 
tance state (off state) to a lower one as a function of the voltage (on state). This ef- 
fect was nonvolatile and the devices were cycled thousands of writing cycles (10°), 
showing long retention times (more than 24 h, only limited by measurements). 

The reproducibility of the switching behavior was very promising. One of the more 
coherent explanations was related to the fact that the conductance occurred along a 
small graphene ribbon that bridges the contacts. To better understand the underlying 
physical mechanism, researchers studied the time-resolved behavior of the switch from 
the OFF to ON states, shown in Figure 1.2a. The conductance I/V showed well-defined 
steps, with magnitude ~Go. Here Gg is 2e^/h, the conductance quantum, where e is the 
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Figure 1.1: (a) SEM image of the device before (left panel) and after breakdown (right panel). 

The arrows indicate the edges of the break in the graphene layer. (b) Repeatable programming over 
hundreds of cycles. Upper left panel: Voltage applied to the junction versus time. A ramp with a peak 
value of ~4 V corresponds to an ON pulse, while a ramp with a peak value of ~6 V corresponds to an 
OFF pulse. A small sawtooth-shaped read-out pulse is applied after each write to determine the 
junction conductance. Lower left panel: Current flow through the junction, with the low-bias 
conductance labeled above each read-out pulse. Right panels: Similar data taken after approximately 
500 intervening cycles, demonstrating the reproducibility of the switching behavior [27]. 


electron charge and h is Planck's constant. Since Go is the conductance of a spin- 
degenerate one-dimensional conductor, for example, a linear chain of gold atoms [28], 
observation of steps in the conductance suggests that the conductance states of a de- 
vice are likely multiples of highly transmitting quantum channels. Authors proposed a 
model for device operation based on the formation and breaking of carbon atomic 
chains that bridge the junctions (see Figure 1.2a, right figure). This mode and the re- 
sults seem to put in evidence the potential for multiresistive states. 

The same physical explanation can be evoked for the results obtained in 2009 
by Sinitski et al. [29] at Rice when they tested analogue carbon-based structures 
(see Figure 1.2b). The sensitive layer was consists of 10 nm thick graphitic disks 
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(a) 


(b) 


Figure 1.2: (a) Left side: Change in //V as a function of V pointing out the quantic behavior linked to 
the formation of atomic chains. Right side: Proposed physical mechanism at the base of the 
nonvolatile resistive behavior (formation of atomic carbon chains). (b) Tilted-view SEM images of 
two different devices before and after applying the voltage difference that is necessary to break 
the layers defined as Vpreak: The arrow shows the fracture across the stripe due to Vpreak [29]. 


grown by chemical vapor deposition (CVD) process. Researchers evaluated the volt- 
age that is necessary to break the layer. Calculating the Joule heating, researcher 
obtained that the breakdown voltage, the applied voltage difference necessary to 
break the layer, was 


Voreak[(C py Tpreak)/T| m l=A (1.1) 


where C, p, y are, respectively, the specific heat, the density of CVD-grown graphitic 
sheets and resistivity, Tpreax is the breakdown temperature of the graphitic layer, A 
is independent of the devices and l the length of the device. Therefore, the voltage 
necessary to break the layer is directly proportional to the length. This is very posi- 
tive because reducing the dimensions will allow to reduce the energy consumption 
during the writing step. 

In 2009, the same team from the University of Rice [30] exploited another 
approach and grew graphitic layers on a freestanding silicon oxide nanowire 
(see Figure 1.3a). Indeed they used the nanowire to preserve the mechanical integ- 
rity of the device after applying a specific voltage breaking the surficial graphitic 
layers. The team from Rice fabricated two terminal devices consisting of discontinuous 
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Figure 1.3: (a) Silicon oxide nanowire with graphitic layers grown on its surface and (b) physical 
principle at the base of the fabrication of nonvolatile graphitic memories exploited by Rice 
University [30]. 


5-10 nm thin films of graphitic sheets grown by CVD on nanowires having stable, re- 
writable, nonvolatile and nondestructive read memories with on/off ratios of up to 10” 
(see Figure 1.4a) and switching times of up to 1 ps (tested limit). 

All the samples were not tested, as in the previous examples, at ambient con- 
ditions but in ultra-vacuum condition (5 x 10? torr) and showed an extremely 
good stability of the performances after several days and also after x-ray exposure 
(Figure 1.4b). 

An interesting work by Wu et al. [31] in 2012 dealt with the electrical properties 
of graphene sheet suspended on a patterned indium tin oxide (ITO) electrode pair. 
They measured very interesting on/off ratio of 6 orders and retention time of at less 
10^ s in ambient conditions. In this case, data storage was achieved by applying 
voltage bias and rewritten after a simple heat treatment. Authors stated that the 
switch happens only in ambient atmosphere when oxygen changed the injection 
barrier at the contact. The degassing using a thermal process allowed coming back 
to the high conduction state. 
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Figure 1.4: (a) The logarithmic /-V behavior of a G-SiO; nanoribbon; the “write” or “erase” regions 
can be achieved through a voltage pulse in the noted regions. After writing or erasing, the high- or 
low-conductivity state of the system is recorded in the “read.” (b) Data retention by two other 
G-SiO; nanocable devices [30]. 


1.2.2 Nonvolatile resistive memories based on GO and R-GO 
oxide layers 


Other teams focused their study on GO. GO is a wide band gap material (6 eV) with 
potential for modulation, thanks to oxidation/reduction process providing tunabil- 
ity of the electronics, mechanical and optical properties. GO [32] is commonly ob- 
tained by the oxidation of graphite using the modified Hummer's method, where 
the long oxidation time is combined with a highly effective method for the purifica- 
tion of reaction products [33, 34]. A great advantage of this material is its cost, 
which is at least one-tenth the price of graphene flakes and it is stable in water sus- 
pensions. In GO, the memory configuration exploited is similar to cross-bar memories 
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Figure 1.5: Cross-bar memory devices based on GO [39, 40]. 


[35, 36] (vertical junctions composed by piling up the different layers, sandwiching 
the active one as metal insulator metal-type structures). We can mention the pio- 
neering work of He et al. in 2009 [37], where reliable and reproducible resistive 
switching behaviors were observed in GO thin films prepared by vacuum filtration 
method [38]. They fabricated Cu/GO/Pt structures showing an on/off ratio of about 
20, a retention time of more than 10^ s and switching threshold voltages of less 
than 1 V. The suggested physical principle at the origin of the switching effect was 
the desorption/absorption of oxygen-related groups on the GO sheets as well as 
the diffusion of the top electrodes. One of the major works on GO-based memories 
was issued in 2010 by Jeong et al. [39, 40]. This team demonstrated the nonvolatile 
effect on the resistance of a 70 nm layer thick GO in a layered structure composed 
of Al/GO/Al (Figure 1.5) in a cross-bar configuration. 

Jeong et al. were one of the first teams that pointed out the performances of these 
kinds of layers on flexible substrates. The researchers from KAIST demonstrated that 
using GO layers they were able to achieve flexible nonvolatile resistive memories 
with impressive performances under traction and compression [they kept three orders 
between on (LRS) and off (HRS) levels]. Concerning the physics explanation of the 
phenomenon, Hong et al. [41, 42] performed a deep analysis of the switching mecha- 
nism for this kind of devices and demonstrated that these structures had performan- 
ces dependent on the origin of the top contacts. For example, in case of Au-based top 
electrodes, there was no oxygen migration in opposition to Al electrodes. The effect 
of bottom contacts was also quite important. Indeed if the roughness was too high, 
the GO layer had such cracks and rough surface that the top electrode material could 
easily penetrate and then build filaments that hindered the switching operations. 
This drastically reduced the lifetime of the device (only 100 cycles) because of the 
formation of a permanent conductive path between the two electrodes. These results 
were highlighted by X-ray photoelectron spectroscopy measurements that pointed 
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out the permanent presence of Al near the bottom electrode as the main reason for 
failure. Evaluating the effect of the migration of oxygen versus the formation of Al- 
based filaments, Hong and coworkers stated that the formation of conducting fila- 
ments was a local phenomenon and that the oxygen migration was the dominant 
mechanism. Indeed, they discovered that the leakage current between the bottom 
and lower electrodes was dependent on the cell dimension. In fact when the cells 
were larger the effect of the oxygen migration or reduction on the conduction was 
enhanced proportionally to the surface. In order to avoid the failure mechanism 
linked to the formation of permanent conductive filaments through the material 
cracks, one suggestion is to exploit deposition techniques that allow a more uniform 
distribution of the deposited material. One technique is the deposition by spray coat- 
ing or heating the substrate to avoid the so-called coffee-ring effect [43, 44], which is 
responsible for the nonuniformity of the deposition [45]. This deposition could poten- 
tially prevent the formation of cracks linked to the roughness in the GO layers. 

Another interesting study was performed by Vasu et al. at the University of 
Bangalore (India) in 2011 [46]. They reported very simple unipolar resistive switch- 
ing device using ultrathin (-20 nm) films of R-GO with on/off ratios reaching five 
orders of magnitude. The thin films were formed at room temperature simply by 
drop-casting R-GO from suspension on ITO glass electrode, followed by aluminum 
or gold deposition (no difference was pointed out after the measurements). A very 
interesting result was also the switching time that could attain 10 us with an on/off 
ratio of 100. In this case, the formation of nanofilaments of carbon atoms was 
evoked as the main physical working principle. 


1.3 Other approaches to achieve nonvolatile 
memories using graphitic layers 


In the previous sections, we have discussed about the resistive NVM. Other teams 
have decided to use different approaches to exploit the properties of graphene- 
based layers to fabricate NVM and have exploited the hysteretic electrical behavior 
of graphene, or multilayered graphene (MLG), when used as channel in a field ef- 
fect transistor (FET)-type configuration. This is substantially based on the chemical 
modification of graphene under an electric field that will enhance environmental 
interactions. However, the more promising results, in terms of performances but 
also in terms of potential industrial exploitation, have been obtained in case of the 
fabrication of flash-type memory [47-48] using graphene or MLG as the floating 
gate (FG). Basically, a flash memory is an electronic nonvolatile computer storage 
medium that can be electrically erased and reprogrammed. In flash memory, each 
cell resembles a standard MOSFET [85], where the transistor has two gates instead 
of one. Indeed there is the traditional top gate, also called control gate and an FG 
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embedded in an oxide layer; hence, it is electrically insulated. In FG, which is lo- 
cated between the top gate and the MOSFET channel, any electrons placed on it, 
activating the top gate, are trapped there and will not discharge for many years. If 
FG holds a charge, it screens (partially cancels) the electric field from the top gate; 
therefore, it modifies the threshold voltage (Vr) of the MOSFET. This operation is 
nonvolatile and is reversible if another adequate voltage is applied on the top gate. 
The present technology used for flash-type memories employs doped polysilicon as 
material for the FG. One of the main advantages of graphene/MLG is the reduction 
of voltage to achieve a correct memory window, which is the shift of the threshold 
voltage of the transistor when the memory is switched from O to 1 state. In the pres- 
ent flash memories, a voltage difference of +20 V [49] is necessary to program/erase 
to achieve a memory window of 1.5 V, which is the industrial standard. Some works 
have demonstrated, as shown in the next paragraphs that thanks to the higher den- 
sity of states (DOS) of graphene compared to degenerately doped polysilicon this 
difference can reduced up to 6 Volts. Another advantage is the higher work function 
of graphene [50-52] that is directly linked to the larger barrier height between the 
FG and the oxide [53-56], which embeds the FG. This potentially allows achieving 
longer retention time, potentially more than ten years. Finally, the possibility to 
achieve device using 2D devices with a lower dimensionality that will potentially 
allow an easier 3D implementation and so a higher memory density. 


1.3.1 Graphitic-based nonvolatile memory using a transistor 
configuration 


One of the first paper highlighting the interest of using graphene-based FET for NVM 
was issued by Wang et al. at the Nanyang Technological University in Singapore in 
2010 [57]. Indeed they observed exactly the same phenomenon observed in case of 
FET with one carbon nanotube as a channel: unhysteretic phenomenon where they 
observed the formation of a loop of the current value between drain and source, Ias» 
as a function of the gate voltage, V,,, during the voltage sweep cycle from negative to 
positive values of V,, and back. This phenomenon is caused by the presence of mois- 
ture interacting with the transistor dielectrics. This effect, it has been demonstrated, 
can be removed in vacuum [58, 59]. Moreover the switching time in case of activation 
of the hysteretic phenomenon in FET using graphene is of the order of some seconds, 
which is too high for memory applications. 

In 2008, a pioneering work of Echtermeyer et al. of the Advanced Research 
Center in Aachen (Germany) and of A. K. Geim at the University of Manchester (UK) 
was performed adopting another point of view compared to the previous studies 
[60, 61]. Considering the difficulty to achieve nanoribbons with dimensions smaller 
than 10 nm to open a gap, they exploited another phenomenon. Indeed they 
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fabricated two-gate transistor structures as shown in Figure 1.6, where the graphene 
has a width of 10 pm and was covered by an SiO, layer. 
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Figure 1.6: (a) Schematic of a double-gated graphene field effect device (FED) and (b) current 
between drain and source as a function of the top gate voltage value [60, 61]. 


They observed that applying a voltage spanning from —5 to +5 V, for positive volt- 
age there was a strong enhancement in the current of 5/6 order, hence, a reduction 
in the resistance that was nonvolatile. Echtermeyer and coworkers suggested that 
considering that the measurements had been performed in ambient conditions, the 
water molecules were split in H* and OH" (see Figure 1.6) by the electric field and 
tended to attach to the graphene, creating graphene or GO (that are insulating ma- 
terials), passing through the highly porous silicon dioxide layer surface, and so 
drastically changing the conductivity of the channel in a nonvolatile way. This phe- 
nomenon allowed obtaining a strong on/off ratio, quite unusual for graphene (and 
not graphene nanoribbons) based FET. However, even if the results seem to be very 
promising, more studies have to be done on the switching time, which intuitively 
appears to be quite long considering the effects involved. We can add that the gra- 
phene layer has been obtained by mechanical cleavage, which is a technique that 
allows fabricating high-quality samples but that it is not industrially suitable. It is 
necessary to test this approach, for example, using CVD graphene in order to 
achieve a potential parallel fabrication of this kind of devices. It has also to be veri- 
fied if it is fundamental to obtain a single graphene layer or if more stacked layers 
of graphene can reach the same performances. 
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1.3.2 Nonvolatile flash-type memories based on graphene/ 
multilayered graphene 


Another potential utilization of graphene is FG in flash-type memories. An ex- 
tremely innovative approach has been proposed by Bertolazzi et al. of the Ecole 
Polytechnique de Lausanne (EPFL) in Switzerland. They were able to fabricate a 
flash memory using exactly the same architecture of the existing ones, exploiting 
only 2D materials [62]. 

Bertolazzi and coworkers were able to build a flash-type memory using MoS; as 
the channel and few graphene layers as the FG (see Figure 1.7), embedded in HfO; 
oxide, where the charge is trapped. It is clear that this kind of structure is extremely 
original but some performances especially considering the switching time (the low- 
est value is 100 ms) are not adequate. Some doubts can be raised on the fabrication 
technique, which includes three quite complex transfer steps and that are not suit- 
able at the moment for an industrial exploitation. Hong et al. at IBM in Watson 
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Figure 1.7: MoS;/graphene heterostructure memory layout. (a) Three-dimensional schematic view 
of the memory device based on single-layer MoS;. (b) Schematic view of the flash-type memory cell 
with a single-layer MoS; acting as a semiconducting channel and graphene contacts and multilayer 
graphene (MLG) as the FG [62]. 
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(USA) fabricated a flash memory [63] using some graphitic layers, grown by CVD on 
Si as FG. In this case, the technological approach appears potentially scalable, 
CMOS compatible and so industrially suitable (see Figure 1.8) compared to the pre- 
vious one. IBM researchers observed a potential retention time of 10 years, with a 
loss of only 896 of the charges in the FLG but also simulating the potential crosstalk 
between neighbor cells, who stated that this kind of memories showed negligible 
interference down to 10 nm (in case of common polysilicon FG, the interference 
raises dramatically under 25 nm). Briefly, considering the power reduction (as told 
previously, thanks to the higher DOS compared to common FG fabricated using 
polysilicon) and the increase in the storage density, Hong and coworkers estimated 
a potential reduction of 7596 in the operating energy of this kind of memory. 
However, in this work there are no data on the switching time and therefore we can- 
not compare their performances with existing memories. 
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Figure 1.8: Retention characteristics of GFM. (a) Retention measurement of MLG-FM showing only 
8% of charge loss in 10 years at room temperature. (b) Energy band diagram of MLG/SiO;/Si 
junctions before programming. (a) Energy band diagram after programming, which shows three 
possible mechanisms for charge loss during retention state: SE, FNT and DT. All three mechanisms 
depend exponentially on the barrier height (pB) between the work function of graphene and the 
SiO, tunnel oxide [63]. 
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1.4 Conclusions and potential applications 


Graphene-based memories is a completely new technology that needs some time to 
demonstrate its utility especially in the difficult field of NVM for space application. 
These memories could constitute a real breakthrough compared to existing technolo- 
gies also considering the dramatic potential reduction of energy consumption. These 
memories can be developed to bear harsh environments with a focus on radiation- 
resistant components. The stability (also in temperature), resilience, nonvolatility and 
on/off ratio that, in certain cases, can attain 10’ or the switching time of 1 us (compet- 
itive with flash-type memories), potential switching voltages of 3-4 V, combined with 
predictable fabrication-controlled I-V behavior, simple two-terminal geometry (no 
need for a gate and so dramatically reducing circuitry) and access to mass fabrication 
for the two approaches make them extremely attractive structures for NVM. 

Indeed, thanks to thin-film technologies, electronic functionality can be fore- 
seen in very large quantity and at very low cost on substrates such as plastic and 
paper. Additional functionality is also an attractive feature of carbon-based memo- 
ries. For example, by using these kinds of materials, we open up the route for mem- 
ories on flexible substrates, a key building block to enable the success of flexible 
electronics. Carbon-based resistive memories should also offer the capability for 
multilevel storage and *memristive-like" behavior, as seen in other resistive mem- 
ory materials. Multilevel storage allows storing more than one bit per cell, so in- 
creasing data storage densities, while memristive-like behavior can be exploited to 
provide a remarkable range of signal processing/computing-type operations, in- 
cluding implementing logic, providing synaptic and neuron-like *mimics," and per- 
forming, in a very efficient way, analogue signal processing functions (such as 
multiply accumulate operations). 

In case of graphitic-based memories based on other approaches, the transistor- 
type configuration exploiting the hysteretic effect, in our opinion, has a limited po- 
tential compared to new flash-type memories using graphene or MLG as the FG, 
because their architecture can be unlikely miniaturized and their performances es- 
pecially considering the switching time are difficult to evaluate. In case of flash- 
type memories, these last exploit intrinsical characteristics of the materials such as 
the higher DOS (which allow reducing the voltage for the memory window) or the 
higher work function (which allows reaching more than 10 years of charge storage) 
that really allows improving the performances of this kind of memories. Moreover, 
it has been demonstrated to achieve high-density architecture with limited cross- 
talk. It is necessary to make an effort on the miniaturization of these devices but the 
recent works, for example, [51], make us think that this technology can be quite rap- 
idly in competition with existing flash memories, also because it is CMOS compati- 
ble and easy to implement in existing production lines. 
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2 Plasma surface activation and 
functionalization of carbon-based 
materials 


2.1 Introduction 


Carbon-based materials have been extensively researched for an increasing number of 
applications because of their abundance and low cost, ease of processability, high- 
temperature stability and corrosion resistance, high electrical conductivity and good 
mechanical properties, particularly in the microscopic and nanoscopic scale [1]. The 
uses of carbon-based materials are very diverse, going from mechanical reinforcement 
in polymer matrix composites [2] to electrodes in supercapacitors and fuel cells [1]. 
However, their chemical inertness poses a common challenge in most applications, 
and several surface activation and surface functionalization methods have been devel- 
oped to improve their performance, their compatibility with other materials or their 
interaction with the service environment [3]. 

The chemical inertness of carbon-based materials is attributed to the presence of 
graphene planes on the surface, having a low density of dangling bonds [5]. Graphene, 
carbon nanotubes (CNTs) and buckyballs are extreme examples of this microstructure, 
as they are composed of a single or a small number of graphene layers with a regular 
arrangement [6]. On the other hand, carbon fibers and amorphous carbon exhibit a tur- 
bostratic structure, with a higher degree of disorder and a larger number of dangling 
bonds on the surface (Figure 2.1). Several surface engineering techniques have been 
developed to activate and/or functionalize carbon-based materials [7]. These processes 
aim at increasing the number of dangling bonds and attaching functional groups to 
them to make the carbon-based material more compatible with the polymeric matrix in 
a composite, or to improve the wettability by an electrolyte in a supercapacitor and 
increase the charge exchange rate at the electrode. 

The conventional surface engineering methods include wet chemical processes 
[8, 9], thermochemical [10] or electrochemical treatments [11]. Even though these 
methods are effective, they can be difficult to control and, in some cases, they raise 
environmental concerns. Therefore, recent studies have focused on the development 
of plasma processes for surface activation and functionalization of carbon-based 
materials [12, 13]. Plasma processes are conducted under vacuum or at atmospheric 
pressure and use highly reactive ionized gas molecules to create active sites and 
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Figure 2.1: Turbostratic microstructure of carbon fibers, with inert graphene basal planes on their 
surface. Reprinted with permission from [4]. 


attach functional groups on the surface of carbon-based materials. In general, these 
processes are easier to control than conventional treatments, and they offer techni- 
cal and environmental benefits. 


2.2 Active screen plasma surface engineering 


Active screen plasma (ASP) is a low-temperature and low-pressure plasma technol- 
ogy of surface engineering, which was developed at the University of Birmingham in 
the early 2000s [14-16]. The technology was initially applied to the thermochemical 
treatment of steels, but its potential for the functionalization of other materials was 
recognized immediately [17]. The schematic diagram in Figure 2.2 illustrates the main 
features of the ASP reactor in contrast to a conventional DC plasma reactor. The ben- 
efits of the ASP technology are mainly attributed to the creation of ions and active 
species on an auxiliary surface, the so-called active screen, as opposed to the surface 
of the treated material [18]. This allows more control over the type of plasma-surface 
interactions and the severity of the physical and chemical processes [3]. 

The active screen is a metallic mesh negatively biased to form ions and other 
active species from gas molecules, most typically N2, H5, Ar, CH,. The surfaces to be 
treated can be simply left at floating potential’ and exposed to the flow of active 


1 Floating potential is the electric potential that an insulated object develops when exposed to the 
plasma. An object immersed in the plasma will be exposed to ions and electrons, but the electrons 
move much faster than the ions because of the difference in mass, so the number of electrons arriv- 
ing at the surface of the object is higher than the number of ions. If the object is electrically insu- 
lated, the negative charge of the electrons accumulates, developing a self-induced bias that attracts 
positive ions and repels the less energetic electrons, until an equal number of ions and electrons 
arrive at the surface at any given time, that is, zero net current. 
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Figure 2.2: Schematic diagram of an active screen plasma reactor. Reprinted with permission 
from [19]. 


species, in which case the plasma-surface interactions are mainly of chemical na- 
ture, for example, attachment of functional groups. In other cases, a moderate bias 
is applied on the treated surface in order to attract the active species with increasing 
energy, making the physical action of the ions more pronounced. This ion bombard- 
ment can create new active sites and change the surface roughness. Alternatively, 
the material from the active screen can be sputtered? and deposited on the surfaces 
exposed to the plasma. This approach is particularly useful to functionalize surfaces 
with metallic materials. 

In summary, the ASP technology has an additional degree of freedom compared 
to other plasma reactors, making it possible to create highly reactive ions and excited 
species on the active screen, and attracting them toward the work table with varying 
levels of energy (bias) to control the type of plasma-surface interaction, thus favoring 
physical or chemical phenomena. The following sections illustrate some of the recent 


2 Sputtering refers to the ejection of material from a solid surface by the physical action of impact- 
ing ions. 


20 —— Hanshan Dong, Santiago Corujeira Gallo 


uses of the ASP treatments to activate and functionalize carbon-based materials with 
a wide range of microstructures, from the highly ordered structure of graphene to the 
disordered microstructure of carbon paper. 


2.3 Reduction of graphene oxide for optoelectronic 
devices 


Graphene has attracted significant attention for optoelectronic devices, and graphene 
oxide (GO) is an interesting precursor for graphene coatings of high-quality and over 
large-surface areas [6, 20]. GO has a nanosheet-layered structure similar to graphene 
but contains different oxygen groups [21, 22]. The GO coatings can be produced by 
chemical methods [23], dissolving graphene powder and other salts in sulfuric acid, 
subsequently dispersing and exfoliating the GO nanosheets by sonication, and finally 
casting it onto the substrate. A reduction process is then required to remove the oxy- 
gen functionalities and obtain a graphene coating of high quality, and the ASP tech- 
nology was used for this purpose [24]. The ASP treatments were conducted at 400 Pa 
in pure H; or N7- H; mixtures, at temperatures between 100 and 200 °C for 1 h. 

The appearance of coatings changed from a light brown color to silver-gray after 
the plasma reduction in pure hydrogen and nitrogen-hydrogen mixtures, and this 
change was accompanied by a decrease in optical transparency and electrical resis- 
tance of the coating (Figure 2.3). The reduction process with the plasma technology 
was more effective than the heat treatment of GO in a reducing atmosphere, and the 
results varied with the treatment temperature (100-200 °C) and the gas mixture (H3 
or H;-N;). The X-ray diffraction (XRD) traces revealed that the GO peak disappeared 
after the reduction process, and the Raman spectra indicated an increase in the ratio 
of disordered to graphitic bands (D/G), which was attributed to the introduction of 
defects in the graphene planes during the reduction process. The X-ray photoelectron 
spectroscopy (XPS) data showed a marked decrease in oxygen functional groups at- 
tached to the carbon network, with the carboxylic oxygen (O—C=O) being the most 
unstable functional group compared to other oxygen bonds. The XPS data also re- 
vealed the presence of Fe, Cr and Mo doping elements in a partially oxidized state. 
These metallic dopants are transferred from the mesh of the active screen equipment 
to the GO-coated substrates, particularly when using nitrogen containing gas mix- 
tures. A small amount of nitrogen was also incorporated into the coatings. 

The oxygen to carbon ratio (O/C), reported in Table 2.1, reflects the efficiency of 
each process to reduce the GO coating. The plasma treatment in pure hydrogen was 
the most effective, resulting in an extensive decomposition of carboxyl functionali- 
ties and the creation of carbonyl groups. On the other hand, the use of nitrogen in 
the gas mixture increased the doping of the GO coating with metal atoms from the 
active screen, with a marked drop in electric resistance. The level of doping can be 
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Figure 2.3: Properties of the untreated and reduced GO coatings: (a) optical transparency and 
(b) electrical resistance. Reprinted with permission from [24]. 

PET: poly(ethylene terephthalate) substrate; GO: graphene oxide coating; ASP: active screen 
plasma surface treatment. 


controlled by varying the processing conditions (composition of the active screen, 
gas mixture, electric bias applied on the treated surface and duration of the pro- 
cess), opening new pathways to reduce and functionalize GO coatings. 
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Table 2.1: Chemical analysis from XPS: O components and O/C ratio. Reprinted with permission 
from [24]. 


Sample Treatment Gas Temperature O-Si  O—Me  O—C-O  O-C 0/C 


GO Untreated N/A N/A 2.6 0.0 9.2 9.0 0.235 
HT200 Heat treated H2 200 °C 4.3 0.0 7.4 5.7 0.154 
ASP150 ASP H2 150 °C 5.6 6.7 6.9 1.5 0.122 
ASPN150 = ASP H2-N2 150°C 2.1 7.4 7.4 1.3 0.125 
ASP200 ASP H5 200 °C 1.7 6.2 6.2 1.1 0.102 


ASPN200 ASP H5-N5 200 °C 2.4 6.8 6.8 1.1 0.115 


2.4 Functionalization of carbon nanofibers 
for supercapacitors 


Carbon nanofibers (CNF) are good candidate materials for electrodes in energy stor- 
age devices because of the intrinsic properties of carbon (thermal and chemical sta- 
bility, electric conductivity) and the ability to form controlled porous structures with 
the nanofibers [1]. This combination allows researchers to create stable electrodes 
with a very large surface area, and recent efforts have focused on increasing the sur- 
face conductivity and the charge exchange efficiency with the electrolyte. Silver and 
nickel are among the proposed electroactive materials, capable of increasing the 
pseudocapacitance and double-layer energy storage mechanisms,’ and the ASP 
technology was used for this purpose [25]. CNFs were dispersed in isopropanol, the 
suspension was placed in Petri dishes and the solvent was left to evaporate before 
the plasma treatments. 

The ASP treatments did not change the size or the structure of the CNFs, but their 
surface was covered with nickel and silver nanoparticles, respectively. The metal 
nanoparticles were approximately 10-50 nm in diameter (Figure 2.4), and their size 
and density increased with the ASP treatment time, which was reflected on the nickel 
and silver contents detected by energy-dispersive X-ray spectroscopy and XRD. The 
ASP functionalization with silver produced a rough surface, even after 120 min of 
plasma treatment, which is beneficial to increase the surface area. On the other hand, 
nickel tended to form a continuous coating on the CNFs after only 45 min of ASP treat- 
ment, thus limiting the amount of material and the morphology of the particles that 
could be deposited on CNFs. 


3 Electrical double layer capacitors are based on non-Faradaic mechanisms, in which the electric 
charge is stored in the double layer formed at the interface between the inert electrode material and 
the electrolyte, without involving redox reactions. On the other hand, the pseudocapacitors rely on 
reversible redox reactions (Faradaic), which take place on the surface of the active electrode mate- 
rial in contact with the electrolyte [1]. 
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Figure 2.4: Surface morphology of carbon nanofibers: (a) untreated and (b) ASP functionalized with 
nickel. Reprinted with permission from [25] Copyright (2017) American Chemical Society. 


The cyclic voltammetry (CV) tests conducted in 1 M sodium sulfate solution re- 
vealed a significant increase in the loop area after the ASP treatments, particularly 
with silver (Figure 2.5). More importantly, the CNFs functionalized with silver retained 
the rectangular shape of the CV curves, which was attributed to the high surface con- 
ductivity and ion transport efficiency. On the other hand, the nickel-functionalized 
CNFs showed severe distortions, particularly at higher scan rates, between 100 and 
200 mV/s. The differences were also observed in the cyclic charge-discharge curves, 
with the nickel-functionalized CNF electrodes showing a twofold increase in capaci- 
tance and the silver-functionalized CNF electrodes exhibiting a fourfold increase 
in capacitance compared to the untreated CNFs. The electrochemical impedance spec- 
tra indicated that the increase in capacitance was mainly due to the double layer 
mechanism, rather than pseudocapacitance. However, recent reports have indicated 
that vanadium and other transition metal nitrides can improve both Faradaic and 
non-Faradaic charge storage mechanisms, which make them promising electrode ma- 
terials for supercapacitors. This finding opens good prospects for the ASP technology, 
considering that it has been extensively used to form nitrides and carbides for surface 
functionalization and hardening. 


2.5 Surface activation of carbon paper gas diffusion 
layers in fuel cells 


Carbon paper is used as gas diffusion layer (GDL) in direct methanol fuel cells 
[26]. The electrocatalysts, which consist of reaction layers of precious metals with 
nanometer-size features, are frequently deposited directly on the GDL. These depo- 
sition processes are typically conducted in an aqueous environment and the poor 
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Figure 2.5: Cyclic voltammetry of carbon nanofibers functionalized with silver (a) and capacitance 
(b) calculated from charge-discharge experiments. Reprinted with permission from [25] Copyright 
(2017) American Chemical Society. 


wettability of the hydrophobic carbon paper results in aggregation and uneven 
growth of the electrocatalyst, which reduces efficiency and increases cost [27, 28]. 
Therefore, the ASP technology was used to activate the surface of SIGRACET® 
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GDL 35BC carbon paper [19]. The activation treatment was conducted in a gas mix- 
ture of 25% N; and 7596 H; at a pressure of 400 Pa. The temperature ranged from 
100 to 210 ?C and the duration of the treatment was between 10 and 30 min. After 
the activation treatment, the carbon papers were rinsed in ethanol and immersed 
in an aqueous solution of H;PtClg and formic acid to grow the platinum nanowires 
(2 mg/cm’), and finally rinsed in water/ethanol/water and dried at 40 °C for 24 h. 
These carbon papers were then used as integrated cathodes in single cell tests, with a 
Nafion 117 membrane and the commercial JM anode (PtRu 4 mg/cm?). The electrodes 
were then hot-pressed against the Nafion 117 membrane at 135 ?C for 2 min under a 
constant pressure of 5 MPa. The reference cells were fabricated with Nafion 117 mem- 
brane, commercial JM cathode (Pt 4 mg/cm?) and JM anode (PtRu 4 mg/cm’). These 
electrode assemblies were sandwiched between two graphite flow field plates to form 
a single cell with an active electrode area of 5 cm". 

The CV tests conducted on a single cell revealed that all the carbon papers acti- 
vated with the ASP technology for 10 min had higher electrochemical surface area 
(ECSA) than the untreated carbon paper. What is more, the ECSA increased with the 
ASP treatment temperature up to 150 ?C, and then decreased with further tempera- 
ture increments (Figure 2.6). The maximum ECSA value of 25.8 m?/g corresponded 
to the carbon paper activated at 150 ?C with the ASP technology, which was 6596 
higher than the ECSA value of the untreated carbon paper. The polarization curves 
also revealed maximum power densities for moderate ASP treatment temperatures 
in the order of 120 and 150 ?C. The highest power density measured in these experi- 
ments was 63.9 mW/cm?, which was three times higher than the power density ob- 
tained with the untreated carbon paper. 

This remarkable improvement in the performance of the single cells was attrib- 
uted to the different surface morphology observed at high magnification. The sur- 
face of the untreated carbon paper was only partially covered by Pt nanowires, and 
they formed rather large cubic-shaped clusters. The uneven surface morphology 
was a result of the super hydrophobic nature of the untreated carbon paper and its 
poor wettability (Figure 2.7). In contrast, the surface of the ASP-activated carbon 
paper showed a more uniform distribution of the Pt nanowires, with a decreasing 
number of aggregates as the ASP treatment temperature increased up to 150 °C. At 
temperatures higher than 150 °C, the morphology of the Pt nanowires changed, be- 
coming shorter and, consequently, reducing the ECSA and the cell performance. 
Longer ASP activation times, that is, 30 min, had a detrimental effect on the surface 
morphology of the catalyst layer and, consequently, on the cell performance. In 
summary, the ASP activation at moderate temperature (120—150 ?C) and for a short 
time (10 min) can effectively promote the growth of the electroactive layer of Pt 
nanowires on the carbon paper, improving the utilization of this expensive catalyst 
and, ultimately, increasing the ECSA and the power density of the fuel cell. 
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Figure 2.6: (a) Electrochemical surface area (ECSA) and (b) power density from polarization tests. 
Reprinted with permission from [19]. 


2.6 Improving the bond strength of carbon 
fibers in composite materials 


Carbon fibers are exceptional reinforcements for polymer matrix composite materi- 
als, particularly because of their high strength, high stiffness and lightweight [29]. 
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Figure 2.7: Surface morphology of the electrocatalyst layer deposited on the (a) untreated and 
(b) ASP functionalized carbon paper (ASP treatment conducted at 120 ?C for 10 min) [19]. 


However, the graphitic basal planes present on the surface of the carbon fibers make 
them relatively inert and this reduces the wettability and the bond strength with the 
polymeric matrix, frequently an epoxy resin [3]. Plasma treatments have been used to 
introduce polar functional groups on the surface of carbon fibers. These functional 
groups help to modify the bond strength with the polymeric matrix and to tailor the 
mechanical properties of the composite material, particularly the strength and tough- 
ness [30, 31]. The ASP technology offers good control over the treatment conditions to 
avoid damaging the fibers, thus retaining their strength. 

Carbon fibers were ASP treated in a gas mixture of 75% H, and 25% N; at a pres- 
sure of 50 Pa. This atmosphere was selected to introduce -NH, functional groups 
on the surface of the carbon fibers [32]. These functional groups are known to create 
strong chemical bonds with the epoxy resins used in these composites, and the 
physical etching creates a microroughness on the carbon fibers, which increases 
the mechanical interlocking with the matrix. Optical emission spectroscopy studies 
conducted during the ASP treatments revealed NH, species in the plasma phase, 
and NH groups were also identified on the surface of the functionalized carbon fi- 
bers using XPS (Figure 2.8). The relative intensity of the NH species in the plasma 
decreased with increasing applied power, and this change correlated with a reduc- 
tion in the nitrogen content on the surface of the functionalized fibers. A similar 
trend was observed on the contact angle of glassy carbon disks, functionalized 
alongside with the carbon fibers. All the ASP-treated carbon disks showed contact 
angles significantly lower (<60°) than the untreated disks (80°), but the lowest con- 
tact angle (40?) corresponded to the disk treated at low power and for short time. 

Further wettability measurements were conducted on carbon disks immedi- 
ately after the ASP treatment and over the subsequent 4 weeks. It was found that 
the contact angle dropped remarkably after the plasma treatments to values as low 
as 10°. However, the contact angle increased rapidly over the subsequent 7 days 


28 —— Hanshan Dong, Santiago Corujeira Gallo 


Nitrogen content from XPS [at%] 
OES relative intensity [cts] 


200V 300V 400V UT 
Potential/treatment condition 


Figure 2.8: Correlation of active species in the plasma phase with functional groups on the surface 
of the carbon fibers as a function of the ASP treatment conditions (potential). Reprinted with 
permission from [30]. Copyright 2017, American Vacuum Society. 


and remained relatively constant for the remaining 3 weeks (Figure 2.9). The chemi- 
cal analysis by XPS showed a small but consistent reduction in the nitrogen content 
(from 5.8 to 5.1 at.96) over a period of 28 days. However, the initial XPS analysis was 
only conducted 3-4 weeks after the ASP treatment, so the wettability measurements 
could not be correlated with the XPS data. The change in the nitrogen content and in 
the wettability was consistent with an exponential decay, which could be attributed 
to a degradation of the nitrogen functionalities by oxidative scissions and evaporative 
loss upon their exposure to air. The elucidation of the degradation mechanism re- 
quires further study. 


b) c) 
Untreated disk AS plasma treated AS (aged 2 weeks) AS (aged 4 weeks) 


Figure 2.9: Evolution of the contact angle of an ASP-treated carbon surface with time. Reprinted 
with permission from [30]. Copyright 2017, American Vacuum Society. 


2.7 Synthesis and deposition of carbon 
nanotube films 
CNTs exhibit extraordinary electrical, thermal and mechanical properties, similar to 


CNFs. CNTs are frequently synthesized by arc discharges, laser ablation or chemical 
vapor deposition (CVD) [33-35]. The latter method is capable of producing vertically 
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aligned CNT (VACNT), but the process requires prior deposition of metal nanopar- 
ticles on the substrate to act as catalysts and nucleation sites. In addition, the high 
temperature required for conventional CVD processes can be detrimental for the 
substrate. It is known that plasma-enhanced CVD (PE-CVD) can reduce the deposi- 
tion temperature, but the energy of the ions needs to be carefully controlled to en- 
hance rather than suppress the growth of CNT [36]. A multistep process to deposit 
metal nanoparticles and grow CNTs at low temperatures has been recently devel- 
oped based on the ASP technology [33]. 

Silicon wafers (100) were used as substrates for the PE-CVD process with RF excita- 
tion (13.56 MHz, 30 W) and a modified active screen experimental setup (Figure 2.10). 
Metal nanoparticles with different compositions were deposited on the silicon wafers 
by changing material of the active screen, in this case nickel and stainless steel. The 
size, morphology and uniformity of the metal particles were critical for the success of 
the process, the density and vertical alignment of the CNTs. VACNTs were grown in the 
same active screen experimental arrangement using a gas mixture of acetylene, argon 
and hydrogen at 100 Pa (Figure 2.11). Transmission electron microscopic observations 
revealed a tip-growth mode, in which the metal nanoparticles were detached from the 
silicon wafer and remained at the tip (growing front) of the CNTs. For a given treatment 
time (100 min), the length of the CNTs increased with the processing temperature from 
0.6 um at 400 °C to 6 um at 500 °C, and their diameter decreased from 28 to 16 nm, 
respectively (Figure 2.11). The structure of the CNTs did not change significantly with 
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Figure 2.10: Schematic diagram of the active screen setup used for PE-CVD of vertically aligned 
carbon nanotubes. Reprinted with permission from [33]. 
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Figure 2.11: SEM images of CNTs grown at (a) 400 °C and (b) 500 °C. Reprinted with permission 
from [33]. 


the deposition temperature, but the ratio of D and G bands (Ip/Ic) in the Raman spectra 
changed from 0.7 to 0.97 with increasing content of acetylene in the gas mixture. 

The multistep PE-CVD process, based on the ASP technology with RF excitation, 
could successfully deposit a layer of metal nanoparticles on a silicon wafer and subse- 
quently grow vertically aligned CNTs. The size, morphology and microstructure of the 
VACNTS could be controlled by changing key processing parameters, such as size, dis- 
tribution and uniformity of metal nanoparticles, gas mixture and processing tempera- 
ture. This process represents a simple and cost-effective method to functionalize the 
surface of a wide range of substrates with CNTs. 


2.8 Summary and future perspectives 


This chapter presents a summary of the research undertaken at the University of 
Birmingham on the functionalization of carbon-based materials using the ASP technol- 
ogy. This research exploits the typical physical and chemical phenomena associated 
with the plasma technology, further optimized by the additional control and degrees of 
freedom offered by the ASP technology. In this way, the highly reactive species in the 
plasma were used to reduce GO films, activate the surface of carbon paper membranes 
or functionalize the surface of carbon fibers. In addition, the sputtering of material 
from the active screen and its deposition on the surfaces exposed to the plasma was 
used to increase the conductivity and surface area of CNFs and to create nucleation 
sites for CNTs. In all cases, the treatment conditions could be optimized, targeting spe- 
cific plasma-surface interactions, to achieve the desired properties. 

Future efforts will focus on expanding the range of gases and active screen materi- 
als used to functionalize temperature-sensitive substrates, as well as combining the 
ASP technology with pre- and post-treatments. There are also challenges related to the 
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uniformity of the ASP treatment of porous and three-dimensional structures due to 
shadowing effects that need to be addressed. In addition, the interaction between ac- 
tive species in the plasma and different microstructures of carbon-based materials de- 
serves further study, together with the evolution of the ASP functionalized surfaces 
over time. Lastly, there seem to be good opportunities to create patterned surfaces 
using the ASP technology, in combination with other processes with potential applica- 
tion in multiple fields such as functional materials, sensors and biosensors, electrical, 
electronic and optoelectronic devices. Without question, the ASP technology remains 
to be a valuable enabling technology for the surface functionalization of engineering 
materials, showing broader prospects of practical applications than ever before. 
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3 Modeling carbon-based smart materials 


3.1 Introduction 


Smart materials show changes in their properties in a predictable manner upon ap- 
plication of external stimuli, such as chemical, electrical, thermal, mechanical, 
magnetic and light ones. Our knowledge about smart materials exists from the era 
of Roman Empire: there are several evidences that building materials used at that 
time, mainly mixtures of several compounds including limestone, had self-healing 
capabilities, since they were able to repair cracks by the simple action of rainwater 
[1]. Modern smart materials are opening up new opportunities in several fields, for 
instance, healthcare, defense, waste management and packaging. In fact, these ma- 
terials are widely used in applications such as sensors, actuators, robots, artificial 
muscles and drug delivery. Smart materials can be classified based on their re- 
sponses to the external stimuli, that is, shape memory, piezoelectric, photorespon- 
sive, electroresponsive, magnetoresponsive and thermochromic materials. 

Shape memory materials have the unique capability of remembering their origi- 
nal shape after deformation, thereby returning to it upon application of external 
stimuli such as thermal, mechanical or magnetic variations [2, 3]. These materials 
have found application in automobiles, aircraft, biomedical devices, robots, civil 
structures and textiles [3-8]. Shape memory materials can be made out of alloys, 
ceramics, polymers, hybrid materials or gels. The first shape memory material dates 
back to 1932, when Arne Olander first discovered a solid cadmium-gold alloy that, 
when deformed in a cold state, returned to its original shape upon heating [9]. 
Nickel-titanium (NiTi) was then a popular shape memory alloy (SMA) discovered 
by Buehler et al. [10] in 1963, which spread the use of SMA to the first industrial 
applications. NiTi alloys are preferable in most applications due to their good duc- 
tility, corrosion resistance, mechanical stability and biocompatibility [11, 12]. In re- 
cent years, carbon-based shape memory polymers attracted attention due to their 
capability of handling high strain and temperature, their lower density and recov- 
ery temperature and their flexibility and recyclability [13-15]. All these superior 
properties make them interesting shape memory materials, particularly for aero- 
space and biomedical fields. 
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For example, Zhang et al. [16] fabricated multistimuli-responsive shape memory 
materials based on polymer composites reinforced by graphene nanoplatelets (GNPs), 
which could be triggered by electrical and infrared stimuli. Figure 3.1 shows the shape 
memory mechanism actuated by photothermal effect: after a deformation at low tem- 
perature (Shape 2), the original shape of the GNP-based polymeric material (Shape 1) 
can be recovered upon exposure to infrared radiation [16]. Recently, carbon-based 
shape memory polymer nanocomposites seem promising candidates for electrical ac- 
tuators. In fact, polymer nanocomposites filled with carbon nanotubes (CNTs) have 
shown low electrical resistivity while fast thermal response to applied voltage. These 
composites can be tailored to specific actuation requests, since a desired shape recov- 
ery can be attained by printing CNTs with a specified pattern throughout the poly- 
meric matrix [17]. 
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Figure 3.1: Schematic illustration of infrared radiation-actuated shape memory effect of GNP-based 
polymer blend. Reproduced from Ref. [16] with permission from the Royal Society of Chemistry. 


Piezoelectric materials are capable of converting mechanical deformation into elec- 
tricity, and vice versa [18, 19]. These materials are widely used as sensors, actuators, 
structural health detectors and energy generators [19-22]. Particularly, there is a 
growing demand of wearable sensors for monitoring physiological and biomechani- 
cal signals in sports and medical applications. However, bulky batteries are cur- 
rently needed to power these sensors; therefore, smart garments generating their 
own power by human body motion would be highly beneficial (see Figure 3.2). This 
requirement can be, for instance, fulfilled by polymeric fiber-based piezoelectric 
generators, which are flexible and durable as well [23-25]. 
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Figure 3.2: Sketch showing power generation from human movement using fiber-based piezoelectric 
generator. Reproduced from Refs. [23, 26] with permission from the Royal Society of Chemistry. 


Zhong et al. [24] fabricated a piezoelectric fiber using cotton threads coated with 
CNTs and dipped into an aqueous suspension of polytetrafluoroethylene. This carbon- 
based polymeric piezoelectric fiber utilizes the biomechanical energy from human mo- 
tion to generate electrical energy, with an average power output of ~0.1 uW/cm?, which 
would be enough to energize wearable sensors or smart shirts (see Figure 3.3). 

Photoresponsive materials show changes in their physical and chemical behav- 
ior when exposed to light. Usages of these smart materials include anticounterfeit- 
ing applications [27], chemical sensors [28], optical switches [29], flexible rewritable 
optical memories [30] and toys (e.g., doll whose skin get tanned under the sun [31]). 
The most widely used photoactive agents are azobenzene, spiropyrans, spirooxa- 
zines, diarylethenes and fulgides, which are summarily depicted in Figure 3.4. 

Recently, Renuka et al. [29] developed an optical switch using multiwalled 
CNTs (MWCNTs) and azo-based chromophores (Figure 3.5). The electrical conduc- 
tion of the composite is controlled by azo-based system, which shrinks when ex- 
posed to UV irradiation (cis isomer), thereby increasing the electrical conduction, 
because of the reduced tunneling barriers through the MWCNTs. The conduction 
drops down again when UV irradiation stops, because azo molecules return to their 
trans isomer configuration. 
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Figure 3.3: Schematic diagram of a smart shirt. Reprinted from Ref. [24] with permission from 
American Chemical Society. 


Lately, electroactive polymers (EAPs) are gaining interest because of their unique 
property to change geometrical/structural properties in response to the application 
of an electric field. Due to this capability, EAPs can be employed in artificial 
muscles, sensors, robots, biomedical devices and energy storage systems [33, 34]. 
Magnetoresponsive materials, instead, can be shaped upon the application of a mag- 
netic field. Magnetoresponsive fluids, for instance, can be used to construct dampers 
with adaptive response and thus enhanced suppression of vibrations. These can be 
fitted to buildings and bridges to mitigate the detrimental effect of strong winds or 
earthquakes. Among other applications, magnetoresponsive materials can be used 
also in nanomedicine [35, 36]. Finally, thermochromic materials change color in re- 
sponse to different temperatures. They have been used, for example, in bath plugs 
that change their color when water is too hot. 

Carbon-based polymer nanocomposites seem to be among the most promising 
candidates for the development of smart materials. However, significant efforts are 
still needed to design smart materials with specific thermal, mechanical and electri- 
cal properties. Since these properties derive by mechanisms coupled from nano- to 
macroscale, multiphysics simulations could support a more rational research of 
smart materials with tailored properties. For example, the overall piezoelectric re- 
sponse of carbon-based polymeric smart materials is influenced by several mecha- 
nisms, including tunneling of electrons and formation of percolating networks of 
carbon fillers in the polymer composite [37-39], and should therefore be addressed 
by linking atomistic to mesoscopic models. In the following sections, the atomistic, 


3 Modeling carbon-based smart materials —— 


N, 
UV light Cy N 


Visible 


cis-Azobenzene 


UV light 
Visible 


Nt 
| 


Spiropyrans, (Open) 


Oo 


UV light 
Visible 


UV light 


Visible 
Diarylethenes, (Open) Diarylethenes, (Closed) 
[9] o 
" UV light 
| o Visible " à 
O oO 
Fulgides, (Open) Fulgides, (Closed) 


Figure 3.4: Chemical structure of photoresponsive molecules [32]. 
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Figure 3.5: Schematic diagram of an optical switch based on MWCNT/azo composite. Reprinted 
from Ref. [29] with permission from American Chemical Society. 
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mesoscopic and macroscopic simulation tools that allow to explore the main properties 
of nanocomposites are extensively presented, and future perspectives of advanced 
modeling of carbon-based smart materials are discussed. 


3.2 Multiphysics modeling 


A better understanding of the thermophysical behavior of carbon-based polymer 
nanocomposites at different length- and timescales could ease the discovery of novel 
smart materials, and could be facilitated by proper material modeling approaches 
[40-47]. Advanced modeling approaches to the properties of nanocomposites are 
based on linking and/or coupling various model types that are applicable at different 
length- and timescales (see Figure 3.6). 
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Figure 3.6: Schematic diagram of advanced modeling of polymer nanocomposite across 
multiple scales. 


This situation is ubiquitous in many sectors, where modeling can support the effec- 
tive development of advanced materials, but the heterogeneity of modeling techni- 
ques may prevent to highlight the computational synergies. For this reason, the 
European Commission published a Review of Materials Modelling (RoMM), now in 
its sixth edition [48], which provides a classification of materials modeling that en- 
ables a coherent description of model types and a standardized documentation of 
modeling data (called *MODA") of materials. Furthermore, the RoMM contains sev- 
eral examples of MODA documentation applied to a compendium of applications 
from EU H2020 LEIT NMBP Materials projects. Based on the above review, the 
European Materials Modelling Council (EMMC, http://www.emmc.info) proposed a 
CEN Workshop Agreement about “Materials modelling — terminology, classification 
and metadata" [49], endorsed by more than 15 European organizations with the ob- 
jective of standardization of terminology, classification and documentation of materials 
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modeling and simulation. Moreover, the EMMC already proposed a European Materials 
Modelling Ontology (EMMO) [50]. The EMMO aims at addressing the granularity levels 
of materials modeling (atomistic, electronic, mesoscopic and continuum), which do 
not really depend on the scale they are applied to. This is why multiscale modeling is 
somehow confusing and it would be better to talk about multiphysics modeling. 

In particular, at atomistic/molecular level, molecular dynamics (MD) simula- 
tions allow understanding the behavior of the constitutive elements of nanocompo- 
sites (such as carbon fillers and polymer matrices) and, most importantly, their 
interaction at the interface [51]. The structural and interfacial characteristics of the 
constitutive elements of nanocomposites at molecular level are successively em- 
ployed in coarser modeling descriptions (mesoscopic models). The mesoscopic 
structure of a smart composite is modeled by a representative volume element 
(RVE) of the material, and the properties computed at mesoscale can be finally 
homogenized to evaluate the effective thermophysical properties at macroscopic 
(continuum) level. 


3.2.1 Atomistic models 


Thermophysical properties of carbon-based nanocomposite materials are mainly 
governed by the interaction between nanofillers and surrounding matrix at molecu- 
lar level. MD simulations are typically the most adequate tool to study such atomis- 
tic phenomena. In this section, the process of building the molecular geometry of 
nanocomposite is first discussed; then, the MD protocols for assessing thermal, me- 
chanical and electrical properties are described, and main results from the litera- 
ture are presented. 


Building the atomistic setup 


While the atomistic structure of carbon nanofillers (e.g., nanotubes, nanoplatelets 
and nanoribbons) can be conveniently generated by tools available in the literature 
[52-54], one of the major challenges in simulating polymer-based nanocomposites at 
the molecular level is to accurately describe the structure of polymeric matrix. Curing 
is a chemical reaction that happens between the components of thermosetting plas- 
tics (e.g., epoxy resins), namely the resin and hardener, which result in a final cross- 
linked (cured) structure. MD methods possess an inherent advantage over continuum 
approaches, since they are able to accurately simulate the mechanisms underlying 
curing reaction. Schulz and Frisch [55] first attempted to model curing procedure at 
the molecular level. They employed a lattice Monte Carlo simulation method to un- 
derstand the reaction kinetics of curing process in terms of degree of polymerization 
and molecular weight distribution. Later works also included topological information 
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of polymer networks in the description of curing process, therefore allowing the sim- 
ulation of thermomechanical properties of epoxy resins by MD simulations [56-64]. 

Protocols for simulating the curing process and creating the atomistic setup 
of epoxy resins can be classified into two groups: (1) single-step methods, where 
all the potential chemical bonds between reactive sites within a specified cut-off 
distance are formed at once; (2) multistep approaches, where reactive pairs within 
a spatial distance are cross-linked together and the resulting structure relaxed by 
MD, iteratively, until the desired curing degree is met. Generally, single-step ap- 
proaches lead to large artificial strains in the cured structure; whereas, multistep 
methods — even though they are more computationally expensive — provide struc- 
tures closer to real curing reactions, since they operate through intermediate MD 
relaxations. 

Jang et al. [64] compared the performance of these two approaches on a system 
consisting of bisphenol A diglycidyl ether (DGEBA) epoxy and poly(oxypropylene) 
diamine (POP) curing agent. In the single-step method, a Monte Carlo algorithm 
was adopted to allocate cross-linking bonds. For the multistep algorithm, an initial 
cut-off radius of 10 A was assigned, and covalent bonds between all pairs of poten- 
tially reactive atoms within the cut-off radius generated. After that, the partially 
cured system was relaxed in the NPT and NVE ensembles, respectively. In addition, 
the partial charges of atoms were adjusted to neutralize the net charge in the sys- 
tem. Finally, in the case that no reactive pairs remained within the current cut-off 
radius, the reaction cut-off radius was increased iteratively. The thermal (thermal 
expansion, heat capacity and glass transition temperature), structural (dihedral 
angle distributions, radial distribution functions, minimum inter-nitrogen contour 
length distributions and fragment molecular weight distribution) and mechanical 
(Young’s modulus) properties of the systems cured by either single- or multistep ap- 
proach were then compared. Results showed that the fragment molecular weight 
distribution was different in the two cases, since fewer and larger fragments were 
obtained with the single-step approach; nevertheless, the other properties had no 
significant dependency on the choice of cross-linking method. 


Thermal properties 


Thermal conductivity 

The thermal conductivity (A) of nanostructured materials can be estimated by either 
equilibrium MD (EMD) or nonequilibrium MD (NEMD) simulations. It is noteworthy 
to mention that EMD simulations are best to study homogeneous systems, whereas 
NEMD techniques are more suitable for heterogeneous systems [65]. Green-Kubo 
formulation underlies EMD simulations, since it relates A with the fluctuations of 
thermal current via the fluctuation-dissipation theorem [66, 67]: 
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where kg is the Boltzmann constant, T and V are the temperature and volume of the 
system, respectively. In addition, J(t) is the heat current vector at time t and 
(J(0)(0)) is the heat current autocorrelation function. In terms of MD entities, the 
heat current vector is commonly defined as 


d N 
J(t) - 5: En (3.2) 
dt 


being r; and Ej the position and total energy of atom i, respectively. E; can be com- 
puted by summing up kinetic and potential energies of the atom, namely 


Ej- simi? s ly Uu) ri), (3.3) 


ir 


where m; and v; indicate mass and velocity of atom i, respectively. U,(rj) is the 
total potential energy of atom i, which depends on the type of bonded and non- 
bonded interaction potentials used in the simulations and rj is the distance be- 
tween atoms i and j. 

For example, Kumar et al. [68] used EMD to study the correlation between A of 
cross-linked epoxies (DGEBA resin, 4,4’-diaminodiphenyl sulfone - DDS hardener) 
and temperature. Results demonstrate that the inclusion of long-range Coulomb in- 
teraction corrections leads to a better agreement between modeling and experimental 
results. Sirk et al. [69], instead, employed MD simulations to study the effect of epoxy 
composition on its thermal, structural and volumetric properties. Different cross- 
linked networks including DGEBA resin and various mixtures of flexible (POP) and 
stiff (MCA, 4,4’-Methylenebis(cyclohexylamine)) curing agents were built. They 
found out that the flexibility of cross-linker has considerable impact on the thermal 
and volumetric characteristics of the cross-linked network, especially at temperatures 
close to glass transition. Fasanella and Sundararaghavan [70] adopted the EMD ap- 
proach to study the effect of temperature and functionalization on A of DGEBA + DDS 
epoxy, possibly reinforced by single-walled CNTs (SWCNTS). In particular, the ther- 
mal conductivity was computed both perpendicularly and along the main axis of the 
nanotube. Results showed that, for what concerns the direction perpendicular to the 
CNT axis, the A of pristine SWCNT/epoxy nanocomposites is lower than the neat 
epoxy one, whereas an improved thermal performance is shown for functionalized 
SWCNT/epoxy nanocomposites (see Figure 3.7). In addition, a significant enhance- 
ment in the A along nanotube axis has been observed in the case of pristine SWCNT/ 
epoxy in comparison to both neat epoxy and functionalized SWCNT/epoxy nanocom- 
posite (see the inset in Figure 3.7). 
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Figure 3.7: Average thermal conductivity in the temperature range of 220-420 K in the 
perpendicular and axial directions with respect to the nanotube axis for different simulated 
nanocomposites. Adapted with permission from Springer Nature [70]. 


NEMD methods for computing thermal conductivity of nanostructured materi- 
als are based on Fourier’s law. For example, the thermal conductivity of an isotropic 
system with temperature gradient along z direction can be found as [71] 


dz 


Mee dT /dz’ 


(3.4) 
where q,=Q,/A is the specific heat flux, namely, the heat flux Q, transferred 
through the surface A along the direction z. In addition, dT/dz is the steady-state 
temperature gradient along z-axis. The temperature profile throughout the simula- 
tion domain can be computed after dividing the box into N slabs (depending on the 
system size) along the desired direction. Equation (3.4) can be exploited in nonequi- 
librium simulations to compute A according to two possible approaches. 

The first approach is typically called NEMD (see an example in Figure 3.8), where 
a temperature gradient is imposed across the extremities of computational system. 
This is achieved by a proper design of regions with either high (red atoms in 
Figure 3.8) or low (blue atoms in Figure 3.8) controlled temperature. To sustain the 
temperature gradient during the simulation, energy is constantly supplied and ex- 
tracted from the high- and the low-temperature regions by thermostats, respectively. 
Such addition or extraction of energy in the hot or cold regions is generally carried 
out by modifying the kinetic energy of atoms through a velocity rescaling procedure. 
The instantaneous local kinetic temperature (T;) in each slab can be computed as 
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where N; is the number of atoms in the ith slab and m; the mass of the kth atom. 
The energy exchange between hot and cold regions at steady state causes a specific 
thermal flux (q;) through the system, which can be calculated as 


Q; 1 1n% p 2 
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where N; is the number of atoms in the thermostated region. Besides, v; and v’x rep- 
resent the atomic velocities in the thermosetted regions (hot or cold regions) before 
and after rescaling, respectively. Note that, in the considered geometry (Figure 3.8), 
the heat flux in eq. (3.6) is divided by 2 due to the symmetrical simulation box. 

The second approach is named as reverse NEMD (RNEMD), and it was originally 
conceived by Müller-Plathe [72]. The main difference between RNEMD and NEMD 
methods is that the former generates nonequilibrium conditions through the compu- 
tational domain by fixing heat flux (which gives rise to temperature gradient), while 
the latter by fixing the temperature gradient (which gives rise to heat flux). In 
case of RNEMD, the heat flux is induced by interchanging the atomic velocities 
between hot and cold regions. In the hot region, the coldest atom is selected and 
its velocity (Veola, red region in Figure 3.8) is exchanged with the velocity of hottest 


Figure 3.8: Specification of hot and cold baths for computing the thermal conductivity of 
a carbon-based nanocomposite by NEMD simulations. 
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atom in the cold region (vy, blue regions in Figure 3.8). This process eventually 
leads to an artificial specific heat flux (q,) flowing from cold to hot regions. All velocity 
interchange procedure (i.e., transfers) follows the total energy conservation law; there- 
fore, at steady state, equal amount of energy per area and time is induced from hot to 
cold regions via heat conduction, namely 


1 m 
dz Sant 5 5 vna IVcoial). (3.7) 


transfers 


The most significant advantage of RNEMD over classical NEMD is the fact that heat 
flux, which is a slowly converging quantity with respect to temperature, is fixed 
from the beginning of simulation. Therefore, in general, RNEMD approaches pro- 
vide faster convergence of À computation with respect to classical NEMD ones. 

Interestingly, Varshney et al. [73] evaluated the thermal conductivity of a cross- 
linked network of EPON-862 and curing agent W (diethyltoluenediamine, DETDA) 
using both EMD and NEMD approaches. In the case of EMD, a cubic simulation box 
(53.6 x 53.6x 53.6 À?) was employed; in the case of NEMD, a thin slab with large as- 
pect ratio along the heat flow direction (21.4 x 21.4 x 373.3 A?) was considered instead. 
The results obtained by both methods were in good agreement with experiments, 
even though EMD led to a slight overestimation of A. Mortazavi et al. [74] employed 
NEMD to estimate the thermal conduction in graphene-reinforced nanocomposites. 
In their study, epoxy consisted of DGEBA resin and DETA (diethylenetriamine) or 
DDS curing agents. They found that the thermal conductivity of a single-layer gra- 
phene immersed in the epoxy matrix at atmospheric pressure and room temperature 
decreases by around 3096 with respect to the value of pristine graphene, while the 
type of hardener has no considerable effect. Furthermore, they also investigated the 
impact of pressure and epoxy-graphene covalent bonds on À of graphene: no signifi- 
cant change in the À of graphene was observed by increasing the pressure up to 14 
GPa, while the formation of 596 covalent bonds led to 6596 reduction in the thermal 
conductivity of graphene. Clearly, such methodologies can be successfully adopted 
in a large variety of polymeric matrices: for example, Alaghemandi et al. [75] used 
RNEMD to investigate the À of SWCNT-polyamide (PA) nanocomposites along both 
parallel and perpendicular directions with respect to the nanotube axis, finding a good 
agreement with experiments [76, 77]. 


Thermal boundary resistance 


The thermal boundary resistance, also known as Kapitza resistance or interfacial ther- 
mal resistance (R;) [78], at filler-filler and filler-matrix interfaces is one of the most 
significant bottlenecks in the overall thermal conduction through polymer nanocom- 
posites [79, 80]. This thermal resistance arises mainly from the weak van der Waals 
interactions between adjacent fillers, and between fillers and the surrounding matrix. 
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Such weak interactions considerably hinder the specific heat flux (g) across the inter- 
face due to phonon scattering, which in turn results in a localized interfacial tempera- 
ture jump (AT). The inverse of thermal boundary resistance, also known as thermal 
boundary conductance, Gx = 1/ Ry, is defined as 


q 


Gi = AT 


(3.8) 
MD simulations are typically used to investigate thermal boundary resistance, since 
it arises from atomistic details of interfaces. Generally, two NEMD approaches have 
been used in the literature to study interfacial thermal properties of nanocompo- 
sites [81]. 

In the first approach, the simulated system is divided into multiple slabs ac- 
cording to the geometrical characteristics of the box. On the one side, a constant 
heat flux is added to the atoms of nanofiller (heat source); on the other side, a con- 
stant heat flux is removed from the outermost region of the computational domain 
occupied by the matrix (heat sink). This thermal power could be extracted or intro- 
duced by means of a velocity rescaling procedure, which must conserve the total 
linear momentum of the system. As an example, this method was adopted to com- 
pute the thermal boundary resistance between the CNT and octane liquid depicted 
in Figure 3.9. In those simulations [81], the continuous addition (CNT) and extraction 
(outermost region occupied by octane) of a constant heat flux leads to a steady-state 
temperature profile throughout the system, where the average temperature in each 


20 
i, 3 15 
Liquid mos 5 
TENE $ 
Nanotube oA 
5 
40 
45 
2 all - ‘ N29 


Figure 3.9: Schematic of division of carbon nanotube/octane system into concentric cylindrical 
slabs to investigate thermal interfacial effects. Reprinted from Ref. [81] with the permission of AIP 
Publishing. 


46 — Rajat Srivastava et al. 


slab can be computed by eq. (3.5). An example of temperature profile measured in 
this system is shown in Figure 3.10: the average temperature of the CNT is around 
575 K, while the liquid temperature ranges from 320 K (close to the CNT surface) 
down to 300 K (heat sink). The significant temperature drop at the nanotube-octane 
interface proves the existence of Kapitza resistance, which is then computed by 
eq. (3.8) since both q and AT can be extracted from simulations. 
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Figure 3.10: Relation between the temperature gradient and distance from the center of CNT at 
constant heat flow (see configuration in Figure 3.9). Reprinted from Ref. [81] with the permission of 
AIP Publishing. 


In the second approach, the whole system is first equilibrated at the initial tempera- 
ture T(0); after that, the filler is heated up to a predetermined temperature T;(0) = 
T(0) + AT(0), while the matrix temperature is kept constant at the T,,(t) = T(0) value. 
Such a temperature gradient can be determined by a proper distribution of thermo- 
stats in the computational domain. The system is then allowed to relax without any 
heat source or sink. As depicted in Figure 3.11, the temperature difference between 
the filler and the surrounding matrix (AT(t) = T;(t) - Tm(t)) tends to decay exponen- 
tially with time according to the Newton's law of cooling, namely 


AT(t) -AT(o)exp[ - j. (3.9) 


If the conduction resistance of the filler is considerably lower than the matrix one 
(which is typically the case of carbon nanofillers), the characteristic decay time (T) 
can be related to the area of filler-matrix interface (Aint), the filler heat capacity (c) 
and the (R,) at filler-matrix interface as [82] 
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Ry = Am (3.10) 
c 

Due to the paramount importance of heat transfer across nanoscale interfaces on 
the effective thermal properties of nanocomposites, several researches have used 
MD simulations to achieve a mechanistic understanding of the parameters that 
could affect thermal boundary resistance, such as curing degree of polymeric ma- 
trix, surface functionalization of fillers and presence of covalent or noncovalent 
bonds between fillers and the surrounding matrix [83-100]. 
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Figure 3.11: Instantaneously heated nanotube at a desired temperature is allowed to relax without 
heat source or sink, then (a) the temperature of the CNT, (b) liquid temperature and (c) difference 
between the temperature of CNT and liquid as a function of time are shown (see configuration in 
Figure 3.9). Reprinted from Ref. [81] with the permission of AIP Publishing. 


As an example, Varshney et al. [83] analyzed by MD simulations a nanocomposite 
made of epoxy polymer (EPON-862/DETDA) reinforced with SWCNTs. They reported a 
20% enhancement in the Gx at filler-matrix interface by increasing the curing degree 
of epoxy. Such behavior was explained by the increased structural rigidity and non- 
bonded interactions of the cured epoxies, since the curing process generates a denser 
matrix structure. Huang et al. [94], instead, investigated the impact of covalent bonds 
(cross-links) between SWCNT and PEK (poly ether ketone) on the thermal boundary 
conductance at their interface. Considering a simulation domain similar to the one in 
Figure 3.9, the number of cross-links between SWCNT and PEK was changed in the 
range 0-80 (0% — 6.25% functionalization degree), and the Gy was then computed in 
each case. Results in Figure 3.12 show that the temperature jump at CNT-PEK interface 
decreases with the number of cross-links; thus, Gi increases from 23.58 (no cross- 
links) to 236.44 MW/m? K (80 cross-links). Finally, Wang et al. [89] employed RNEMD 
to determine how different surface functionalizations of graphene could affect the 
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Figure 3.12: (a) Radial temperature distribution in the system under steady-state condition for 
seven simulated cases: no cross-links (red, 096), 1 cross-link (blue, 0.0896), 5 cross-links (black, 
0.39%), 10 cross-links (cyan, 0.78%), 20 cross-links (dark green, 1.56%), 40 cross-links (brown, 
3.12%) and 80 cross-links (green, 6.25%). (b) Temperature difference and thermal boundary 
conductance at the CNT-PEK interface for different studied cases. Reprinted from Ref. [94] with the 
permission of AIP Publishing. 


interfacial thermal resistance at graphene- polymer interface. To this aim, various types 
of covalent (butyl, carboxyl and hydroxyl) and noncovalent (1-pyrenebutyl, 1-pyrenebu- 
tyric acid and 1-pyrenebutylamine) functional groups were considered at different cov- 
erages, namely, the ratio between the number of functional groups to the total number 
of carbon atoms in the graphene sheet. Furthermore, the effect of isotope doping and 
acetylenic linkage in graphene was also assessed. Simulation results for covalent and 
noncovalent functionalizations are reported in Figures 3.13 and 3.14, respectively. It is 
noteworthy to mention that the thermal boundary resistance values in these figures are 
normalized using the value of R; in the case of no functionalization on graphene, that 
is, Ryo = (0.713 + 0.036) x 10 -8 m? K/W. While both increasing covalent and noncova- 
lent functionalizations led to Ry decrease, no significant effects were instead ob- 
served with isotope doping or acetylenic linkages in the graphene filler [89]. 


Mechanical properties 


MD simulations are also useful to compute mechanical properties of composite ma- 
terials, such as the elastic moduli and the whole stress-strain response of a nano- 
metric specimen [101, 102]. Discrepancies between MD and experimental results 
may be due to the faster strain rate simulated, the effect of dislocations at higher 
scales, the local high variability of temperature and the possible unwanted residual 
stresses from cross-link process [103]. 
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Figure 3.13: Dependence of the relative thermal boundary resistance on coverage percentage of 
different covalent functional molecules. Reprinted from Ref. [89] with the permission of American 
Chemical Society. 


—9—1-Pyrenebutyl 
7—1 -Pyrenebutyric acid 
79-1 -Pyrenebutylamine 


Rul Ryo 


“0.0 0.5 1.0 15 2.0 2.5 
Functionalization coverage (96) 


Figure 3.14: Dependence of the relative thermal boundary resistance on coverage percentage of 
different noncovalent functional molecules. Reprinted from Ref. [89] with the permission of 
American Chemical Society. 


In general, the elastic properties of a material along one specific direction can be 
quantified by the Young's modulus (Em) as 


0 - Eg£, (3.11) 
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Figure 3.15: Longitudinal and shear deformation of a simple volume of material. 


where the stress (c) is a measure of the average force applied over an equiva- 
lent area (F/A), and the strain (e) is the ratio of deformation with respect to the ini- 
tial length (see Figure 3.15). Such relation can be generalized to three-dimensional 
space to relate all the components of the stress and strain tensors (generalized 
Hooke's law): 


gi- Cia Ex, (3.12) 


where C’ iig are the components of the fourth-order stiffness tensor of the considered 
material. Note that the C’ tensor is symmetric (oj; = 0j; and thus Cj; = Cim, similarly 
for the strain tensor), and that the Voigt’s notation [104] is typically adopted to sim- 
plify eq. (3.12): 


Oj = Cije;, (3.13) 


with i,j=1,2, ...,6. Therefore, for any linearly elastic material, the stiffness matrix 
can be written as 


Cu Co Ch Cu Cis Cie 
Cu Cy C5 Ca Cs C% 
wa |S Ce Ge Gy Ge Ce ie 
Cá Ca Cas Ca Cas Cie 
C5 Co Css Cu Css Cu 
Cor Cez Ces Ce Cos Cee 


MD simulations can be used to apply a deformation & =[£1 £2» €33 2€23 2£13 2€12] = [£1 
E2 €3 E4 E5 Es] to a representative unit cell of the material under investigation, and 
then to measure the respective residual stress in each direction [105]. Given these 
values, eq. (3.13) allows then to obtain the 6 x 6 elastic stiffness matrix (C^) of the 
simulated material. 
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The expression of C’ is simpler in materials characterized by symmetry, since 
some of its coefficients are equal to zero. For example, the elastic stiffness tensor of 
orthotropic materials can be reduced to 


Cu Cp Cp 0 
Cn Cy Cy 0 0 
Cy C5 C$ 0 0 
0 0 0 C4 O0 
o D 0 0 “Ge 0 
o 0 0 0 0 Cg 


C= (3.15) 


namely to nine different stiffness coefficients (Cj; = Cx, C13 = C31, C5 = C32). Instead, 
the elastic properties of homogeneous isotropic linear elastic materials are deter- 
mined by two moduli only, namely Àm and um: 


1 " 
Hm 3 (Cha + Css + Cic) (3.16) 
and 


e J J 2 Je "å 
Am = = (Cu + Cz + C33) - 5 (Cha + Css + Co), (3.17) 


Wile 


and thus 
1 / 
Àm s 2g = 3 (Cir + Cop + C33). 


In this case, the elastic stiffness tensor simplifies to [103] 


Àm * 24g, Àm Àm 0 0 O 
Àm Àm + 24, Am o 0 O 
do Àm Am Àat2u, O0 O 0 (338) 
0 0 O umm 0 0 
0 0 0 O um 0 
0 0 0 O 0 n, 


Once Àm and up are determined, the following relations for homogeneous isotropic 
linear elastic materials hold: 


Gm - Hg, (3.19) 


Km = Àm + sem (3.20) 
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Àm 
v= = 3.21 
Wn + Hr) = 
3Am + 24 
Ey = — en, :22 
m= bn TR (3.22) 


where Gn is the shear modulus, which regards forces parallel to one surface of the 
material; Km is the bulk modulus, which regards the isobaric compression of the 
material; v is the Poisson’s ratio, which relates the expansion of the material in one 
direction while it is compressed in another; Em represents the Young’s modulus, 
which is a measure of the linear elastic stiffness of the material. 

MD can be also used to mimic the experimental tensile test of nanomaterials. 
Typically, the system is first equilibrated; then, a small strain is applied along one 
longitudinal direction of the computational box with a ramped function of time 
(see Figure 3.16), and the resulting strain measured [106, 107]. 


Figure 3.16: Longitudinal deformation of a simulation box. 


Electrical properties 


Polymeric materials are typically characterized by very low (or no) electrical con- 
ductivity; however, the insertion of carbon fillers — such as graphene sheets, CNTs 
or carbon fibers (CFs) - can improve the electrical properties of the resulting poly- 
mer nanocomposites. 

A limited volume fraction of CFs has been observed to significantly improve 
electrical conductivity of polymer nanocomposites [108, 109]. This is due to the for- 
mation of continuous conductive networks of CFs in the polymer matrix, thus al- 
lowing a transition from nonconductive to conductive electrical properties. The 
percolation threshold is the critical volume fraction of conductive fillers in the 
nanocomposite above which conductive pathways in the system always occur. This 
threshold determines the transition from electrical insulating to conducting behav- 
ior, and it is therefore the crucial figure of merit to understand the electrical features 
of nanocomposites. Furthermore, researchers have reported that the conductive be- 
havior of nanocomposites depends on the spatial distribution of CFs as well. On the 
one hand, high volume fractions of fibers poorly dispersed in the polymer matrix 
may lead to fiber agglomeration, which possibly results in electrical conductivity en- 
hancement (but degraded mechanical properties). On the other hand, spatially uniform 
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distributions of fibers at low volume fractions are unfavorable for the electrical conduc- 
tance of nanocomposites. 

Molecular simulation is an efficient tool to study the electrical properties and 
percolation mechanism in nanocomposites, since it allows determining whether 
CFs are forming conductive networks within the polymeric matrix or not. According 
to the electric tunneling theory, two carbon fillers can be considered as electrically 
connected to each other if their shortest distance is less than a critical gap, which is 
known as tunneling distance. The electrical resistance offered by the thin (i.e., less 
than the tunneling distance) insulating polymer matrix between the conductive fi- 
bers is therefore known as tunneling resistance. 

The generalized formula for electric tunnel effect through a potential barrier be- 
tween electrodes separated by a thin insulating film was derived by Simmons in 
1963 [110]. Hence, if the thickness of insulating film between the conductive fibers 
is treated as uniform, the curvature of fibers is neglected, the potential barrier is 
taken as rectangular and the image forces are included, the current density (j, con- 
sidered in A/cm?) flowing through the thin film of polymer matrix between the fi- 
bers is given as [110] 


. 62x10? 
I= — As 
-($ + U) exp(- 1.025As($ + U)"®)], (3.23) 


[p exp(- 1.025Asp"”) 


where 


S2(S— S1) 


scm (3.24) 


5.75 
$= bo - (U/2s)(s1 +52) - In 
s is the thickness of the insulating polymer film (in A), As = (s2 - s1) is the difference 
of the limits of barrier at Fermi level (in A), x is the dielectric constant of insulating 
polymer film, U is the voltage across the insulating polymer film (in V) and @, is 
the height of potential barrier (or work function) (in V). If U < $,, then 


$1 = 6/KPo 
s2 - s|1 - 46/(3@ Ks + 20 - 2Uks)] + 6/x$. (3.25) 
Moreover, the voltage across the insulating polymer film can be calculated as 


es 


ame AckEo’ 


(3.26) 


where e, C, Ac and £o are electron charge, capacitance, area of the insulating film in 
contact and permittivity of free space, respectively. The electric tunneling resis- 
tance can be finally computed as 


54 —— Rajat Srivastava et al. 


U 

Reunnel = j Ac . (3.27) 
Li et al. [111] performed Monte Carlo simulations to understand the effect of tunnel- 
ing resistance between CNTs on the electrical conductivity of the whole polymer 
nanocomposite. Their results indicate that the tunneling resistance increases with 
the thickness of insulating polymer film (see Figure 3.17). They also reported that 
the conductivity of nanocomposites drops below 10? S/m, namely the nanocom- 
posite behaves as an electrical insulator, when the tunneling resistance is larger 
than 10” Q. The thickness of polymeric film corresponding to this Q value is found 
to be 1.8 nm, which is then considered as the maximum tunneling distance. 
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Figure 3.17: Tunneling resistance as a function of the thickness of insulating polymer film between 
carbon nanotube fillers in a nanocomposite. Different diameters of nanotubes and matrix materials 
are reported. Reprinted from Ref. [111] with the permission of AIP Publishing. 


Hu et al. [112] studied the effect of CNT aspect ratio on the electrical properties of 
nanocomposites by means of three-dimensional statistical percolation modeling. 
They reported that the higher aspect ratio of CNTs leads to larger aggregates of fillers 
in the nanocomposite, which ease the creation of paths with high electrical conduc- 
tivity. As a result, both lower percolation threshold and higher electrical conductivity 
of nanocomposites are reported for larger aspect ratio of fillers [112]. Furthermore, the 
electrical percolation threshold is found to be significantly influenced by the inter- 
tube van der Waals interaction and tunneling resistance for nanocomposites with 
low CNT aspect ratio (i.e., length/diameter ratio «200), whereas negligible effects are 
observed with high aspect ratio [113]. 

Researchers agree that the tunneling resistance, which shows values several orders 
of magnitude larger than the intrinsic resistance of carbon fillers, plays a dominating 
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role in controlling the electrical behavior of nanocomposites [111-114]. Thus, to en- 
hance computational efficiency and simulate larger domains, the intrinsic resistance of 
nanofillers is often neglected in simulations of nanocomposites with well-dispersed 
fillers [114]. 


3.2.2 Mesoscopic models 
Thermal properties 


The most popular approaches to study effective thermal properties of composites at 
mesoscale are off-lattice Monte Carlo, dissipative particle dynamics (DPD) and lat- 
tice Boltzmann method (LBM). 

Duong and coworkers [115] developed an off-lattice Monte Carlo simulation ap- 
proach to investigate the effective thermal conductivity (gy) of CNT-reinforced 
nanocomposites, taking into account also the thermal boundary resistance (Ry) at 
CNT-matrix interface. In their model, heat transfer comes from the random motion 
of a large number of discrete thermal walkers. The proposed model has been widely 
utilized to investigate the impact of fibers distribution, morphology, volume frac- 
tion as well as Ry value at fiber-matrix interface on the Agg of nanocomposites 
[116-119], nanofluids [120, 121] and aerogels [122] based on CNTs, with a good capa- 
bility to reproduce experiments. Kui et al. [123] modified Duong's model to evaluate 
the thermal properties of composite materials reinforced by graphene sheets. Later, 
Gong et al. [124—127] extended the initial model to evaluate the thermal conductivity 
of multiphase composite materials as well. As an example, they studied the Agg of a 
polymer composite consisting of polyether ether ketone (PEEK) matrix reinforced 
by SWCNT fibers and tungsten disulfide (WS) nanoparticles [125]. They explored 
the impact of interfacial thermal resistance at CNT-PEEK and WS;-PEEK interfaces, 
CNT morphology (i.e., diameter and aspect ratio of SWCNTs) and SWCNT orienta- 
tion (i.e., random, parallel or perpendicular to heat flux) on the effective thermal 
conductivity of composites. According to their findings, Agg increases with SWCNT 
concentration, while it decreases with larger Kapitza resistances at the different 
interfaces. Regarding SWCNT orientation, fibers parallel and perpendicular to the 
heat flux direction show the best and the worst effect on the overall thermal perfor- 
mance, respectively. They also reported that SWCNTs with higher aspect ratio and 
smaller diameter generate higher thermal conductivity enhancements at given 
mass fraction of fillers. 

In DPD, the simulated particles (beads) represent whole molecules or fluid re- 
gions [128]. For instance, Zhou et al. [129] coupled DPD and smoothed particle hy- 
drodynamics (SPH) to study, at mesoscopic level, the Aep of polymer composites 
reinforced by CNTs. This numerical model was used to investigate the impact of dif- 
ferent parameters, such as filler dispersion, volume fraction, length and matrix 
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characteristics, on Aeg of nanocomposites. Their results show that Ag changes 
quadratically with the volume fraction of fillers, with both random and aligned dis- 
tributions of fibers. 

LBM is another method that can be used to study the thermal properties of com- 
posite materials at mesoscopic level. The most significant advantage of LBM is the 
easy implementation of different interparticle interactions and complex boundary con- 
ditions [130]. Wang and coworkers [131] developed a three-dimensional LBM model to 
calculate the thermal conductivity of composites reinforced by CFs, in good agreement 
with experimental evidence. Chiavazzo and Asinari [132], instead, computed through 
LBM the Aeg of composites made of polypropylene polymer (matrix) and graphite par- 
ticles (fillers). Fang et al. proposed a multiple relaxation time LBM to predict numeri- 
cally the Agg of anisotropic heterogeneous materials, such as polymers reinforced by 
braiding yarns [133]. The model has been afterward applied to compute longitudinal 
and perpendicular Agg of computational domains with different volume fractions of 
fillers, considering the interfacial thermal resistance between components as well. 


Mechanical properties 


Atomistic simulation methods are well established for computing the interfacial 
properties of polymer composites at the molecular level (nanoscale). However, the 
considerable computational resources required limit the maximum length- and 
timescale of feasible simulations [134, 135]. While sophisticated strategies for speed- 
ing up complex MD simulations of systems that undergo significant configuration 
changes have been recently suggested on the basis of the equation-free method 
[136, 137], coarse-grained (CG) models have been specifically developed to simulate 
polymers [138-140] and carbon fillers (CNTs and graphene) [141, 142]. 

CG MD overcomes the complexity of atomistic simulations by reducing the degrees 
of freedom of particle- particle interactions, so that larger length- and timescales (from 
nanoscale to mesoscale) can be afforded [138, 143]. For the development of CG molecu- 
lar models, chemically connected atoms are grouped together to form superatoms or 
CG beads that contain sufficient information to reproduce the chemical features of the 
molecular structure (see, for instance, Figure 3.18) [139]. Moreover, several methods 
have been proposed to parametrize the force field among CG beads starting from their 
respective atomistic details [134, 138, 139]. Arash et al. [143, 144] developed a CG model 
of CNT/polymer systems, with the aim to compute the mechanical behavior of polymer 
nanocomposites. Their results confirm that CNT reinforcement significantly enhances 
the Young's modulus of the composite, and that CG models provide comparable accu- 
racy and lower computational cost than atomistic simulations. Later, Mousavi et al. 
[145] introduced nonbonded interactions between polymer chains, CNTs and polymer 
matrix in the CG model to investigate the influence of cross-linking, weight fraction 
and distribution of CNTs on the elastic properties of polymer nanocomposites. 
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Figure 3.18: Coarse-grained model of polystyrene (red and gray beads) from atomistic details 
(white and cyan sticks). Reprinted from Ref. [139] with permission of American Chemical Society. 


The effect of dispersion, alignment and morphology of nanofillers on the proper- 
ties of polymer composites can be studied by DPD as well. Zhou et al. [146] studied 
the dispersion and alignment of CNTs using DPD simulations. Numerical experiments 
demonstrated that CNTs tend to align themselves in the polymer matrix with increas- 
ing nanotube length and volume fraction, especially in well-dispersed systems. 
Instead, a DPD study on graphene/PMMA (Poly(methyl methacrylate)) composite 
showed that, while pristine graphene and PMMA are immiscible, high concentrations 
of graphene functionalizations lead to better filler dispersions in the PMMA matrix 
[147]. Furthermore, when the volume fraction of fillers is larger, graphene needs to be 
functionalized more to reach good dispersion. Similar results were obtained by Lin 
et al. [148], who reported better dispersions of covalently functionalized graphene fill- 
ers in PMMA matrix compared to nonfunctionalized ones. 


Electrical properties 


The RVE concept is also applicable to predict the mesoscale electrical properties of 
nanocomposites. The electrical conductivity of carbon-based polymer nanocompo- 
sites mainly depends on [149—151]: 
— the intrinsic electrical conductance of CFs; 
— the normal transport through conductive networks of carbon fillers in contact 
with each other within the polymeric matrix; 
— the hopping transport due to the tunneling of electrons between CFs separated 
by polymer matrix of thickness less than the tunneling distance. 
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To compute the electrical conductivity of nanocomposites at mesoscale, CFs should 
be first randomly distributed within a RVE of polymer nanocomposite. These simu- 
lations typically assume that fibers have low electrical resistance and thus act as 
conductors, whereas polymeric matrix is considered as a highly resistive material 
with low electrical conduction. Therefore, the most conductive path for a given com- 
posite is the one minimizing the relative distance between CFs [151, 152]. Different al- 
gorithms are available to find the shortest path in a given configuration (e.g., Dijkstra 
algorithm [153]), which can be thus employed to calculate the electrical conductivity 
of the simulated nanocomposites at mesoscale. This procedure is typically repeated 
over different random distributions of fibers in the matrix to reduce statistical fluctua- 
tions in terms of composite conductivity. 

In a cuboid RVE with dimensions (Ly, Ly, Lz) containing randomly distributed 
fibers in a polymer matrix [154], each fiber can be modeled as a line segment in 
spherical coordinates (see Figure 3.19): 


x; =X; + li sin 6; cos vi, 
yt =y} + li sin 0; sin qi, 
zi =Z} + li cos 0;, (3.28) 


(x?, yf, 27) 


Figure 3.19: Fiber modeled as a line segment in spherical coordinates. 


where i is the index of the i^ fiber, (x7, y3,z$) and (x5, y5, 2?) are the start and end 
point coordinates of the fiber. The length, polar and azimuthal angle of the i“ fiber 
are denoted as l;, 0; and vj, respectively. 

The start point coordinates, the polar angles and the azimuthal angles of each 
fiber can be generated randomly as 


Xi-L.& yi-Ly& Zale 
0;-2n£, i= cos !(2£ -1) (3.29) 
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where £ is an uniformly generated random number in the range [0, 1]. To better rep- 
resent experimental variability, the simulated length of fibers can follow a Weibull 
distribution [155]. It should be noted that the end points of fibers might initially lie 
outside the RVE cuboid: in that case, periodic boundary conditions are applied to 
relocate them inside the cuboid. 

The electrical conductivity of nanocomposites can be estimated once the series 
of resistances in the percolating network is computed, including both intrinsic re- 
sistance of fibers (Rj) and their contact resistance (Reontact). According to the Drude 
model, the Rj of a fiber (i,j) with lj length, D diameter and ofiber intrinsic electrical 
conductivity (see Figure 3.20) can be calculated as [154] 


(3.30) 


Electrodes 


Figure 3.20: Network of electrical resistors in the RVE of a nanocomposite reinforced by nanofillers, 
for example, CNTs. Reprinted from Ref. [151] with permission from Elsevier. 


Considering a marginal effect of temperature, the contact resistance Roos between 
two fibers can be instead estimated as [150] 


h 


T uL 31 
2e?MTp ee) 


Reontact = 
being h the Planck constant, e the electron charge, M the number of conduction chan- 
nels and T, the transmission probability for the electron to tunnel through the polymer 
layer between fibers. According to Bao et al. [150], the transmission probability Tp can 
be estimated by the Wentzel-Kramers-Brillouin approximation [156] as follows: 


exp( - quit) OxdzD-dyaw 
tunnel 


Tp= (3.32) 


exp(- L2) D+ dyaw <d<D+ dai 


tunnel 
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where d is the minimum distance between adjacent fibers, dygw is the van der 
Waals separation distance and dunner is the tunneling length. The latter can be de- 
scribed as [150] 


h 
tunnel = /8m.AE > 
e 


where fi = h/2m is the reduced Planck constant, me the mass of electron and AE the 
energy barrier. Note that in eq. (3.32), tunneling effects are neglected if the thick- 
ness of polymer matrix between contiguous fibers is larger than deut, namely, the 
considered cut-off distance. 

Bao et al. [150] studied the effect of CNT alignment in polymer matrix on the critical 
percolation threshold, and their results revealed that the maximum electrical conduc- 
tivity is obtained in nanocomposites with partially aligned CNTs. Similarly, Monte Carlo 
simulations by Zeng et al. [157] reported a clear dependence of the electrical percolation 
threshold on CNT length, waviness, distribution anisotropy and volume fraction. 
Improved three-dimensional percolating network theories can also be used to under- 
stand how the deformation of CNTs at crossed nanotube junctions influences the elec- 
trical conductivity of nanocomposites [158]. In detail, the study by Gong and colleagues 
[158] suggests that CNT deformation at crossed nanotube junctions determines a large 
increase in the intrinsic resistance of the nanotube while only a limited decrease in the 
CNT-CNT contact resistance, therefore, leading to enhanced overall resistance at the 
junction. Furthermore, numerical simulations based on the Dijkstra algorithm showed 
that increased contents of MWCNTs in polymer nanocomposites result in higher electri- 
cal conductivities, in good quantitative agreement with experiments (see Figure 3.21) 
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Figure 3.21: Simulation and experimental results of electrical conductivity of CNT-based polymer 
composites as a function of weight percentage of MWCNTs. Reprinted from Ref. [152] with 
permission from Elsevier. 
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[152]. Finally, Rahatekar et al. [159] performed DPD simulations to explore the influ- 
ence of aspect ratio and fiber alignment on the electrical percolation threshold of 
nanocomposites. 


3.2.3 Continuum models 
Thermal properties 


The effective thermal conductivity (Ag) of polymer nanocomposites at macroscale 
can be predicted by two approaches, namely analytical or simulation ones. Analytical 
models, for example, effective medium approximation, generally provide quick evalua- 
tions of the properties of composites, but cannot take into account neither the interac- 
tions between adjacent inclusions nor the actual morphology of composite. Further 
details on the analytical models for effective thermal conductivity of nanocomposites 
can be found in a recent review article by Zhai et al. [160]. 

Due to the current improvement of computational tools and numerical algo- 
rithms, finite element method (FEM) has been increasingly utilized to compute the 
thermal conductivity of composite materials at continuum level. FEM calculates the 
Aeg of nanocomposites by solving numerically the Fourier's law for conduction at 
steady state, and this procedure is - in principle — capable to consider different 
morphologies (e.g., fiber distribution) of the composite material by proper meshing 
strategies. FEM simulations of composites are typically limited to an RVE of the sys- 
tem. The RVE is built in such a manner that the smallest constituent that can influ- 
ence the first-order macroscopic behavior of the system is consistent with the 
simulated length scale; then, the results obtained from the RVE are extended to de- 
velop a full-scale model. For example, Figure 3.22 depicts 3D cubic RVEs of compo- 
sites with different inclusions, namely cylindrical, platelet-like and spherical fillers. 
The major drawback of FEM in comparison with analytical approaches remains the 
high computational cost, as well as the modeling complexity given by systems in- 
cluding fibers with high volume fraction or aspect ratio. 

Ramani and Vaidynathan [162] carried out an automated finite element analysis 
to compute Àe of composite materials. The employed FEM model allowed to study 
the effect of various parameters, such as fiber aspect ratio, volume fraction, disper- 
sion, orientation and fiber-matrix interfacial thermal resistance on the resulting 
thermal behavior of composites. Numerical results were compared to analytical 
models and experiments, finding a good agreement especially at high volume frac- 
tions. Ahmed and Masud [163] employed FEM to study the effect of geometrical pa- 
rameters on the Ag of polymer composites reinforced by MWCNTs. Their study 
revealed that the À of nanocomposites increases with the aspect ratio of nanotubes, 
whereas interface thermal resistances do not have a prime contribution in lowering 
the thermal conductivity at the macroscopic scale. Instead, Li et al. [164] quantified 
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Figure 3.22: Temperature fields in 3D representative volume elements of different composites 
simulated by FEM method. Reprinted from Ref. [161] with permission from Elsevier. 


the influence of SWCNT-matrix interfacial thermal resistance (Renr-m) on Aeg of 
SWCNT/polymer composites. Their results showed that the Ag of composite can be 
lower than the one of pristine polymer matrix if Renr-m is greater than a critical value. 


Mechanical properties 


The FEM can be used to understand the macroscale mechanical properties of nano- 
composites, for example, the elastic moduli. At molecular level, the RVE of a ge- 
neric nanocomposite is a heterogeneous medium consisting of polymer matrix, 
carbon nanofibers and interface region. However, at continuum level, such hetero- 
geneous RVE should be homogenized to compute the effective material properties 
through micromechanical analyses (see Figure 3.23). Notice that the homogeniza- 
tion procedure of RVEs to calculate the effective elastic moduli of nanocomposites 
is equivalent to average the mechanical properties of heterogeneous nanocompo- 
sites [165]. 

The constituents of the RVE can then be assumed as isotropically symmetric 
within the nanocomposite [166], and the constitutive relation given by the general- 
ized Hooke's law: 
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Figure 3.23: Homogenization of a representative volume element of nanocomposite, which 
includes the nanofiller (CNT), the CNT/epoxy interface and part of the surrounding epoxy matrix 
into a representative, homogeneous fiber of the composite to be considered for micromechanical 
analysis. Adapted from Ref. [166] with permission from Elsevier. 


Gij = Cij Ekb ijklz123 (3.34) 


being Cia the elastic stiffness tensor of the RVE, oj and €, the stress and strain ten- 
sors, respectively. The effective stiffness coefficient (Cha) of the homogenized multi- 
constituent nanocomposite can thus be written as 


Oy = Cia Ext (3.35) 


being 0; and £j the stress and strain tensors averaged over the RVE volume (Vnyz), 
respectively, that is: 


1 
Ojj- — ojdV, 
OF jdV. (3.36) 
Verve Jvgvg 
£g = ——— £ydV. (3.37) 
Vave Jvpyg 


The effective elastic moduli of nanocomposite are finally calculated using the aver- 
aged stress and strain tensors as follows: 


Ci = ~ i-j-k-l (Young's modulus) (3.38) 
" 0i . . 
E XE (i=k)+(j=1) (Shear modulus) (3.39) 


Instead, the Poisson’s ratio of nanocomposite is computed as 


ii 
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Odegard et al. [167] proposed a method to transform the discrete molecular struc- 
tures of nanotube fillers, polymer matrix and CNT/polymer interface at nanoscale 
into an equivalent continuum model of composite. Liu and Chen [168] computed 
the effective mechanical properties of CNT-reinforced composites using a nanoscale 
RVE and the FEM. They showed that the stiffness of composite increases several 
times with volume fractions of CNT nanofillers up to 596. Instead, an equivalent 
continuum modeling study by Huang and Rodrigue [169] found that both aspect 
ratio and length of CNT fillers have a significant influence on the mechanical prop- 
erties of polypropylene-based composites. In a similar study, Kumar and Srivastava 
[170] reported that the elastic stiffness of polymer matrix was improved by carbon- 
based filler reinforcements. This study also concluded that graphene reinforce- 
ments provided the best in-plane stiffness properties of nanocomposite, whereas 
CNTs the best out-of-plane ones. Finally, a continuum simulation using multiscale 
homogenization procedure was applied by Golestanian and Gahruei [171] to under- 
stand the effect of CNT waviness on the elastic properties of CNT-reinforced polymer 
composites. The simulation study revealed that wavier CNTs reduce elastic modulus 
of nanocomposite along the longitudinal direction, while they increase elastic mod- 
ulus along the transverse one. 


Electrical properties 


The effective electrical properties of nanocomposites at the micro- and macroscale 
are substantially affected by phenomena at molecular level, such as tunneling 
mechanism and formation of conductive networks. 

In case of a 1D conductor with / length and A, cross section, the effective electri- 
cal conductivity of composite (oe) can be computed by measuring the electric current 
(I) flowing through the computational domain, given a certain electric potential (U), 
that is 


Il 


= s 41 
A.U ap 


Ge 


Instead, in case of 3D anisotropic materials, a more general expression should be 
considered: 


j=0eE, (3.42) 


where j being the current density vector, E= — VU the electric field vector under 
electrostatic conditions and Ge the effective conductivity tensor. Typically, FEM anal- 
ysis with micromechanical corrections is adopted to compute effective electrical prop- 
erties of nanocomposites at continuum level as, for example, depicted in Figure 3.24. 
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Figure 3.24: Electrical potential within a nanocomposite evaluated by FEM analysis. Reprinted from 
Ref. [172] with permission from Elsevier. 


Seidel and Lagoudas [173] developed a micromechanical model to estimate the 
electrical properties of CNT-polymer nanocomposites. They reported that the large 
increment in the de of MWCNT-based nanocomposites at low nanotube concentra- 
tions is mainly due to formation of conductive networks, whereas in SWCNT-based 
nanocomposites, both electron hopping and conductive network formation play a rel- 
evant role [173, 174]. The macroscale piezoresistive properties of CNT-polymer com- 
posites along both transverse and axial directions were studied instead by Ren and 
Seidel [37], considering both electric tunneling and inherent piezoresistivity of CNT in 
the selected RVEs. In that study, the macroscale piezoresistive response of different 
mesoscale CNT dispersion scenarios in the nanocomposite was predicted using 
computational micromechanic techniques based on the finite element analysis. The 
authors found that the electric tunneling mechanism dominates the overall piezore- 
sistive properties of CNT-polymer composites along both transverse and axial direc- 
tions, whereas the inherent piezoresistivity contribution is more evident along the 
axial direction rather than along the transverse one [37]. The influence of electric 
tunneling was found to be negligible in composites with well-dispersed distributions 
of CNTs, where the intertube distance was larger than the tunneling distance; con- 
versely, agglomerated CNT dispersions eased electric tunneling mechanism. 
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3.3 Perspectives 


The idea to combine different materials to develop an engineered material having 
unique capabilities, distinct from their constituent materials, is very fascinating, 
and it has been well known by long time. Nowadays, the possibilities offered by 
materials processing and, sometimes, nanotechnology have revamped the empha- 
sis on composite materials. In particular, carbon-based composites have already 
proved to provide remarkable improvements in terms of stiffness, durability, strength 
and lightweight characteristics in some industrial sectors, for example, automotive. 

In spite of the previous successes, some industrial challenges remain for compos- 
ite materials. According to a recent commentary by Materials Today [175], mass pro- 
duction of composites is a key aspect that must still be developed before composites 
will prove to be useful for many new applications. The lack of mass production/auto- 
mation techniques sometimes reveals the lack of knowledge of the composite 
manufacturing process, which forces to rely on try and error in developing new prod- 
ucts [175]. Clearly, this approach, which is still dominant in many applications, can- 
not be properly standardized. Hence, the lack of standardization, design standards 
and good technical support/backup from suppliers tends to reduce the composite ac- 
ceptance in an industrial context [175]. Composite materials without a warranty are 
not well accepted by engineers, who have to deal with high-fidelity standards in the 
design and production of complex products. Things become even worse in case of 
recycling CFs, as envisioned by circular economy. End-of-life fiber-reinforced compo- 
sites are particularly challenging for implementing circular economy, even though 
some large-scale demonstration has been already funded by the European 
Commission [176]. All the above factors significantly increase the price of carbon- 
based composites with regard to their competitors, namely steel, concrete and ther- 
moplastics, and this induces a loss of market share [175]. Price reduction can hardly 
be achieved by combining raw materials, because this intrinsically requires some ad- 
ditional processing costs. However, one promising alternative option consists in in- 
creasing the value chain by adding new functionalities for composite materials. This 
leads to the so-called smart materials (see Table 3.1 and references therein). Materials 
with self-sensing and structural health monitoring capabilities, devices with fast trig- 
gering shape memory, thermally driven actuators and high-performance all-solid-state 
supercapacitors are good examples of this strategy. On the other hand, this strategy 
leads necessarily to even more scientific challenges, because high precision and repro- 
ducibility of the arrangement of CFs into the composites are required. 

The previous discussion reveals that the exploitation of carbon-based compo- 
sites requires addressing many scientific challenges, which could benefit from sys- 
tematic and state-of-the-art materials modeling. Essentially, modeling should 
address two main issues: (i) elucidating the atomistic details of interface between 
CF and polymer and (ii) taking into account properly the extreme variability in the 
CF orientation within the composite. 
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Table 3.1: Modeling methods for computing thermophysical properties of carbon-based 
nanocomposites of interest as smart materials. PDMS is polydimethylsiloxane; PVDF is 


polyvinylidene fluoride. 


Material Modeling Properties Applications as smart 
method material 
CNT/PDMS composite FEM Structural, mechanical Thermally driven actuator 
and thermal properties — [181] 
CNT/polymer composite Multiphysics Electromechanical Strain sensor [182] 
modeling properties 


CNT/polymer 


Percolation 
network model 


Electrical conductivity 


High accuracy sensor 
[183, 184] 


PVDF/carbon fiber Multiphysics Mechanical and Self-sensing and structural 
composite micromechanics piezoelectric health monitoring [39] 

and FEM properties 
CNT sponge/shape FEM Electrical properties Fast triggering shape 
memory polymer memory material [185] 
nanocomposite 
Graphene oxide (GO)/ Quantum/MD Capacitance High-performance all-solid- 
reduced GO/water simulations state supercapacitors [186] 


The first issue is very well known and is also largely explored from the experi- 
mental point of view. For example, CF-reinforced thermoplastic composites gained 
recent interest due to ease of production and recycling compared to that of thermo- 
set composites, but surface treatments are used to increase the concentration of 
surface functional groups and hence the fiber-polymer interfacial adhesion 
[177]. Similarly, different methods of preparation for thermoset composites rely 
on functionalizing the interface between various nanofillers and epoxy, which is 
typically characterized by several experimental techniques [178]. Remarkably, 
nowadays both experimental and modeling techniques are used in order to in- 
vestigate the relationships between interfacial characteristics and composite 
properties, for example, with regard to the advantages and disadvantages of co- 
valent and noncovalent functionalization of CFs [179]. 

The second issue where materials modeling can provide some useful insights is 
due to the extreme variability in the microscopic structure of composites. Different 
arrangements of CFs within the composite, shape and orientation of the CF aggre- 
gates and dynamics response of the microstructure due to the applied loads are just 
examples of phenomena that require modeling at scales much larger than that of 
individual fiber. This may appear simply as a geometrical problem, but it leads im- 
mediately to (i) a significant increase in the computational demand and to (ii) a prob- 
lem of model reduction in postprocessing the modeling results. The second problem is 
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particularly severe, because it may hinder the opportunity to use materials modeling 
for increasing the understanding of composites under realistic conditions. In fact, un- 
derstanding requires squeezing a large amount of information in very simple design 
guidelines. Moreover, this second problem has an indirect effect on the first problem 
about affordable computational demand as well. In fact, the possibility to develop 
simple interpolation models, which can summarize the outcome of more complex 
models, is essential for any coupling and/or linking strategy. Concerning the latter 
point, analytical simplified models, which can describe the transfer of energy, me- 
chanical stresses and electricity between neighboring fibers, are highly desirable for 
developing reliable predictions about percolating networks within the composites. 
Nowadays, machine learning techniques and artificial intelligence offer solutions for 
developing analytical interpolating functions, which can be tuned on underlying de- 
tailed models and are suitable as input for up-scaled models for predicting macro- 
scopic properties of the materials. These analytical functions keep some insights into 
the investigated phenomena and, at the same time, are very fast and transferable in 
up-scaled models, which can describe large portions of the material under investiga- 
tion. This is a clear advantage with regard to previous techniques based on neural net- 
works. Hence, we envision a huge impact of such techniques in dealing with realistic 
simulations of carbon-based composites for engineering applications. 

Finally, last but not least, carbon-based composites still raise some concerns in 
terms of their toxicity, particularly when nanofibers are involved. These concerns 
are about the manufacturing processes of composites, the release of nanoparticles 
during operation and the end-of-life waste management. Clearly, the fact that some 
nanofibers, for example, CNTs, are similar in shape and size to asbestos raises fur- 
ther concern [180]. Even though materials modeling may help in assessing some 
guidelines and in clarifying some scenarios, quantitative and predictive computa- 
tional nanotoxicology is still far from being mainstream. However, it represents a 
very interesting approach, which is currently under development in many different 
fields for assessing preliminarily the impact of new nanotechnologies, and we envi- 
sion that also carbon-based composites will benefit from it. 
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4 Carbon nanotube-based materials 
for space applications 


4.1 Introduction 


Carbon nanomaterials embedded in carbon fiber-reinforced composite materials can 
improve specific performance characteristics such as electrical conductivity, resis- 
tance to mechanical fatigue and crack propagation, structural damping, Passive 
Intermodulation (PIM) reduction, Electromagnetic compatibility (EMC) shielding, re- 
machining and coating capabilities. Some of these improvements were investigated 
and demonstrated in the ESA projects *NACO-1" and *NACO-2" performed between 
2007 and 2013. This chapter outlines the work performed during these activities for 
the development of carbon fiber-reinforced polymer (CFRP) composites integrating 
carbon nanotube (CNT) buckypapers (BP) to be used in space applications. 

The ranking of space applications potentially benefitting CNT technology has been 
reviewed throughout these activities. The final activity is aimed at finally selecting an 
application with highest potential for performance improvements. For that, an optical 
CFRP mirror with comparable characteristics of conventional Zerodur-based mirrors 
featuring the inclusion of 3D nanotube structures was selected and manufactured. 

Although CFRP composites are being used for reflectors in Space applications, 
due to its outstanding mechanical and thermal properties at low mass, optical mir- 
roring properties are limited by the surface properties of the CFRP. One should note 
that each composite layer is composed of 1D and/or 2D patterns of continuous car- 
bon fibers, each one of them having about 10 um of diameter, which are combined 
in tows of thousands of carbon fibers leading to a tow a few millimeters large. This 
configuration creates an undesired composite surface roughness, in a phenomenon 
usually known as print-through effect (Figure 4.1). 

CFRP surface roughness limits its reflective properties. To afford further mirroring 
characteristics, polishing and metal coatings strategies can be followed. However, the 
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Figure 4.1: Representative image of print-through effect in a CFRP. 


former is inherently limited by the continuous form of the carbon fibers, whereas the 
latter is also limited by the print-through characteristics of the CFRP. To overcome this 
limitation, inclusion of carbon nanotubes (CNT) at the top surface of the composites 
has been proposed as a solution to reduce its print-through behavior. In addition, the 
nearly isotropic nature of the CNT top layer is expected to ease polishing and metal 
coating processing. A material configuration in which a top layer of resin rich in CNT 
is added in the CFRP lay-up allows the reduction of the print-through effect on the 
manufactured CFRP and improves the reflective properties of the top surface. 

The potential of CNTs is considered as tremendous; their intrinsic properties 
make them very attractive for Space applications where the characteristics of the ma- 
terials often represent a limit to the development of technologies. Much of the devel- 
opment efforts have been focused on retaining CNTs large aspect ratio, achieving 
good interfacial stress transfer (surface modification), good dispersion and alignment 
in the polymer matrix during processing, all for maximizing the nanotubes effect. As 
a result, a growing number of companies already supply CNT-based masterbatches or 
compounds to be used in plastic/composite processing. Nevertheless, these materials 
are expensive, limited in the type of dispersing resin systems and on the CNT weight 
loading, and very difficult to incorporate in the conventional processes due to the 
very high viscosities. 

In fact, one of the main difficulties that have been associated with the manufac- 
turing of CFRP materials containing nanotubes is the poor processability of the resin 
systems modified with CNT, due to the extremely high viscosities. It is practically im- 
possible to obtain a surface with quality when using high viscosity systems. 

CNT BP consist of a 3D structure (skeleton) of CNTs either being single-walled 
(SWNT) or multiwalled or carbon nanofibers (CNF), which show a unique combina- 
tion of inherent CNT/CNF properties while maintaining their lengths in a mat-like or 
3D structure. These structures can be an alternative to overcome dispersion issues 
usually associated with the use of CNTs in resin masterbatches. In addition, higher 
CNT contents can also be attained by the use of three-dimensional CNT structures; 
therefore, CNT in the form of BP is a promising choice for boosting CFRP properties. 
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4.2 Preparation of epoxy nanocomposites based 
on high CNT content buckypapers 


The main objective of the first study was to develop and characterize highly loaded 
CNT-based nanocomposites, by infiltrating nanotube skeletons with either organic 
or inorganic matrix. We herein present the main results obtained with an organic 
thermoset (epoxy) resin. 

A new approach was used to produce these CNT nanocomposites, where the 
CNT skeletons or papers were infiltrated with a low-viscosity epoxy resin in a closed 
mold. The CNT were first produced in thin papers called BP to form networks of 
elongated nanotubes. The impregnation was performed by infiltrating the resin 
along all the 3D CNT structure. The level of impregnation of the final structure was 
assessed by physical measurements and observations under scanning electron mi- 
croscopy (SEM). The results showed that the proposed processing approach is capa- 
ble of producing well-dispersed nanocomposites with high CNT loading (more than 
12 wt% in resin weight), which are important for developing high-performance 
structures based on CNT with good thermal and electrical conductivity. 


4.2.1 Materials 


CNT thin papers, with 80 mm of diameter, functionalized (-COOH) and solvent (hex- 
ane or ethanol, EtOH) were produced by Future Carbon, GmbH. The properties of the 
nanotubes used for the preparation of the skeletons are described in Table 4.1. 


Table 4.1: Properties of the CNTs used for 
the development of buckypapers. 


Property 

Diameter (nm) 10-20 
Length (um) 20-40 
Purity (min wt% C) >98% 


Number of walls 5-15 


The epoxy resin used to infiltrate the CNT skeleton was SR 8100 with SD 8824 hard- 
ener from SICOMIN. This system has low viscosity at room temperature (110 mPa s), 
facilitating the resin flow and impregnation of the skeletons. The gel time of about 
1 h at 40 °C is suitable for a low-pressure infiltration. The recommended cure and 
postcure cycles for this system are 24 h at room temperature and 8 h at 60 °C, 
respectively. 
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4.2.2 Sample manufacturing 


The impregnation of the CNT BP was performed at room temperature using a closed 
mold, as represented in Figure 4.2. For this mold, two spacers are available: one of 
0.5 mm used for one layer composites and the other of 1 mm used for two layer com- 
posites. Three runs were performed using CNT papers with different chemical treat- 
ments, namely hexane or ethanol. 


Figure 4.2: (Left) CNT skeleton positioned in the closed mold, before resin infiltration; (right) resin trap, 
placed before the vacuum pump, used to store eventual excess resin from the infiltration process. 


The procedure used for the infiltration process can be roughly divided into the fol- 
lowing steps: 

Prepare the mold with the application of release agent 

Put spacer in place 

Place the CNT BP in the bottom mold 

Close the mold 

Connect the resin entry and the vacuum pump outlet 

Turn on vacuum pump and verify if mold is completely sealed 

Prepare the resin 

Infiltrate the resin 


po DO» wu. de wm E 


The vacuum level used during the infiltration process was of 0.2 bar; the resin was 
injected in the mold at room temperature and the mold itself was kept at room tem- 
perature. The skeletons were allowed to cure 24 h at room temperature inside the 
mold, followed by demolding and postcure cycle in an electrical oven. 
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4.2.3 Characterization 


CNT papers and composite samples were visual inspected, measured and weighted 
before and after infiltration. Fractured surfaces of samples were observed under a 
high-resolution SEM for the level of impregnation and surface quality assessment, 
as well as for measurements of the nanotube and resin layer thicknesses. 


4.2.4 Results and discussion 
4.2.4.1 Physical characterization 


Figure 4.3 shows an SEM image of the BP surface (a), appearance of a CNT paper 
before (b) and after (c) infiltration. CNT papers have good strength and flexibility to 
allow handling like carbon fiber mat. The SEM image (a) shows that CNT possess 
diameters in the range of 10-20 nm and lengths in the range of 20-40 um. CNT are 
produced solely by catalytic chemical vapor deposition as described elsewhere [1]. 
The number of walls is between 5 and 15. 

Table 4.2 summarizes the main physical characteristics of the CNT papers used 
and of the nanocomposites obtained, namely the thickness and mass values. The 
CNT BP papers used in the manufacturing of samples 1 and 2 are hexane treated 
and the one used for sample 3 is ethanol treated. The thickness of the CNT papers is 
uniform, ranging from 350 to 375 um. 

After resin infiltration, the CNT BP thicknesses increase more than doubling its 
value. This increase is, however, very close to all samples, in the range of 405—430 um. 
The mass of the obtained nanocomposites ranges between 3.11 and 3.98 g and the CNT 
weight content calculated directly (from the variation of mass before and after infiltra- 
tion) is between 11.6% and 15.4% [2]. The higher CNT content was obtained for sample 
1, with hexane-treated nanotubes. From these results we cannot conclude about a rela- 
tion between solvent treatment and CNT content in the composite. Sample 2 is the one 
that has higher mass and consequently lower CNT content. This is due to the higher 
amount of resin that infiltrated the paper. The infiltration process shall be controlled 
in order to achieve highly loaded CNT composites, the reproducibility of the process 
and thus be able to conclude on the effect of the paper solvent treatment on composite 
properties. 
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Figure 4.3: SEM image of the buckypaper surface (a) and the aspect of the buckypaper before (b) 
and after infiltration (c). 


Table 4.2: Summary of the CNT papers and nanocomposite characteristics. 


Sample CNT paper Nanocomposite 


Solvent Thickness (mm) Mass (g) Thickness (mm) Mass (g) CNT wt% 
1 Hexane 0.350 0.4798 0.760 3.11 15.4 


2 Hexane 0.370 0.4616 0.803 3.98 11.6 
Ethanol 0.375 0.4576 0.780 3.27 14.0 


4.2.4.2 Structure observations 


In Figure 4.4, two microscopic pictures are presented, which are obtained under two 
magnifications and from different locations along the thickness of the sample 1. It is 
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(a) 


Figure 4.4: SEM pictures of sample 1 (a) in the bottom side of the composite and (b) in the middle 
of the composite. 


(a) (b) 


Figure 4.5: SEM pictures of sample 2 (a) in the bottom side of the composite and (b) in the middle 
of the composite. 


possible to observe a thick layer of resin (more than 0.2 mm) located in the bottom side 
of the sample and nanotubes well impregnated by resin in the middle of the sample. 

Figure 4.5 presents two microscopic pictures for sample 2. As for sample 1, a 
rich resin layer is observed in the bottom side (lower than 0.2 mm) and upper side 
(lower than 0.100 mm) of the sample, and a good matrix infiltration in the middle 
is witnessed. 

Figure 4.6 shows the pictures obtained for the third run sample. It is observed 
that the BP does not maintain the initial form (appearing slightly curved in the im- 
ages) and their thickness is variable along the diameter. Once again, a thick layer of 
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(a) (b) 


Figure 4.6: Composite SEM pictures of the third run sample: (a) along the thickness of the 
composite and (b) in the middle of the composite. 


resin is covering the paper. In fact, it is possible to see that the thickness of the 
resin layer is much higher then the CNT BP one. Once more, a good matrix infiltra- 
tion in the middle is observed. 

Table 4.3 summarizes the measurements of the top and bottom resin-rich layers 
and CNT paper thicknesses, after the infiltration process, measured from the SEM ob- 
servations. It is observed that the thickness of the paper does not remain unchange- 
able, since it swells for samples 1 and 2 and shrinks considerably for sample 3. 


Table 4.3: Resume of the characteristics of CNT papers and CNT nanocomposites. 


Sample Nanocomposite 
CNT thickness before CNT thickness after Thickness of CNT wt% 
resin infiltration resin infiltration (mm) resin layers (eliminating 
(mm) (top + bottom) resin-rich layers) 
1 0.350 0.428 0.366 34.8 
2 0.370 0.462 0.282 17.4 


3 0.375 0.220 0.624 47.4 


Two possible explanations for this can be prompted: the first one is that the solvent 
treatment has an effect on the resin behavior — contraction or expansion - during the 
cure [3]; the second is that the solvent treatment increases or decreases the infiltration 
of the resin into the paper, either by changing the wettability or the permeability of 
the CNT papers [4]. P.E. Lopes and coauthors [5] studied the wettability (by contact 
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angle measurements) of CNT BP with different treatments. However, the CNT treat- 
ments studied were heat, carboxyl or nonfunctionalized, which are different than the 
ones reported here. 

A rough estimation of the possible CNT content that could be achieved, if the 
resin-rich layers were eliminated, calculated and presented in Table 4.3. The most 
surprising results are the differences between samples 1 and 2, which have the 
same CNT treatment. The second sample has the lower CNT content and confirms 
the high portion of resin in the composite. Once more, we conclude that the infiltra- 
tion process should be well controlled, and a release film must be used in the surfa- 
ces of the CNT papers intended for reducing the thickness of top and bottom layers 
of resin and the permeability of the CNT papers accessed. The resin content is a crit- 
ical aspect of this CNT composites, since the excess of resin can affect negatively 
some properties, namely the thermal and electrical conductivity as well the me- 
chanical properties [6-8]. 

These extrapolated results of CNT content (presented in Table 4.3) are more in 
agreement with those reported in literature [9]. Wang and coworkers have prepared 
SWNT BP/epoxy resin nanocomposites by using very thin papers (0.1-0.5 mm) 
stacked together in a hot press machine. The resin was infiltrated in papers after 
being diluted with acetone due to its high viscosity (2,700 mPa.s). They produced 
nanocomposites with SWNT loadings as high as 39% with nanopore sizes in the BP 
in the range of 100-200 nm. 


4.2.5 Conclusions 


It was the purpose of this experimental work to produce CNT/epoxy resin compo- 
sites by infiltrating a low-viscosity epoxy resin in membranes of nanotube networks 
called BP, supplied by Future Carbon. The CNT content achieved in the composites 
is between 11.696 and 15.496 in weight. In all the SEM observations of the compo- 
sites produced, there is evidence that there is a good level of impregnation, yet the 
resin layer around the CNT BP is very thick. The process must be optimized in order 
to obtain composites without excess of resin in the surfaces, since the resin affects 
negatively some properties of the composites, namely the thermal and electrical 
conductivity as well the mechanical properties. 


4.3 Manufacturing process scale-up and integration 
into CFRP composite structures 
The main purpose of project NACO2 was the scale-up processing of CNT BP and their 


integration on CFRP composites. In order to find the most suitable manufacturing 
method, a number of tests were planned taking into account critical aspects of the 
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manufacturing parameters (such as mold materials, vacuum and pressure levels during 
curing and resin impregnation temperature) [9]. Based on this, several manufacturing 
runs were experimentally performed, aiming to obtain high-quality CNT-CFRP compo- 
sites. From these runs, it was concluded that the use of a steel mold coated with Teflon 
release films, combined with manual resin impregnation at room temperature and cure 
in the autoclave under vacuum and pressure lead to improved materials. At the end of 
the project, a CNT-skeleton composite demonstrator was obtained; herein we present 
the manufacturing details of that demonstrator. 


4.3.1 Materials 


Oxidized CNT skeletons in paper containing binder were used. The epoxy resin sys- 
tem used was: Araldite LY556, Aradur HY906 and Accelerator DYO70, manufactured 
by HUNTSMAN, Switzerland. A Style 450/T300 Plain Weave 140 g/m? carbon fiber 
fabric of Cramer Fabrics, Germany, was used for the manufacturing of the demon- 
strator. Table 4.4 presents the characteristics and suppliers of the raw materials used 
for manufacturing the polymer demonstrator. 


Table 4.4: Resume of the materials used in CNT-polymer demonstrator. 


Materials Main characteristics Supplier 

CNT skeleton Oxidized CNTs, paper form, with binder Future Carbon 

Resin Epoxy system (LY556 + HY906 + DYO70); pot-life: Huntsman 
50—55 h at 25 ?C; viscosity: 1,900-2,100 mPa S 

Carbon fiber CF Style 450/T300 Plain Weave 140 g/m? CC Fabrics 


The CNT skeleton was produced by tape casting and has a binder amount of 2996 in 
weight. The thickness is about 70 um, as measured by optical microscopy and the 
size was 107 x 11 cm’. 


4.3.2 Sample manufacturing 


The demonstrator was produced by hand lay-up followed by vacuum bagging in the 
autoclave and using Teflon as release film. The resin was applied at room temperature 
between all the reinforcement layers, and the CNT BP was applied as the top layer. 
The “striped” surface of the CNT skeleton (BP side that was cast in the release film) 
was faced toward the CFRP material during the impregnation. No resin in the top layer 
(on the BP) was applied. The pressure used in the autoclave was of 3 bar, the vacuum 
level was of 150-250 mbar and the temperature program was 2 h at 120 °C + 2h at 
180 °C. Figure 4.7 shows the different stages of the demonstrator manufacturing. 
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Figure 4.7: Different stages of the manufacturing process of the demonstrator. 


4.3.3 Results and discussion 
4.3.3.1 CNT skeleton handling 


Since the CNT skeleton used had a very low thickness (70 um), it was important 
that the material is very well conditioned. Also, due to its large size, it is crucial to 
be easy to manipulate so the composite manufacturing process is feasible. The 
CNT skeleton used in the demonstrator was packaged between two polymer sub- 
strate foils. The upper substrate foil was very easy to remove; however, the lower 
foil was completely “glued” to the BP. The cutting of the structure was easily per- 
formed with the use of supporting plates to maintain the CNT BP immobile and 
fully supported. Figure 4.8 shows two pictures of the conditioning and handling 
of the skeleton. 


b ——8àÀ—— 


Figure 4.8: Conditioning, handling and aspect of BP used for manufacturing the demonstrator. 
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4.3.3.2 CNT skeleton preparation 


As mentioned earlier, the BP was “glued” to the supporting foil. The procedure for re- 

moving the CNT paper, used in the demonstrator, from the casting substrate is given 

as follows: 

1. cutting of the edges with the aid of supporting plates; 

2. facethe supporting film upward; 

3. release a film tip with the aid of a tweezer; 

4. Separate the film from the BP using a metallic roller, while pulling the film up- 
ward and 

5. continue until both materials are completely separated. 


Figure 4.9 shows the sequence of procedures performed to separate the CNT skeleton 
from the supporting foil. It can be seen that the face of the BP that was “glued” to the 
polymer foil showed a different aspect, where a "striped" pattern is clearly visible. This 
stripped surface was faced toward the CFRP material during the impregnation run. 


Figure 4.9: Procedure used for cutting and removing the supporting film from the buckypaper. 


4.3.3.3 Visual inspections 


Figure 4.10 shows the aspect of the whole sample after curing cycle and demolding, 
as well as samples already cut and ready for characterization tests. Some cracks in 
the BP were visible in the corners of the samples, after the cutting procedure. This 
can be due to either (i) lack of adhesion of the BP to the CFRP and/or (ii) the poor 
impregnation of the BP. 


4.3.3.4 Structure observations 


Figure 4.11 shows the resulting SEM pictures of the demonstrator. The SEM images 
show that, in some extent, the BP is separated from the CFRP laminate. The level of 
impregnation of the fibers seems good, and only a small resin layer (50-90 um) is 
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Figure 4.11: Polymer demonstrator photographs. 


visible on the bottom surface. As previously mentioned, it is visible that the BP is not 
bonded to the laminate. 


4.3.4 Conclusions 


A small demonstrator was manufactured with the main following objectives: 


1. 


2. 
3. 
4 


to evaluate the quality of packaging; 

to assess the hand ability of large and thin BP; 

to determine the effect of the autoclave cycle on the BP's binder; 

to verify the level of wettability and quality of impregnation of the BP of drop of 
resin and 

to demonstrate that the manufacturing of large composites with large and very 
thin CNT layers is feasible. 
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The main conclusions from this activity can be summarized as follows [10, 11]: 

— The CNT paper packaged between two polymer substrate foils and rolled inside 
a metallic container is effectively protected. 

- BP can be easily manipulated when confined between the substrate foils. 
However, the removal of the lower substrate requires additional caution. 

- The most significant change on the BP after autoclave cycle is the appearance 
of a layer of polymer (binder) in the surfaces of the skeleton. 

—- The resin spreads quite easily in the surface (nonstriped) and is able to enter 
into the material. 

- The manufacturing of large composites with large and very thin CNT layers is 
feasible. However, some improvements or additional studies must be done, 
since the BP and the laminate were poorly bonded. 


4.4 Technology demonstration in CFRP mirrors 
for space applications 


The aim of NATAP activity was to identify, develop and test two selected demonstrators 
with functional CNT layers. Due to the high variety of options for adapting material 
properties, an initial trade-off was conducted to identify various space products, where 
the effect and the possibility of implementation would be maximized. Based on previ- 
ous development activities, a large number of applications have been investigated in 
detail with respect to the likelihood to be used in the near future and with potential for 
improvement of weight, cost, performance and resistance to environmental impacts. 
From these, an optical mirror with comparable characteristics of conventional Zerodur- 
based mirrors was selected [12]. The realized demonstrator (28 cm diameter) consists of 
a sandwich structure, composed of CFRP honeycomb. CFRP face sheet and an infil- 
trated CNT BP, which serves as a “grindable surface" for highest accuracy. This surface 
was finally coated with a metal layer (by vapor deposition of aluminum coating and top 
protective SiO/SiO, coating layer) for optical reflectivity in the appropriate wavelength. 


4.4.1 Materials 


FutureCarbon was in charge of providing a dense CNT-based paper or nonwoven struc- 
ture to be used as top carbon layer on the concave (optical) side of the mirror surface. 
From the manufacturing point of view, the CNT paper was the first layer on the mold 
surface. The CNT papers were manufactured using a filtration process, which resulted 
in 4—5 wt% of -COOH groups on the nanotube surface, aiming to provide covalent 
bonding with the resin system. CNT papers with a thickness of approx. 200-300 um 
were manufactured, with diameters in the range of 31-32 cm [13]. 
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TOSEDA supplied cyanate ester-based prepregs, using a pitch carbon fiber fab- 
ric (UD) [Granoc PF(S)-YSH 50A-75 (Mitsui, JP)]. The prepreg sheets used have Fibre 
Volume Fraction (FVF) of about 50% and were manufactured using a solvent-based 
process [14]. 

The honeycomb used was UCF-126-3/8-2.0 RS3 with YSH50 fibers and cyanate 
ester resin, supplied by Ultracore (Table 4.5). 


Table 4.5: Resume of the materials used [15]. 


Materials Main characteristics Supplier 
CNT skeleton -COO groups, 6,96-7,42 mg/cm?, 31-32 cm diameter FutureCarbon 
Resin CE Primaset DT4000 (included as prepreg matrix) Lonza 
Carbon fiber Granoc PF(S)-YSH50A-75 75 g/m? Nippon 
(included as prepreg reinforcement) 
Honeycomb 4 pieces 100 x 200 mm?, UCF-126-3/8-2.0 RS3 Ultracore 
Adhesive film Redux 312L Hexcel 
Prepregs Average carbon fabric volume concentration in prepreg 51 vol% Toseda 


The mold used for manufacturing the CNT-CFRP mirror was made from Invar, with a 
nickel surface treatment. The active surface of the mold, however, was machined 
with sharp tolerances to give the mirror the best surface accuracy possible before the 
polishing operations. 
The prepregs manufacturing process includes the following steps: 
— Impregnation of Granoc fabric with cyanate ester resin diluted to 35% in acetone 
- Air-drying of the impregnated prepregs to remove the solvent (acetone concen- 
tration « 396) 
- Cutting the air dried pre-preg sheets into the requested size followed 
- Marking one corner of the prepreg sheets to keep the same orientation of the 
carbon fabric fibers (warp and weft) 
- Placement of the prepregs between PP foils for protection against bonding dur- 
ing storage and transport. 


4.4.2 Sample manufacturing 


After several manufacturing methodologies used (mainly deriving from the pres- 

ence of cracks on the CNT surface, which are presented below), the procedure for 

manufacturing the mirror demonstrator was defined as follows: 

1. Cutthe CNT BP using a clean sharp tool, with the aid of metal plates, if necessary 
(Figure 4.12). 

2. Dry the CNT BP during 16 h at 110 °C before stacking the materials. 
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Figure 4.12: Cutting the CNT-BP with the proper size. 


If necessary, cut the prepregs in the appropriate dimensions (the warp and weft 
directions of the prepregs are identified). 

Register the weight, area and thickness of the materials. 

Place the Teflon release film in the mold and hold the materials by applying an 
adhesive tape in the corners (Figure 4.13). 

Prepare the vacuum bag. 

Place the CNT-BP in the Teflon release film inside the vacuum bag (Figure 4.14). 
Lay up the two first prepregs layers, taking into account fabric warp and weft 
directions (Figure 4.15). 

Close the vacuum bag and apply vacuum during 10 min with a vacuum level of 
800 mbar with a Teflon release film between the prepreg layers and the vacuum 
bag (Figure 4.16). 


. Repeat the prepreg lay-up procedure until completion. 
. Place the demolding release film as top layer, use breather to eventually absorb 


excess resin and close the vacuum bag. 


. Apply the vacuum pressure (Figure 4.17). If necessary, the materials can remain 


under vacuum overnight. 
Place the materials in the autoclave with the proper curing schedule. 


pm wo 


Figure 4.13: Teflon release film. 
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Figure 4.14: CNT-BP facing the mold. 


Figure 4.15: Prepreg orientations. 


Figure 4.16: Vacuum application. 
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Figure 4.17: Final vacuum-bagged sample, before curing in the autoclave. 


4.4.3 Results and discussion 
4.4.3.1 Requirement definitions 


For specifying the optical mirror properties, the nominal requirements for optical sur- 
faces (not superpolished) were used as baseline. The usual base materials for classical 
optical mirrors are materials like Zerodur (lithium aluminosilicate glass-ceramic pro- 
duced by Schott AG), Astrositall glass, ULE (ultralow expansion titaniasilicate glass), 
fused silica, SiC or nickel Canigen. Table 4.6 lists out the requirements [16]: 


Table 4.6: CNT-CFRP optical mirror requirements. 


Property Requirement 

Mass At least 20% less than conventional mirrors 
Reflectivity R> 0.93 

Surface form tolerance RMS « 70 nm 

Operating temperatures -65 °C/+45 °C 

Nonoperating temperatures +120 °C 

Outgassing According to ECSS [17] 

CME 24 ppm 

CTE <1 x 1076 

Thermal conductivity »0.75 W/mK 


Electrical conductivity »10? S/m 
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4.4.3.2 Mirror design 


The design of the mirror demonstrator was performed taking into account both the re- 
quirements and material properties. The active surface of the mirror is parabolic, circu- 
lar symmetric, in order to ease the machining operations. The design (Figure 4.18) is 
based on a sandwich with a flat skin and a front skin with low curvature, rotationally 
symmetrical with a radius of curvature of 3,360 mm [18]. 


Mirror active surface 
| | CFRP + CNT skin 


CFRP honeycomb 
thickness 15 mm 


| | | CFRP + CNT skin 


Flat backside skin 


Figure 4.18: CNT-CFRP overall mirror design. 


The detailed design also includes the following considerations: 


The reflector sandwich is fully symmetrical in thickness direction. 

Each skin is composed of 16 layers fabric prepreg with surface weight of 75 g/m 
(without resin) together with a CNT layer. 

The honeycomb is made of CFRP using YSH50A carbon fibers. 

The CNT-CFRP skins are further metalized with vacuum deposit aluminum and 
protected with a SiO; layer. 

The adhesive film bonding the CFRP and the honeycomb is based on cyanate 
ester resin. 

The CTE of the combined skin + CNT + adhesive film shall be zero. 

The CTE of the CFRP honeycomb shall be zero. 

A design of the complete sandwich as close as possible to CTE zero is necessary. 


2 


In order to create a very accurate quasi-isotropic lay-up, the following lay-up of the 
CNT-CFRP layers is proposed. All orientations are based on warp fabric direction, 
explained in Figure 4.19. 
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Figure 4.19: CNT-CFRP mirror design: detail of CNT-CFRP skins. 


4.4.3.3 Manufacturing trials 


NATAP sample manufacturing was performed in several phases. These phases and 
main objectives of each one are presented in Figure 4.20 [19]. 

From batch 1 and 2 sample manufacturing, the following conclusions were 
taken [20]: 

- The preimpregnated fiber sheets and BP are easy to manipulate during the 
manufacturing process — including the production of curved samples - there- 
fore, no major risks are attributed to draping for the expected geometries for 
the final applications. 

— Prepregs prepared for batch 2 with 5196 of FVF were found to be a good ap- 
proach for manufacturing the demonstrator, since the resin amount present is 
enough to guarantee a good BP impregnation without the need to add addi- 
tional resins. The CFRP layers were found to be well compacted, reducing the 
risk of having dry CFRP layers. 

— Theoretical calculations also show that these new prepregs result in a final FVF 
of approximately 6596, which is a good balance between high fiber volume frac- 
tion for increased performance and enough resin to guarantee an homogenous 
structure of the CFRP structure (if previous prepregs with higher FVF were 
used, it was expected that composite was significantly dry — with 71-78 vol% 
of fiber). 
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DEVELOPMENT PHASE 1 
Impregnation tests | determine type of CNT-BP functionalization 
Run 1 evaluate if resin in the prepregs was enough to impregnate the CNT-BP 
evaluate raw materials with different characteristics (lower CNT-BP thickness) and new experimental 
conditions (impregnations of CNT-BP with cyanate ester resin before lay-up) 
Run 3 Scale up of run 2 


Run 2 


Manufacturing of three samples with the methodology used in runs 2 and 3 (development phase 1) 


DEVELOPMENT PHASE 2 
evaluate the influence of new conditions of sample manufacturing (no resin added to CNT-BP, application 
of vacuum pressure during lay-up and increase the post cure time) 
Run 2 Scale up of run 1; test the curvability of material 


Run 1 


Manufacturing of three samples with the methodology used in development phase 2 


ATTEMPTS TO REMOVE CRACKS 


Several approaches were performed with the goal to remove/reduce the CNT-BP surface cracks 


Figure 4.20: NATAP project phases and main objectives of each one. 


— The vacuum pressure applied during lay-up is an important step to improve co- 
hesion between layers. 

- Curved samples can be manufactured without significant additional effort, 
however, they were found to be more difficult to demold. It is recommended to 
coat the curved mold with a Teflon release film with extended length (larger 
than the mold) to ease the demolding process; 

— In batch 2, two main problems were detected: samples warping after cure cycle 
(only detected during demolding) and some marks on the surface. Warping can 
be attributed to the formation of tension forces formed during the cure cycle (at 
high temperatures and due to the geometry of the samples) and, therefore, the 
prepregs lay-up sequence may be reviewed for specific geometries; 

— The existence of a surface pattern is should not be a problem, since the marks 
are expected to disappear after polishing and metalization processes. 


Sample characterizations performed throughout the project allowed to detect some 
issues are shown below. The actions taken to avoid/eliminate those issues are also 
presented. 

As referred before, some changes were made in both the characteristics of the raw 
materials used and experimental conditions during sample manufacturing. Figure 4.21 
summarizes the experimental conditions, changes in the methodology used and prob- 
lems encountered in the different stages of NATAP. 
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Project phase 


Development 
phase 1: Run 1 


Issue 


Dry CNT structure, 
only bonded to 
CFRP at the ends 


Characterization 


Detail showing the adherence only in the edges of the CNT 
layer and the carbon fiber; the CNT structure is dry 


Batch 1 CNT BP thickness 
after impregnation 
was higher than 
expected (BP 
expands with - . C NUN i db 
resin addition) Cross-sectional view toss sertiona view performed by 
performed by INEGI AAC (BP-impregnated average 
(nonimpregnated CNT thickness 430 um)/ 
BP average thickness nonimpregnated CNT BP average 
300 um) thickness 
Batch 1 Presence of cracks 
in the CNT BP 
structure 
SEM images of coated surface mirror sample of batch 1 
provided by AAC, showing the existence of cracks 
Batch 1 Glass transition "RI 
temperature lower " 
than expected . 
T, measurements performed by AAC 
Development Lack of resin 


phase 2: Run 1 


BP surface (before 
impregnation) — 
marks are visible 


Cross-sectional SEM picture, 
showing that the structure is dry 
(low resin content); bright surface 
shows that the BP is well 
impregnated 
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(continued) 


Project phase Issue Characterization 


Batch 2 Cracks on CNT-BP 
structure 


«n uam 


Cross-sectional view performed by AAC. Cracks are visible. 


4.4.3.4 Mirror demonstrator manufacturing 


The main steps taken in the demonstrator manufacturing are summarized as fol- 
lows [15]: 

Mold preparation 

Back skin manufacturing 

Front skin manufacturing 

Bonding honeycomb to front skin 

Honeycomb machining 

Back skin bonding to honeycomb (final assembly) 

Final machining 


I4 OX VA qe Q3 dS qs 


Back and front skins were prepared according to the stacking sequence previously 
presented. Figure 4.22 shows the aspect of both CNT-CFRP skins after demolding 
(Figures 4.23—4.25). 

The first assembly step consisted in bonding the honeycomb to the front skin. 

The CNT-CFRP skin was placed in the mold and prepared for the adhesive and hon- 
eycomb settlement. Afterward, the assembly was placed in the oven under vacuum 
pressure and temperature cycles. 
After the first assembly, the honeycomb was machined on the free side in order to 
have a flat surface, for further bonding with the back skin bonding. A flat disk in- 
stalled in a defined position in the turning machine was used for the in-house ma- 
chining process. 

The final assembly consisted in bonding the back skin to the first assembly. 
The back skin, containing adhesive, was placed on the top of the machined honey- 
comb. The entire assembly was then placed in the oven under vacuum pressure and 
temperature cycles. 
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(a) (b) 


Figure 4.22: (a) Mirror back skin and (b) mirror front skin. 


(a) E 


Figure 4.23: (a) Application of adhesive and (b) view of the honeycomb bonded to front skin. 


Figure 4.24: Honeycomb machining to obtain a flat surface. 
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(b) 


Figure 4.25: (a) Assembly after curing with the back skin in the top. (b) Detail of the back surface, 
where no cracks are visible after all the assembly operations. 


4.4.3.5 Mirror polishing and metalization 


One of the project goals was to determine the most suitable polishing process con- 
ditions applicable on CNC polishing machines and allowing to achieve targeted mi- 
croroughness Rq (Sq) of 10-15 nm on the optical mirror surface. The individual 
work steps are summarized as follows [21]: 


4.4.3.6 Polishing 


- Machining outside diameter 
Smoothing an optical surface 
Surface form correction 

— Roughness correction 


4.4.3.7 Metalization 
- By vacuum deposition using e-beam evaporator unit. 


The selected technologies were: 
- Polishing by two-step process using CeO, polishing slurry and polyurethane 
pad 
- Metalization by physical vapor deposition (PVD) technology (e-beam evaporator) 


Figure 4.26 shows the surface appearance of a CNT-CFRP skin after the polishing 
procedure and the surface microroughness of that surface measured using an inter- 
ferometer, showing that the surface root mean square height (Sq) value is 11.82 nm. 
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(a) (b) 


Figure 4.26: (a) Surface appearance of CNT-CFRP skin after polishing and (b) interferometer results 
after polishing. 


The polished CNT surface of the optical mirror was coated by combination of Al 
(100-120 nm) and SiO; (70 nm) protective layer without loss of adhesion. This com- 
bination is a standard type of coating for optical mirrors. SiO; greatly increases du- 
rability of the surface (protection of the Al coating against oxidation). Resulting 
mirror surface has good abrasion resistance, but the top CNT layer of the composite 
is soft (compared to ceramics and glasses); hence, the optical mirror surface can be 
damaged due to damage of the CNT layer. 

The achieved RMS value of 412 nm, which included the outer border, does not 
reach the target value of 70 nm. A reduction of the effective area to exclude the 
edge effect enhances the RMS to -230 nm, which still is not sufficient for the tar- 
geted quality. 

Nevertheless, the quality is well passing the necessary surface quality for the 
backup option of the radiometer example MetOp, where 1,000 nm is targeted. 

Further improvements on RMS quality will be possible when the print-through 
of the backing structure as well as the cracks in the CNT layer can be reduced or, 
respectively, avoided. 

An enhancement of hardness of the top layer could improve the local rough- 
ness, which has been measured to be in the range of 10 nm [22]. 


4.4.3.8 Mechanical characterization 
The mechanical properties have been determined at ambient temperature before 


thermal cycling (RT) and after thermal cycling (RT, TC), and at the maximal (120 °C) 
and minimal (-120 ?C) nonoperational temperatures foreseen for the respective 
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application. The summary of determined mechanical properties for the mirror is 
shown in Table 4.7. The measured density of the material was 1.493 + 0.015 g/ cm’. 


Table 4.7: Summary of mirror mechanical properties [19]. 


Temperature (°C) 3-Point Young’s 4-Point Young’s Interlaminar 
bending modulus bending modulus shear 
strength (GPa) strength (GPa) strength (MPa) 

(MPa) (MPa) 

-120 282.8 + 12.9 24.9 £ 1.4 324.7 + 16.7 27.7 £ 1.5 29.1 £ 1.6 

RT 284.6 + 9.9 31.9 x 1.6 270.5 € 5.5 29.1 € 1.1 26.2 £ 1.4 

+120 262.2 + 5.9 24.9 + 2.7 236.6 + 6.3 32.4 € 1.3 23.6 £1.1 


RT, TC 251.2 + 11.1 28.6 + 1.6 254.5 9.2 27.0 € 0.5 27.121.2 


Summarizing all the above presented results, the mechanical properties of the ma- 
terials are considered suitable for the application with sufficient safety margin. 
4.4.3.9 Electrical characterization 

A total of four measurement points were performed onto the antenna polished and 
metallized surface; three to be placed at the outer circumference of the antenna and 


one into the center (Table 4.8). 


Table 4.8: Radial positions for electrical measurements. 


Position Code Details 

1 FS-R1 to FS-C Front side radial-1 to front side-center 
2 FS-R2 to FS-C Front side radial-2 to front side-center 
3 FS-R3 to FS-C Front side radial-3 to front side-center 
4 FS-R1 to FS-R2 Front side radial-1 to front side-2 

5 FS-R1 to FS-R3 Front side radial-1 to front side-3 

6 FS-R2 to FS-R3 Front side radial-2 to front side-3 


The DC resistance using a Metrohit Multimeter was measured from each of the three 
radial positions (FS-R1 to FS-R3) to the center FS-C (Figures 4.27 and 4.28). 

The results of this test are shown in Table 4.9. The standard deviation of the 
measurements before thermal cycling is very high. A reason for this can be a capac- 
itor effect of the honeycomb structure under the surface. 

The electrical conductivity after the thermal tests was much higher than before 
the thermal cycling. The reason for this can be the many small cracks of the metalized 
surface; therefore, the insulating silica top layer was probably fractured, exposing 
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Figure 4.27: Electrical surface resistance test. 


Figure 4.28: Surface crack in mirror sample after thermal cycling (left), Cross section of mirror 


sample (right). 


Table 4.9: Summary of electrical measurement results. 


Code Before thermal cycling After thermal cycling 
Average (Q) Standard Conductivity Average (Q) Standard 
deviation (S/m) deviation 
FS-R1 to FS-C 9,166.7 1040.8 9,92E-04 21.7 7.2 
FS-R2 to FS-C 6,166.7 1040.8 1.47E-03 13.0 4.2 
FS-R3 to FS-C 3,600.0 1216.6 2.53E-03 21.3 4.6 
FS-R1 to FS-R2 2,066.7 378.6 2.42E-03 27.1 5.8 
FS-R1 to FS-R3 2,000.0 435.9 2.50E-03 23.3 2.5 
FS-R2 to FS-R3 2,566.7 907.4 1.95E-03 22.7 3.0 
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either the aluminum or the CNT layer, and we could measure the electrically conduc- 
tive surface through the top layer. 

The requirement of 10? S/m, that is, *metallic" electrical conductivity, could 
not be reached. The as-received surface conductivity was 4.28 S/m for the CNT side 
and 2.59 S/m for the CFRP side. Careful roughening the surface with Grit 1200 sand- 
paper improved the surface conductivity substantially to 4.28 S/m for the CNT side, 
and 16.18 S/m for the CFRP side where carbon fibers were apparently contacted. 


4.4.3.10 Thermo-optical properties 
Thermo-optical properties were tested at ZAE Würzburg, each for the CNT and for 


the CFRP side. The thermo-optical properties do not differ from CFRP to CNT side 
(Table 4.10). 


Table 4.10: Summary of mirror absorbance values. 


Sample Visual absorbance Solar absorbance UV absorbance 
CFRP face 0.92 + 0.02 0.91 + 0.02 0.93 + 0.02 


CNT face 0.94 + 0.02 0.94 + 0.02 0.94 + 0.02 


4.4.3.11 Structure observations 


The mirror samples consisted of about 1.1-1.2 mm CFRP, and a CNT layer of about 
350 um. The surface shows cracks after thermal cycling. 

Table 4.11 summarizes the main results achieved with the CNT-CFRP mirror, 
highlighting that mass, CTE and thermal conductivity properties achieved are in 
high-performance levels; however, specific heat, surface accuracy and reflectivity 
values are below the requirement definitions. 


4.4.4 Conclusions 


The selected demonstrator — a CFRP-based mirror — was manufactured using two CFRP 
skins based on cyanate ester prepregs, with BP on the top surface. The CNT structure 
was incorporated as top layer of the skin during prepregs lay-up and then impregnated 
with the resin from prepregs. These CFRP skins were assembled with a CFRP honey- 
comb for obtaining the final part, further polished and metalized. Mirror mass was sig- 
nificantly improved at very good thermoelastic distortion values, but the CNT surfaces 
have many cracks, which further propagated after thermal cycling. 
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Table 4.11: Optical mirror demonstrator — major results. 


Property Results Performance 
level 
Mass More than 2096 mass saving in relation to a conventional ++ 
optical mirror can be achieved. The mirror demonstrator’s 
mass was only 270 g. 
CTE -0.706 x10 $/K in X/Y-direction (requirement:<+1 x 1076 /K) ++ 
Thermal 0.84 W/(m*K), better than required (20.75 W/(m*K)) + 


conductivity 


Specific heat 


0.73 J/(kg*°C) acceptable in relation to Zerodur 0.85 J/(kg*°C) 


Surface accuracy 


0.135 um has been reached, not far away from requirements 
(«0.1 um) 


Surface quality 


Not reached - presence of many cracks 


Reflectivity 0.7 (after cycling), required value (20.93) not reached due to -- 
cracks in surface 

Temperature -120 °C > +100 °C 

range 


Several development steps were performed, aiming to eliminate or reduce the 
formation of cracks, such as the study of the influence of the lay-up, curing condi- 
tions and mold conditions, but none of the proposed approaches was effective. 
After all physical parameters being tested, the chemical influence of the resin was 
investigated. The cracks were eliminated by replacing the CTE for an epoxy resin in 
the CNT layer. 

In fiber-reinforced composites, the CTE is affected primarily by the carbon fi- 
bers, while on the CNT layer, the resin CTE has the dominant effect. CTE dispropor- 
tion can generate fatal stresses. 


4.5 Conclusions 


In this chapter, we presented a brief overview of the activities performed under three 
ESA projects, as well as the main results, targeting the use of CNT structures in the 
space sector. The first project was intended to prove the feasibility of CNT structures — 
BP - to be infiltrated with an organic matrix, potentially as an alternative to poorly 
dispersed and high viscosity nanotube-filled resin systems. Focus was given on as- 
sessing the morphology of the developed materials, by using high-resolution 
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microscopic analyses. The results show the proposed processing approach as capa- 
ble of producing well-dispersed nanocomposites with high CNT loading (more than 
15 wt% CNT by composite weight), which are important for developing high- 
performance structures based on CNT with good thermal and electrical conductivity. 

The main objective of this activity was to scale up the material developments 
and thereby increase the Technology readiness level (TRL) to 2 for CNT skeleton- 
based polymer composites. The developed materials are used for manufacturing 
and testing of demonstrators for typical space applications. Several experimental 
trials were performed, trying to optimize process conditions such as mold tools, cur- 
ing conditions or the use of different methods for BP impregnation. In a first step, 
the manufacturing process for obtaining CNT-skeletons with well-defined and repro- 
ducible characteristics was followed. In a second step, the composite manufacturing 
processes were further developed on sample level and in parallel the dimensions of 
the CNT-skeletons were scaled up to 0.1 m x 1 m with thickness ranging from few tens 
of micrometers up to few centimeters. In the third step, the CNT-reinforced composite 
manufacturing was stabilized on scaled-up dimensions, suitable for manufacturing 
technological demonstrators. It was possible to obtain CFRP/CNT skeleton composites 
based on small papers in a controlled manner, with high CNT content, no defects and 
voids and without excess of resin in top, interlayer and bottom surfaces. 

The final activity targeted a space application, and a demonstrator was manufac- 
tured and tested during the project. The selected application was an optical mirror, in 
which the CNT structures were intended to improve the surface characteristics of the 
material, facilitating the polishing and metalization steps. Extensive mechanical, 
thermal and optical characterizations were performed at sample and demonstrator 
level. The CNT-CFRP mirror resulted in more than 2096 of mass savings, when com- 
pared to a conventional mirror, with a final mass of only 270 g. The measured surface 
accuracy was of 0.135 uum, not much distant of requirements (< 0.1 um). The main 
issue of the optical mirror is related to the presence of several cracks on the CNT 
layer, limiting the surface quality and the reflectivity values obtained. Nevertheless, 
the principle of producing a highly accurate optical mirror based on a CFRP/CNT is 
rated as promising, with good thermoelastic distortion values. It is estimated that a 
TRL level between 3 and 4 has been reached. 
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Emmanuel N. Skountzos and Vlasis G. Mavrantzas 
5 Molecular dynamics simulations of 
graphene-based polymer nanocomposites 


5.1 Introduction 
5.1.1 Carbon-based materials — the case of graphene 


During the last decades, carbon has attracted a great deal of scientific and indus- 
trial attention due to the discovery of several allotropes (graphene, fullerene, car- 
bon nanotubes [CNTs]), which are characterized by unprecedented physical and 
chemical properties such as high mechanical strength, extremely high electrical 
and thermal conductivity, high optical transparency and excellent gas barrier prop- 
erties. Because of these unique features, carbon materials are widely used as nano- 
fillers for the fabrication of composite materials with applications in several fields 
(biology, energy storage, transport and aviation, optoelectronics, pharmaceutics, 
medicine and many others). Fullerenes, for example, are used as electron acceptors 
for the fabrication of organic solar cells based on semiconducting polymers such as 
poly-3-hexyl thiophene, which has significantly improved the efficiency of the cor- 
responding devices, reaching values comparable to those of inorganic solar cells. 
Graphene-based biosensors characterized by remarkable detection efficiency to- 
ward certain target molecules have also been manufactured by exploiting the excel- 
lent electrical and optical properties of graphene even when present at very low 
concentrations. The very high specific surface area of CNTs combined with their 
low electrical resistance and high charge transport capability have allowed the pro- 
duction of high-performance supercapacitors with increased energy storage, power- 
delivery capabilities and rather long-life cycle compared to conventional batteries. 

From the broad family of the carbon-based nanoparticles (NPs), the one that 
has attracted the most significant attention is graphene, a one-atom-thick planar 
sheet of sp2-bonded carbon atoms densely packed in a honeycomb crystal lattice. 
Graphene is the basic building block for graphitic materials of all other dimension- 
alities (Figure 5.1), because it can be wrapped up into OD fullerenes, rolled into 1D 
nanotubes or stacked into 3D graphite [1]. 
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Figure 5.1: Graphene, the building block of all graphitic forms. 


Until 2004, the production of single-layer (SL) graphene sheets (GS) was con- 
sidered impossible, since theoretical calculations had predicted that GS would be 
completely unstable. The successful isolation of SL, defect-free GS using the "scotch 
tape method" [2] led to a rapid expansion of the field. The discovery was made by 
Prof. Sir Kostya Novoselov and Prof. Sir Andre Geim, who were awarded the Nobel 
Prize in Physics in 2010. 

The impact of graphene on the scientific community has been enormous. 
Thousands of scientific articles (experimental and theoretical) are published 
every year and hundreds of these are heavily cited, while the paper describing 
the simple method to produce perfect GS has been characterized as “one of the 
most cited recent papers in the field of Physics" according to the ISI citation 
index. Apart from the great scientific interest in this exotic material, a significant im- 
pact on the largest technological industries of the world has also emerged after its 
discovery. Every year, thousands of patents are filed around the world on the design 
of new materials and new devices with important implications for practically all as- 
pects of modern life (medicine, transportation, electronics, etc.), which explains why 
graphene has been characterized as the material of the future. 


5.1.2 Graphene and graphene-based polymer nanocomposites 


Numerous excellent reviews have been written in the last few years on the produc- 
tion, properties and applications of graphene [3-7] and graphene-based nanocom- 
posites [8-13]; hence, here we will limit ourselves only to a brief overview of its 
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unique structural and physicochemical properties, and for more details we refer the 
interested reader to the original articles. 


Production of graphene 


After a long and tenacious series of unsuccessful attempts to produce SL graphene, 
the publication of a simple method was presented in 2004 [2] known today as the 
“scotch tape method.” By repeatedly cleaving a graphite crystal flake with an adhe- 
sive tape to its limit, and then transferring the thinned-down graphite onto an oxi- 
dized silicon wafer with the appropriate color, a 2D carbon lattice is produced. 
In general, production methods of graphene fall into two main categories: 
— Bottom-up methods 
— Top-down methods 


In bottom-up processes, graphene is synthesized by a variety of methods such as 
1. Epitaxial growth on metal carbides [14-16] 

2. Chemical vapor deposition (CVD) [17-18] 

3. Unzipping CNTs [19-20] 


In top-down processes, it is synthesized from a bulk material (e.g., graphite), which 
is broken down into smaller pieces using mechanical, chemical or other forms of 
energy. Typical examples of the top-down production processes include 

1. the micromechanical exfoliation of graphite [2]; 

2. the direct sonication of graphite [21-22]; 

3. the chemical reduction of organically treated graphite oxide [23] and 

4. the thermal exfoliation/reduction of graphite oxide [24]. 


In principle, with bottom-up approaches, large-sized, defect-free monolayer gra- 
phene is produced, a perfect material for the subsequent studies. The only draw- 
back with these methods is that they afford only the production of tiny amounts of 
graphene. Top-down processes, on the other hand, are suitable for the large-scale 
production of graphene as (for example) required for the fabrication of graphene- 
based polymer nanocomposites. 


Properties of graphene 


The rapid adoption of graphene as the material of interest lies primarily in the ex- 
cellent spectrum of properties (mechanical, electrical, thermal, optical, etc.) charac- 
terizing monolayer, few-layer graphene and graphene oxide (GO). The goal of many 
research efforts nowadays is to exploit these extraordinary properties for applications 
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in nanotechnology by fabricating materials with improved mechanical, electrical, 
thermal and optical performance. 

Carbon-based nanomaterials such as graphite, diamond and CNTs have their 
own record in terms of mechanical strength, hardness or Young's modulus. The 
newest member of the family, graphene, is no exception despite that its mechanical 
behavior has not been investigated as much as its electronic and optical properties. 
The reported stiffness of about 300-400 N/m (with a breaking strength of about 
42 N/m) represents the intrinsic strength of a defect-free sheet [25] while estimates 
of the Young's modulus are on the order of 0.5-1.0 TPa [25]. Interestingly, and de- 
spite their defects, suspended GO sheets retain to a large extent their mechanical 
performance, characterized by a Young's modulus of 0.25 TPa [26]. These features, 
combined with the relatively low cost for the production of thin graphite and the 
ease of processes for blending GO into matrices [11], render these materials as ideal 
candidates for electrical and mechanical reinforcement [27, 28]. 

Regarding electrical properties of graphene, we mention that its electrical resis- 
tivity at room temperature is about 1 uQ cm; as a result, graphene is about 35% less 
resistant than silver, the lowest resistivity material known today at room tempera- 
ture. In semiconductors, a different measure is used to quantify electronic motion 
known as mobility. Mobility is often expressed as the conductivity of the material 
per electronic charge carrier. This implies that high mobility is advantageous also 
for chemical or biochemical sensing applications in which a charge signal (e.g., a 
molecule adsorbed on a device) is translated into an electrical signal, thanks to the 
changing conductivity of the device. Thermal vibrations of atoms set the upper 
limit to electron mobility in graphene, which is -200,000 cm?/V s [29-30] at room 
temperature, which should be compared to ~100,000 cm?/V s [31] in CNTs. 

Chen et al. [32-33] showed that although the room temperature limit of mobility 
in graphene can be as high as 200,000 cm?/V s, in present-day samples the actual 
mobility is lower (around 10,000 cm?/V s) leaving significant room for improve- 
ment. Because graphene is only one-atom thick, current samples must be supported 
by a substrate, typically silicon dioxide. Trapped electrical charges in silicon diox- 
ide (a sort of atomic-scale dirt) can interact with electrons in graphene, which can 
cause a reduction in mobility. 

In addition to being considered as a promising material for applications in 
chemical and biochemical sensing, its low resistivity and extremely thin nature ren- 
der graphene a very promising material also for use in thin, mechanically tough, 
electrically conducting, transparent films that are needed in many applications in 
electronics (ranging from touch screens to photovoltaic cells). 

Apart from its excellent mechanical and electrical properties, graphene exhibits 
amazing performance as a thermal conductor. Balandin et al. [34] reported values 
of thermal conductivity for SL graphene at room temperature in the range of 
4,840-5,300 W/m K. These extremely high values of thermal conductivity suggest 
that graphene can outperform CNTs in heat conduction [35-36]. However, when 
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graphene is in contact with a substrate, its thermal transport properties can be sig- 
nificantly affected. Seol et al. [37] showed experimentally that the value of thermal 
conductivity of monolayer graphene exfoliated on a silicon dioxide support is still 
as high as about 600 W/m K near room temperature, exceeding those of metals 
such as copper. It is lower, however, than that of suspended graphene because of 
phonons leaking across the graphene-support interface and strong interface scat- 
tering of flexural modes, which make a large contribution to thermal conductivity 
for suspended graphene. In general, the superb thermal conduction properties of 
graphene have established it as an excellent material for thermal management. 


Production of graphene-based polymer nanocomposites 


In recent years, a variety of processing methods have been proposed for dispersing GS 
into polymer matrices. Many of these procedures are like those used for other nano- 
composite systems (e.g., CNT/polymer nanocomposites) [38], while others apply only 
to graphene-based polymer nanocomposites. 

A crucial step in the production of any polymer nanocomposite is the disper- 
sion of the nanofiller. A well-dispersed state ensures a maximized reinforced sur- 
face area, which directly affects all properties of the nanocomposite. Efforts are 
therefore focused on achieving a well-dispersed, homogeneous system by develop- 
ing either covalent or noncovalent functionalization of the filler surface, an issue 
that will be discussed in some more detail in the following sections. 

Most polymer/graphene composites are produced today with one of the follow- 
ing three strategies: (1) solvent processing, (2) in situ polymerization and (3) melt 
processing with each one of them having its own advantages and disadvantages. In 
the solvent processing method, GS are initially dispersed in a suitable solvent, a 
procedure typically assisted by ultrasonication [11, 39]. Then, the polymer is added 
in the solvent/graphene blend, and the solvent is finally removed by evaporation or 
distillation. It is a rather simple method, used widely to prepare polymer/graphene 
composites. Its most important drawback is that common organic solvents adsorb 
on GS stronger than most of the polymers. 

In the in situ polymerization method, GS are mixed with the targeted mono- 
mers, and by adjusting parameters (such as temperature and pressure), the poly- 
merization reaction proceeds [40-41]. The advantages of the method are twofold: 
(a) it provides a strong interaction between the polymer matrix and the surface of 
GS and (b) it leads to highly homogeneous dispersions. However, an increase in the 
viscosity of the blend is usually recorded as a side effect, which affects processabil- 
ity (especially at high GS loadings). 

Melt processing is commercially the most attractive method to produce gra- 
phene-based polymer composites. It involves the direct inclusion of GS into the 
melted polymer using a twin-screw extruder by suitably adjusting parameters such 
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as screw speed, temperature and time [42-43]. Drawbacks of the method include 
the low density of thermally exfoliated graphene that makes extruder feeding a 
troublesome task, and the lower degree of dispersion achieved compared to solvent 
blending. Reduced degree of dispersion typically results in poorer mechanical, elec- 
trical and thermal properties. 


Properties of graphene-based polymer nanocomposites 


In the modern literature, significant improvements in the mechanical, electrical, ther- 
mal and barrier properties of graphene-based polymer nanocomposites have been re- 
ported, as summarized in several reviews [8-10, 12, 44]. Fang et al. [40], for example, 
have reported an increase in the tensile strength by 7096 and in the Young's modulus 
by 5796 for graphene-polystyrene nanocomposites with polystyrene (PS) chains 
grafted onto GS by atomic transfer radical polymerization. The experimental studies of 
Ramanathan et al. [39] and Li and McKenna [45] have shown that behind the extraordi- 
nary mechanical properties and the increase in the glass transition temperature (T) of 
poly(methyl methacrylate) (PMMA) nanocomposites filled with functionalized gra- 
phene sheets (FGS) are the enhanced interfacial interactions with PMMA chains as 
driven by oxygen functionalities across the surface of graphene. FGS contain pendant 
epoxy, hydroxy and/or carboxy groups on their surface, which may form hydrogen 
bonds with the ester branches of PMMA. However, Liao and coworkers [46-47] have 
argued that an 80% increase in the PMMA modulus [39] at only 1 wt. 96 loading of the 
nanocomposite in FGS, and an increase in T, by 29 °C [39] at only 0.05 wt.% loading 
seem unrealistically high. Thus, Liao et al. repeated the experiments carried out by 
Ramanathan et al. [39] and found an increase of only 2596 in the Young's modulus 
measured [47] and no change in the T; [46]. These very different observations were 
attributed to the experimental procedure followed by Ramanathan et al. [39]. In 
an earlier study [A8], significant improvements in the Young's modulus and in ul- 
timate tensile strength had been reported for poly(vinyl alcohol) (PVA) samples 
enhanced with GO sheets functionalized with PVA chains (PVA chains had been 
grafted onto the GO surface). The reported enhancement in the mechanical perfor- 
mance reached almost 60% for GO loadings below 0.3 vol.% [48]. For the same 
nanocomposite, Zhao et al. [49] have reported an improvement of ~150% in the 
tensile strength and an order of magnitude increase in the Young's modulus for 
only a 1.8 vol.96 loading in graphene. 

As already mentioned, graphene-based materials are very promising for the de- 
velopment of new devices with applications in electronics, owing to their high degree 
of electrical conductivity. In the past, several carbon-based NPs (i.e., carbon filler, 
carbon nanofibers, expanded graphite, etc.) have been exploited for the production 
of electrically conductive composites; however, the key advantage of graphene is 
that the insulator-to-conductor transition (known as the percolation threshold) can 
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be achieved at significantly lower loadings. Production of electrically conductive pol- 
ymers has been reported [50-51] upon successful GS dispersion in the host polymer 
matrix. Stankovich et al. [11, 52] determined the percolation threshold for a polysty- 
rene solvent blended with GO to be 0.1 vol.96, perhaps the lowest percolation thresh- 
old ever reported. Eda and Chhowalla [53] studied the electrical properties of 
solution-processed, semiconducting thin films consisting of FGS as the filler and 
polystyrene as the host material and found that upon increasing the average size of 
FGS significantly enhanced carrier mobility and thus device performance. This study 
demonstrated how a commodity plastic can be used to develop low-cost, macroscale 
thin-film electronics. 


Simulations 


In addition to experimental efforts, theoretical and computational works have ad- 
dressed several aspects of the structure-property- processing relationship in gra- 
phene or graphene-based nanocomposite materials. Several simulation techniques 
and approaches have been employed, extending from the quantum level to the at- 
omistic to the mesoscopic and finally to the macroscopic. The findings of these the- 
oretical studies have significantly improved our understanding of the microscopic 
mechanisms and interactions governing the macroscopically exhibited properties 
of these new classes of materials. 

Using classical ab initio calculations, Van Lier et al. [54] and Liu et al. [55] re- 
ported values of graphene's Young's modulus equal to 1.11 and 1.05 TPa, respec- 
tively, which are in reasonable agreement with the experimentally measured ones 
[25]. It is worth mentioning that the computational values were reported in 2000 
[54] and 2007 [55], respectively, while the experimental one in 2008 [25]. In a recent 
study that combined density functional theory calculations and classical molecular 
dynamics (MD) simulations, Kalosakas and coworkers [56] proposed a new force 
field specifically for graphene that takes into account only bond stretching (de- 
scribed by a Morse-style potential) and bond bending (described by a nonlinear 
function containing quadratic and cubic terms) interactions between carbon atoms. 
The new potential was employed in simulations with model graphene systems sub- 
ject to uniaxial tension or to hydrostatic compression, yielding a value of 0.95 TPa 
[56] for graphene's Young's modulus, which is also consistent with the one mea- 
sured experimentally [25]. The dependence of graphene's Young's modulus on tem- 
perature and size of GS was systematically studied by Jiang et al. [57], through MD 
deformation simulations of SL graphene using progressively larger GS, and an in- 
crease in the Young's modulus was observed with increasing GS size. The plateau 
value was reached for a GS size equal to 25 À x 25 À beyond which no further in- 
crease was recorded [57]. The MD deformation experiments were carried out at tem- 
peratures ranging from 100 to 600 K; a slight increase in the value of the Young's 


122 —— Emmanuel N. Skountzos and Vlasis G. Mavrantzas 


modulus (within the statistical error) was monitored for temperatures up to 500 K, 
followed by a rapid decrease at higher temperatures. 

Ab initio studies of suspended GS are not restricted solely to the estimation of 
the mechanical properties. Excellent articles have been reported in the literature ad- 
dressing atom-atom interactions between graphene and popular substrates (e.g., SiC 
[58] kat MoS, [59]) that are extensively used for the fabrication of SL graphene. Also 
important are computational and simulation studies of the mechanical, thermal, bar- 
rier and electronic properties of polymers filled with GS. Issues addressed here in- 
clude microscopic structure, chain conformation and local and terminal dynamics of 
polymer matrix chains in the presence of graphene. Earlier works focused on the 
study of the interfacial behavior of graphene-based polymer composites. Awasthi and 
coworkers [60] carried out atomistic MD simulations with the consistent valence 
force field to study nanoscale load transfer between polyethylene (PE) and GS and 
characterize the force-separation behavior between CNTs and a polymer matrix. 
Separation studies were conducted for opening and sliding modes, and cohesive 
zone parameters (such as the peak traction and the energy of separation for each 
mode) were evaluated as a first step toward the development of continuum length- 
scale micromechanical models for tracking the overall material response by incorpo- 
rating information about the underlying interfacial interactions. MD simulations have 
also been employed by Li et al. [61] in their study on the effect of the shape of car- 
bon-based NPs on the viscoelastic properties of a PE matrix. They found that it is the 
surface-to-volume ratio of the NPs that plays the most important role in the struc- 
tural, dynamical and viscous properties [61]. More recent MD simulations [62-63] of 
PMMA/graphene nanocomposite showed strong adhesion of PMMA chains (espe- 
cially of the side groups) on graphene, and considerably slower segmental and chain 
mobility in the interfacial area. The MD simulations suggest that local mass density, 
segmental dynamics and chain terminal relaxation all differ from the bulk behavior 
up to distances equal to several nanometers from the GS surface. Very similar results 
have been reported for a different matrix, PE [64], demonstrating large density inho- 
mogeneities due to strong PE chain adsorption on the surfaces of GS, exactly as was 
reported in the case of PMMA/graphene nanocomposites [62]. Close to graphene, PE 
chains prefer to stand parallel to the graphene surface [64], and all polymer confor- 
mational and dynamic properties [64—65] are significantly affected: (a) the size of 
polymer chains (as measured by their radius of gyration) increases and (b) their dy- 
namics change dramatically because their orientational relaxation time increases al- 
most by one order of magnitude compared to the bulk value. In a very recent work 
[66], the effect of GS on the crystallization process of PE, polyvinylidene fluoride 
(PVDF) and PS oligomers was examined. It was reported that GS tend to act as nucle- 
ation sites for the crystallization of PE and PVDF but not for PS, which remains al- 
most amorphous [66]. It was also reported that at high temperatures (e.g., close to 
600 K), the crystalline structure of PE is destroyed, a result that is in accordance with 
the recent MD study of Gulde et al. [67]. 
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Simulations have also addressed polymer nanocomposites enhanced not with 
pristine GS but with GO. GO is graphene-bearing epoxy, hydroxy and/or carboxy 
groups on its surface. When GO is used as a nanofiller of a polar (e.g., acrylic) poly- 
mer, the interactions between polymer atoms and graphene are intensified due to 
strong attractive forces that develop between polymer and GO oxygen atoms. Lv 
et al. [68] examined two different polymers as the host matrix, PMMA and PE, and 
concluded that with increasing concentration in carboxyl content, the interaction 
energy between polymer chains and GO sheets decreased (it became more attrac- 
tive), followed by a significant increase in the value of the shear stress. However, 
an upper bound in the concentration was found, beyond which no further change 
in the values of these two properties was observed. This was explained as a satura- 
tion effect: high concentrations of functional groups strengthen the interactions be- 
tween the surface functionalization groups themselves, thus no space is left for 
interactions with the surrounding polymer chains. Karatasos and Kritikos [69] have 
reported a strong increase (by 38 °C) of the T, of GO-based poly(acrylic acid) (PAA) 
nanocomposites compared to pure polymer (PAA). This strong increase was ex- 
plained by the strong adsorption of PAA chains onto the surface of GO sheets facili- 
tated by the hydrogen bonds that develop between the hydroxyl groups of GO and 
the oxygen atoms of PAA branches. Earlier, Xue et al. [70] had studied the T, of 
PMMA matrices enhanced either with pristine GS or with GS modified with -COOH 
and -NH, groups. By employing classical MD simulations, they observed a 30 °C in- 
crease in the T, of the PMMA/GS system, whereas for the nanocomposites with the 
modified GS, the shift was higher (40 °C) [70]. More recently, Azimi et al. [71] 
showed that the dynamics of a polar (e.g., PVA) matrix is affected more by the pres- 
ence of GO than by the presence of GS, whereas for an apolar polymer matrix (e.g., 
poly(propylene)) the effect of the two types of graphene is almost the same. That 
less polymer is adsorbed on GO than on pristine graphene can be explained by the 
roughness of GO particles due to OH- and -O- groups on their surface and agrees 
with a recent detailed MD study by Skountzos et al. [72]. However, the strength of 
interactions between oxygen atoms of the polar polymer and of GO particles is so 
strong that despite the smaller adsorbed amount on GO (in comparison to pristine 
GS), the dynamics of the PVA/GO system is considerably slower than the dynamics 
of the PS/GS nanocomposite. 

In the last years, significant progress has been made in predicting the unique 
mechanical properties of graphene-based polymer nanocomposites through de- 
tailed atomistic-level simulations and understanding the underlying molecular 
mechanisms behind these properties [55, 72-76]. A typical example is the atomistic 
simulation work of Skountzos et al. [72] on the effect of pristine graphene and GO 
on the structure, conformation and mechanical properties of a syndiotactic PMMA 
(sPMMA) matrix. The atomistic simulations predicted a significant enhancement of 
all elastic constants (Young’s, shear and bulk moduli and Poisson’s ratio), espe- 
cially for the PMMA/GO nanocomposites, which was attributed to the hydrogen 
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bonds that develop between the oxygen atoms at the branches of PMMA chains and 
the OH- and -O- groups on the surface of GO particles. Wang et al. [76] have studied 
the effect of surface functionalization and graphene size on the interfacial load 
transfer in graphene-PE nanocomposites (by employing the ab initio polymer- 
consistent force field, PCFF) and showed that oxygen-FGS lead to larger interfacial 
shear force (than hydrogen-functionalized or pristine ones) during the pull-out pro- 
cess. Increasing the oxygen coverage and graphene size enhanced the interfacial 
shear force, but further increasing the oxygen coverage to about 796 led to a satura- 
tion in the interfacial shear force. Liu et al. [75] studied the effect of shape of car- 
bon-based NPs on the toughening efficiency of a PE matrix by applying pure tensile 
strain experiments in MD simulations with several model systems. At the same wt. 
96 loading in either CNTs, or fullerenes or GS, they found that the highest toughen- 
ing was observed for the PE/graphene system; this should be related to the in- 
creased surface area that graphene offers for PE chain adsorption. This effect was 
examined in both the rubbery and glassy states (i.e., above and below the glass 
transition temperature of PE) and it was found that the nanocomposite in the rub- 
bery state was tougher than the corresponding glassy material [75]. More recently, 
in a systematic study, Lin et al. [73] investigated how temperature and GS loading 
affect the mechanical performance of model PMMA/GS nanocomposites. The gen- 
eral trend was that Young's and shear moduli increase with increasing GS loading 
but decrease with increasing temperature while the temperature dependence of the 
two moduli is stronger in the nanocomposites with the higher concentration in GS. 
Comparison of their simulation results with the equations based on the simple rule 
of mixing for the mechanical properties of microfiber-reinforced composites re- 
vealed that the latter do not apply for graphene-reinforced nanocomposites; thus, 
further research is needed to address this issue [73]. 


Graphene dispersion 


Clearly, among all NPs considered as polymer reinforcing agents, graphene is the 
most promising due to its outstanding features. However, because of its high spe- 
cific surface area and tendency for self-adhesion (driven by the very strong van der 
Waals forces and rt-7 interactions), GS tend to agglomerate in the form of multi- 
layer graphitic structures. Significant research has thus been undertaken in the last 
years to understand and improve GS dispersion, and thus realize its unique proper- 
ties in practice. The degree of GS dispersion in nanocomposites affects primarily 
their mechanical performance. Indeed, several studies have shown that the me- 
chanical properties of graphene-based polymer nanocomposites increase as their 
content in SL graphene increases [39, 77-78]. Motivated by these findings, in the 
next paragraphs of this chapter we provide a brief overview of the methods employed 
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nowadays, typically based on covalent or noncovalent graphene functionalization, to 
achieve homogeneous GS dispersions [79-80]. 

Techniques based on covalent functionalization make use of small molecules 
or entire macromolecules that are attached on the surface of GS through covalent 
bonding. Then, owing to strong chemical interactions developing between the 
grafted groups and the GS, very homogeneous and stable dispersions can be ob- 
tained. However, covalent functionalization disrupts the perfect crystal structure of 
graphene, which can affect many of its outstanding properties, especially the electrical 
ones. An attractive alternative is noncovalent functionalization. This is based on the 
use of functional groups that are attached on the surface of graphene through favor- 
able n-n interactions without disturbing its electronic network. Noncovalent methods 
are by default nondestructive; however, the forces that develop between the wrapping 
molecules and the graphene surface are much weaker compared to those in the case of 
covalent functionalization; this can limit the degree of dispersion eventually achieved. 

Covalent functionalization of graphene has been exploited quite widely in the liter- 
ature for the improvement of the mechanical properties of graphene-based nanocom- 
posites [81-85]. Using mostly peripheral ester linkages that are extensively used in 
tissue engineering applications, Sayyar et al. [84] managed to successfully link polycar- 
bonate to GO. The fabricated nanocomposite had a very homogeneous and stable GS 
dispersion and exhibited several good properties: improved Young's modulus, im- 
proved tensile strength and high electrical conductivity (-14 orders of magnitude com- 
pared to the conductivity of the pure bulk polymer). In an earlier study, Cheng and 
coworkers [81] had produced PVA nanocomposites by incorporating PVA-grafted GO 
fillers, which resulted in an enhancement of Young's modulus by 15096 and of tensile 
strength by 88%, for 1 wt.% loading of the nanocomposite in GO. 

The tendency of GS to agglomerate and form graphitic structures is witnessed 
most easily in atomistic or coarse-grained simulations. Park and Aluru [86] have re- 
ported GS self-assembly in atomistic MD simulations of GS in water solutions, but this 
was probably an expected result due to hydrophobic nature of graphene. In a later 
study, Li and coworkers [61] found that several carbon allotropes (buckyballs, graphe- 
nes, etc.) agglomerate in a PE matrix, thus lowering the degree of their interaction (con- 
tacts) with the host matrix. Ju et al. [87] studied the degree of GS miscibility in PMMA 
as a function of NP volume fraction with the help of dissipative particle dynamics 
(DPD) simulations using different values of the DPD-repulsive interaction parameter to 
effectively account for surface functionalization. GS agglomeration has also been re- 
ported by Karatasos et al. [88-89] in all-atom MD simulations of linear and hyper- 
branched polymer matrices. Guo et al. [90] have studied the effect of GS agglomeration 
on the mechanical properties of polymer nanocomposites by examining cases where 
several GS were intercalated (polymer chains were “sandwiched” between two conse- 
cutive GS) or stacked (in the form of graphite-like structure) in the sample. By investi- 
gating their mechanical behavior in a series of uniaxial tension experiments, they 
concluded that intercalated systems exhibit a higher storage modulus compared to 
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stacked ones due to stronger interactions between the polymer and the GS resulting 
from the higher surface-to-volume ratio in the intercalated design. More recently [91], 
atomistic MD simulations were combined with advanced chemistry techniques to pro- 
duce graphene-based PMMA solutions and nanocomposites of high degree of GS dis- 
persion. Instead of graphene covalent functionalization or incorporation of small 
molecules as dispersing agents in the solutions, in the new method a fraction of PMMA 
chains were covalently functionalized by a procedure that allows adding pyrene mole- 
cules preferentially to the two ends of the polymer chain [91]. The accompanying MD 
simulations provided very useful information about the interaction of these functional- 
ized PMMA chains with GS and how these are kept well separated in the polymer ma- 
trix, thereby achieving a stable, highly homogeneous dispersion. 


5.2 Simulation methods at the atomistic level: 
the molecular dynamics method 


In the last few decades, molecular simulations have emerged as an excellent tool 
for explaining many of the microscopic mechanisms behind the measured macro- 
scopic properties of complex materials, also for connecting the predictions of theo- 
retical models with experimentally measured properties. Being able, in particular, 
to predict important material properties from the chemical composition and molec- 
ular architecture of the chemical formulation is of paramount importance because 
it can guide design efforts both at the level of materials synthesis and at the level of 
materials processing. With their inherent potential to predict important physical 
properties directly from the underlying atomistic structure, molecular simulations 
are considered today as very useful, virtual experiments that can replace in many 
cases actual laboratory measurements. 

Molecular simulations encompass two main techniques: MD and Monte Carlo 
(MC) [92-93]. MD is based on the solution of Newton's equations of motion in an 
appropriate statistical ensemble. These are numerically integrated, and the result of 
integration yields the positions and velocities of each atomistic unit in the system 
in time (this is also known as the trajectory of the system in phase space). In this 
way, we can monitor how the model system involves in time (under the macro- 
scopic constraints imposed by the statistical ensemble in which the simulation is 
carried out) and thus extract information about its thermodynamic, structural, con- 
formational, dynamic and rheological properties. To run an MD simulation, we 
need initial conditions, namely atomic positions and atomic velocities at zero time. 
The trajectory followed by the system depends then crucially on the potential en- 
ergy function describing intra- and intermolecular interactions. 

In contrast to MD which is a deterministic method, MC is based on the design 
of (artificial or even unphysical) moves to sample new states that are selected 
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randomly with the help of an appropriate acceptance criterion. Thanks, in particu- 
lar, to the design of some very clever moves, MC helps the system tunnel through 
large potential energy barriers, and this can accelerate system equilibration by or- 
ders of magnitude compared to MD. The major drawback of the method is that it 
does not offer any dynamic information about the system. This happens because 
the system evolves stochastically in configuration space based on predefined accep- 
tance criterion for each attempted move. Indeed, in a MC simulation, only atomic 
positions matter. Our discussion in the remainder of this chapter is devoted solely 
to the MD method. 

In general, the complexity of polymers prohibits the analytical solution of the 
corresponding statisticomechanical problem. Molecular simulations overcome this 
by solving the corresponding problem numerically, given a mathematical model for 
the molecular geometry and description of interatomic interactions. Two general 
types of force fields are typically employed for the molecular simulation of poly- 
mers: Explicit atom (EA) models where every atom is considered as an independent 
interaction site or entity, and united atom (UA) models where hydrogen atoms are 
neglected by considering a relatively larger, spherically interacting particle that em- 
bodies the contributions of both the hydrogen atoms and of the atom (e.g., the car- 
bon) to which the hydrogens are bonded. EA models offer a very detailed and 
accurate description of the problem. Their main drawback is that, for the same 
number of total molecules considered in the simulation box, they involve almost 
twice as many interacting sites as the corresponding UA model; as a result, molecu- 
lar simulations with EA models are significantly more CPU time demanding than 
simulations with UA models. 

Given the type of molecular model adopted, an initial configuration for the system 
and initial velocities for all atoms (assigned according to the Maxwell-Boltzmann dis- 
tribution), what remains to be defined next is the potential energy function describing 
interatomic interactions. In general, these are distinguished between bonded and non- 
bonded. Bonded interactions typically include contributions from 

— bond length stretching (Figure 5.2a), 

—- bond angle bending (Figure 5.2b), 

— proper dihedral angles (Figure 5.2c) and 
— improper dihedral angles (Figure 5.2d). 


Nonbonded interactions, on the other hand, include 
— inter- and intramolecular van der Waals interactions (Figure 5.2e) and 
— inter- and intramolecular Coulomb interactions (Figure 5.2f). 
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Figure 5.2: Schematic representations of (a) bond length stretching, (b) bond angle bending, 
(c) proper dihedral, (d) improper dihedral, (e) van der Waals and (f) Coulomb interactions to the 
potential energy of the system. 


A typical mathematical form of the total potential energy function [94] is 


Bonded energetic terms 
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Nonbonded energetic terms (5.1) 


where the first and second sums on the right-hand side correspond to a harmonic po- 
tential for the description of bond length stretching and bond angle bending contribu- 
tions, the third term is a cosine function describing contributions from proper dihedral 
angles (formed by four consecutive atoms along the chain) and the fourth term is a 
harmonic function describing contributions due to improper dihedral (or out-of-plane) 
angles formed by four nonconsecutive atoms along the chain (they help maintain the 
stereoregularity of the fourth atom with respect to the plane of the other three atoms at 
branch points). The fifth term in eq. (5.1) corresponds to the total van der Waals energy 
of the system based on a typical 12-6 Lennard-Jones (LJ) potential while the last term 
is the total electrostatic energy due to Coulomb interactions between charged atoms. 
Force fields described by eq. (5.1) are called classical force fields. Some exam- 
ples are the DREIDING [94], the OPLS [95], the CHARMM [96] and the AMBER [97]. 
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The functional form of eq. (5.1) includes the minimum number of terms (bonded 
and nonbonded), providing a satisfactory description of the total potential energy 
of the system. However, during the last 30 years, more detailed atomistic-level force 
fields have been developed based on ab initio calculations parameterized on the 
basis of a large body of experimental data for several organic compounds consisting 
of carbon, oxygen, hydrogen, nitrogen and sulfur atoms. These second-generation 
potentials are called class II potentials and differ from the classical potentials in 
that they include extra (coupling) terms between the various contributions in eq. 
(5.1), thus affording a more accurate representation of the total potential energy. A 
typical expression reads [98] 
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In general, simulations with class II potentials yield excellent predictions for the 
majority of the physical properties of system. Their only drawback is that they are 
computationally more intensive compared to the simpler, classical force fields. 
Typical examples of class II potentials are the COMPASS [98], the CFF [99] and the 
PCFF. 

Having available a mathematical expression for the calculation of the potential 
energy of the system under study, the next step is to integrate Newton’s equations of 
motion in the relevant statistical ensemble to sample system configurations in phase 
space. From this point of view, MD simulations are in many aspects similar to real ex- 
periments. When a real experiment is performed, a sample of the material is prepared 
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and connected to the instrument (a thermometer, a manometer, a viscometer, a spec- 
trometer, etc.) to measure the desired property. The value of a specific property is 
computed as the average over many different measurements to increase statistical 
accuracy. In a similar way, in an MD simulation of N interacting atoms described 
by the potential energy function Upot, the solution of Newton’s equations of mo- 
tion provides snapshots of the system in real time, from the analysis of which one 
can extract predictions of the relevant physical properties. Newton's equations of 
motion read 


mir; -Fi i- 1, 2, oN (5.3) 
or, equivalently, 
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where i denotes any atom in the system, m; its mass, r; its position vector, F; the 
force acting on it, and t the time. Solving the 3N (where N stands for the total num- 
ber of atomistic units) second-order differential equations (eq. (5.4)), we track the 
time evolution of the system and obtain its trajectory, that is, the positions and ve- 
locities of all its atomistic units in time. Applying then basic principles of statistical 
mechanics [100], we can obtain estimates of several properties of interest (e.g., ther- 
modynamic, structural, conformational and dynamic), which can be directly com- 
pared to experimental data. Figure 5.3 provides the basic structure of a typical MD 
algorithm. 

To solve the equations of motion, a numerical integrator must be employed, which 
should guarantee efficiency, stability and accuracy. To this, the algorithm must: 

- not require an expensively large number of force evaluations per time step 
(thus, very popular techniques for the solution of ordinary differential equa- 
tions such as the fourth-order Runge-Kutta method become inappropriate); 

— permit the use of a large time step; 

— be fast; 

- require little memory; 

— satisfy the energy conservation law. 


Over the years, several algorithms have been proposed that satisfy the above criteria. 
The most popular are the gear predictor-corrector methods [92], the Verlet algorithms 
[93, 101] and the multiple time step methods (such as the reversible REference System 
Propagator Algorithm or rRESPA [102]). In addition, thanks to the rapid growth of 
computing power, the development of robust parallelization techniques (based on 
domain [103-105] and force [106-107] decomposition and CUDA programming) and 
the availability of supercomputing infrastructures (computer machines with millions 
of CPU processors), highly complicated systems containing hundreds of thousands or 
even millions of interacting particles (corresponding to simulation boxes with edge 
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Figure 5.3: Simplified flow diagram of a typical MD algorithm. 


lengths equal to hundreds of nanometers) can be simulated today in full atomistic 
detail for times up to several microseconds. Also, highly parallelized, user-friendly MD 
simulation packages are commercially available or freely distributed; we mention: 

- Large-scale Atomic/Molecular Massively Parallel Simulator, LAMMPS [107] 

- GROningen MAchine for Chemical Simulations, GROMACS [108] 

— Assisted Model Building with Energy Refinement, AMBER 

- NAnoscale Molecular Dynamics, NAMD [109] 

— Chemistry at HArvard Macromolecular Mechanics, CHARMM 


5.3 Atomistic MD simulation of graphene-based 
PMMA nanocomposites 


As a model polymer, we have chosen PMMA whose nanocomposites with GS 
have been studied in detail over the years with several techniques. For example, 
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Ramanathan et al. [39] used sonication to break rigid nanoplatelets of expanded 
graphite apart into thinner platelets, which were dispersed next in a PMMA solu- 
tion using high-speed shearing methods. This led to an increase in the T, by 29 °C 
at 0.05 wt.% loading of the matrix in FGS, and up to an 80% enhancement of the 
Young's modulus at 1 wt.% loading in FGS. Similar observations have been re- 
ported by Li and McKenna [45] for GO/PMMA nanocomposites. One reason for the 
extraordinary mechanical properties of FGS-PMMA nanocomposites is the en- 
hanced interfacial interactions of oxygen functionalities across the surface of gra- 
phene with PMMA chains. FGS contain pendant hydroxyl groups, which may form 
hydrogen bonds with the carbonyl groups of PMMA. Additional enhancement 
comes from the nanoscale surface roughness of FGS, the defects caused during 
thermal exfoliation of the precursor graphite oxide and their wrinkled topology at 
the nanoscale due to their extremely small thickness. These can enhance mechan- 
ical interlocking with the polymer chains, which also leads to better adhesion. 

Atomistic MD simulations [62-63] have shown strong adhesion of PMMA chains 
(especially of its side groups) on graphene and considerably slower segmental and 
chain mobility in the interfacial area. According to simulation data, local mass den- 
sity, segmental dynamics and chain terminal relaxation differ from the bulk behav- 
ior up to several nanometers from the graphene surface. 

In the remaining of this chapter, we will focus on a methodology [110], initially 
proposed for a simpler class of systems (glassy vinyl polymers such as polypropylene 
and polystyrene) [110], which allows the determination of the mechanical properties of 
PMMA nanocomposites filled with GS (functionalized or nonfunctionalized), based on 
small-strain deformation experiments on the computer of microscopically detailed 
model structures. The procedure involves several modeling and mathematical steps 
and allows computing the elastic constants (Young's modulus E, bulk modulus B, 
shear modulus G and Poisson's ratio v) of a polymeric glass under the assumption that 
vibrational contributions of the hard degrees of freedom are not significant; as a result, 
estimates of the elastic constants can be obtained by computing changes only in the 
total potential energy of static microscopic structures subjected to simple deformation 
modes. For glassy atactic polypropylene for which the method was first developed and 
implemented by Theodorou and Suter [110], elastic constants were predicted within 
1596 of the experimentally measured values. 

All the simulations have been performed with an all-atom force field, allowing 
for a direct comparison of the computed with available experimental data. We 
chose DREIDING [95] because it combines simplicity with accuracy (for acrylic poly- 
mers). Since DREIDING [95] does not provide information about the values of partial 
charges of PMMA atoms, these were borrowed by the OPLS-AA [96] force field. 
Additional technical details (such as the parameter values of all bonded and non- 
bonded interactions describing intra- and interatomic contributions to potential en- 
ergy) can be found in two published articles [72, 91]. Figure 5.4a-d provides typical 
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Figure 5.4: Typical atomistic structures of (a) an SPMMA chain, (b) a nonfunctionalized graphene 
sheet (GS), (c) a functionalized graphene sheet (FGS) and (d) a functionalized (py-sPMMA-py) chain. 


atomistic structures of a PMMA chain, a GS, a GO and a functionalized PMMA chain 
with pyrene groups added to its two ends. 


Systems simulated and simulation strategy 


We focus on sPMMA, at P = 1 atm. The simulations were performed with strictly 
monodisperse samples with the model system consisting of 27 chains of degree of po- 
lymerization X = 15 (corresponding to a molecular weight of 1,503.75 g/mol). 
Unfunctionalized and functionalized GS had lateral dimensions 12 A x 12 A. Three 
model systems were studied: (a) the neat sPMMA matrix (no GS added; it will be de- 
noted as sPMMA in the following), (b) its nanocomposite with three unfunctionalized 
monolayer GS (it will be denoted as GS-sPMMA in the following) corresponding to 
5.67 wt.% concentration in GS and (c) its nanocomposite with three functionalized 
monolayer GS or GO (it will be denoted as FGS-sPMMA in the following) correspond- 
ing to 6.54 wt.96 concentration in GO. The surface concentration of GO in hydroxyl 
(-OH) and epoxy (-O-) groups in the latter system was chosen to match as closely as 
possible the experimentally determined concentration reported by Ramanathan et al. 
[39] through elemental analysis. 
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To build initial configurations of all systems we used MAPS [111] and to execute 
the MD simulations we used LAMMPS [107]. All initial configurations were subjected 
to static structure optimization using a molecular mechanics algorithm to remove 
overlaps, and the resulting minimum potential energy structures were annealed to 
500 K for several hundreds of nanoseconds to render them completely amorphous 
prior to quenching them down to room temperature, also to completely equilibrate 
them at all length scales. We used rectangular parallelepiped simulation cells of ini- 
tial sides 40 A x 40 A x 40 A subject to full periodic boundary conditions. Technical 
details regarding the execution of the MD simulations (type of thermostat-barostat 
used, calculation of electrostatic interactions, calculation of LJ interactions and of the 
tail corrections, integration of equations of motion, time step, etc.) can be found in 
the two relevant publications [72, 91]. 


Structural and conformational properties 


From the equilibration runs at T = 500 K, we calculated several properties that pro- 
vided a good picture of the effect of GS and GO on the structural, conformational and 
thermodynamic properties of the polymer matrix. Figure 5.5 shows a typical atomistic 


Figure 5.5: Typical atomistic snapshot from the simulation with the FGS-sPMMA nanocomposite at 
6.54 wt.% loading. The simulation cell contains 27 PMMA chains and three FGS with five hydroxyl 
groups and three oxygen atoms on their surface. Initial cell dimensions 40 À x 40 À x 40 À. Carbon 
(SPMMA), carbon (FGS) and oxygen atoms are represented with gray, yellow and red colors, 
respectively. Hydrogen atoms have been omitted for clarity. 


5 Molecular dynamics simulations of graphene-based polymer nanocomposites —— 135 


configuration of the FGS-sPMMA nanocomposite at the end of the MD simulation with 
this system at T = 500 K. A first quantity that can be easily calculated from an MD 
simulation in the NPT ensemble is the density p. Our predictions are p = 1.065 g/cm? 
for the sPMMA, p = 1.082 g/cm? for the GS-sPMMA and p = 1.091 g/cm? for the FGS- 
sPMMA system. The experimentally determined value for infinite molecular weight 
PMMA at the same temperature (500 K) is 1.072 g/cm? [112]. 

We also calculated (see Figure 5.6) the variation of polymer mass density with 
distance from a GS. 
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Figure 5.6: Local mass density normal to graphene sheets as obtained from the present MD 
simulations (T= 500 K, P = 1 atm) for the cases of pristine graphene sheets (a) and graphene oxide 
(b). The blue dashed line at the zero value of the horizontal axis denotes the average position of 
the graphene sheet midplane. 


The high-density values (up to ~30% compared to the bulk density value) observed 
at distances up to ~5 A from the surface of the GS indicate that sPMMA chains ad- 
sorb strongly on GS (either unfunctionalized GS or GO sheets). An example of a typ- 
ical conformation of an adsorbed sPMMA chain on a GO is displayed in Figure 5.7. 
An interesting point is that the local mass density of SPMMA is enhanced less in 
the FGS-sPMMA system than in the GS-sPMMA system. This should be attributed to 
the relative roughness of GO sheets (compared to the perfectly smooth surfaces of 
pristine GS) because of the presence of the characteristic -O- and -OH groups, which 
leads to the adsorption of less SPMMA molecules on GO than on GS. However, as 
we will see below, the -OH groups present on the surface of GO help the system 
develop a non-negligible number of hydrogen bonds with the oxygen atoms of 
SPMMA ester branches; this will be shown to have a strong impact on the overall 
mechanical performance of the nanocomposite. Figure 5.8 shows an example of a 
GO NP, which has developed two hydrogen bonds with one sPMMA chain on its one 
side and one hydrogen bond with another chain on its other side. The hydrogen 
bonds are highlighted with dashed circles in the figure. We clarify that the formation 
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Figure 5.7: Atomistic snapshot from the MD simulation at T= 500 K showing a typical SPMMA chain 
configuration next to a GO. 


of hydrogen bonds is not imposed directly in our simulations, but it is the indirect 
result of the employed force field (particularly of the partial charges assigned to the 
various atoms). 


Mechanical properties 


For the estimation of the mechanical properties of the simulated systems, we followed 
the methodology first proposed by Theodorou and Suter [110] for an amorphous glassy 
polymer, to which the reader is kindly referred for more details. The method involves 
the selection of several (about 15) completely equilibrated configurations of the sys- 
tem, which are then submitted to deformation experiments, from which one can cal- 
culate in a rigorous way the elastic properties of the sample. At the temperature and 
pressure conditions of interest here (T = 300 K and P = 1 atm), sPMMA and its GS- or 
GO-nanocomposites are in the glassy state, implying that one cannot directly use MD 
to sample well-equilibrated system configurations because the simulations will be 
nonergodic (the system will be trapped in configurations characterized by a local min- 
imum in their potential energy). One way to overcome this is to equilibrate the system 
at a higher temperature (above the melting point), where equilibration is much easier 
to achieve, select a good number of relaxed configurations from this simulation and 
subject them to cooling runs down to the lower temperature (T = 300 K), followed by 
a short MD run for the density and local structure to equilibrate further. The resulting 
glassy structures will then be good candidates to use in the subsequent computational 
deformation experiments for the estimation of the elastic constants (Young's, bulk, 
shear moduli and Poisson's ratio). The results can be significantly improved by aver- 
aging over all configurations subjected to deformation. 
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Figure 5.8: An example of a situation where three hydrogen bonds develop between a GO and the 
surrounding sPMMA chains. We can observe the formation of two hydrogen bonds with the same 
chain on the one side of the GO (a) and of one hydrogen bond with a different chain on the other 

side of the GO (b). 


Applying the methodology [110] to the systems studied here, we obtained the 
results shown in Table 5.1. From the numerical data presented in Table 5.1, we can 
draw several conclusions: 

1 The predicted values of the Young's modulus, shear modulus, bulk modulus 
and Poisson’s ratio for the pure sPMMA system are in an excellent quantitative 
agreement with reported experimental values, which are summarized in ref. 72. 

2) Asignificant enhancement of the mechanical properties of both types of nano- 
composites (GS-PMMA and FGS-PMMA) is observed. This is more pronounced 
for the FGS-based ones, which should be attributed to the development of hy- 
drogen bonds between the filler (GO) and the polar chains (SPMMA). 
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Table 5.1: Predicted elastic constants for all simulated systems from our study. 


System Lamé constants Elastic constants 

p (GPa) A (GPa) E (GPa) B (GPa) G (GPa) v 
Pure SPMMA 1.3 3.4 3.4 4.2 1.3 0.36 
Experimental values = - 2.07-6.2 5.9-6.17 1.78-2.3 0.34-0.36 
GS-sPMMA 1.4 3.7 3.8 4.7 1.4 0.36 


FGS-sPMMA 2.3 3.5 6.0 5.0 2.3 0.31 


To further appreciate the effect of graphene and GO on the mechanical rein- 
forcement of PMMA, we have normalized the predicted values of the four elastic 
constants with the values corresponding to the pure sSPMMA matrix, and the results 
are shown in Figure 5.9. 
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Figure 5.9: Summary of property improvement for the elastic constants of GS-sPMMA and FGS- 
sPMMA nanocomposites (T = 300 K, P = 1 atm). Numerical values have been normalized with the 
corresponding value of the neat PMMA matrix at the same conditions (E = 3.4 GPa, B = 4.2 GPa, 
G = 1.3 GPa and v= 0.36). 


Overall, our results are in a good qualitative agreement with the experimental work of 
Ramanathan and coworkers [39] who reported an increase in the Young's modulus of 
~80% in PMMA samples modified with FGS. However, a direct comparison of our sim- 
ulation results with the experimental data is difficult to make because the molecular 
weight of the polymer matrix and the size of graphene flakes used in the experimental 
measurements are too large to address with atomistic MD simulations. 
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Graphene agglomeration in the PMMA matrix and how to prevent it 


The homogeneous dispersion of GS in a polymer matrix is perhaps the most important 
and challenging issue in the fabrication of graphene-based polymer nanocomposites. 
In principle, GS fine dispersion offers larger areas for the effective adsorption 
of polymer chains, thus also for the better improvement of the properties of the 
nanocomposites. 

All calculations reported in the previous section were carried out with model 
systems characterized by a uniform distribution of GS in the polymer matrix. 
However, several times during the MD equilibrations at T = 500 K, GS agglomera- 
tion was observed to occur at long times, as GS exhibited a strong tendency to come 
close to each other and form graphitic (n-n stacking) structures. To gain a better 
understanding of such a phenomenon, we conducted an additional simulation 
study with a larger system containing 100 atactic PMMA (aPMMA) chains with de- 
gree of polymerization X = 30 and six GS of size 60 A x 60 A (system 1), and we 
monitored the time evolution of the positions of the six GS inside the simulation 
box. The results are shown in Figure 5.10a, where aPMMA chains have been omitted 
for clarity. We see that already from the first two nanoseconds of the simulation, a 
pair of GS have come close to each other to form an agglomerate. 

To overcome the problem of GS agglomeration, we proposed [91] a novel method- 
ology, which relies on the functionalization not of GS but of a good fraction of PMMA 
matrix chains by adding pyrene groups to their ends. The functionalized chains are 
noted as py-PMMA-py (see Figure 5.4d). The key idea is that pyrene groups adsorb 
strongly on the surface of GS due to very favorable n-n interactions developing be- 
tween their four benzene rings and the corresponding hexagonal structures of GS; 
then, the intervening polymer mass between the two GS prevents them from ap- 
proaching one or the other, thus self-assembly is avoided. To test the idea, we re- 
peated the simulation of system 1 at the same temperature and pressure conditions 
(T = 500 K and P = 1 atm) by replacing approximately 40% of the aPMMA chains in 
the matrix with functionalized py-PMMA-py chains (system 2). The time evolution of 
the positions of GS for such a system is shown in Figure 5.10b. To enable a one-to- 
one comparison between the two systems (system 1 and system 2), the initial posi- 
tions and initial velocities of all atoms in GS in the simulation with system 2 matched 
exactly those in the simulation with system 1. Then, according to Figure 5.10b, for the 
entire duration of the simulation (~500 ns), the six GS remained homogeneously dis- 
persed in the polymer matrix. It is also clear that many of the pyrene groups of the 
functionalized py-PMMA-py chains were adsorbed strongly on the surface of GS. 

To shed additional light in our MD findings, we calculated the time evolution of 
the distances between the centers of mass of all GS pairs in the two simulated sys- 
tems. The results are shown in Figure 5.11a and b for system 1 and system 2, respec- 
tively. For system 1, we see that the distances between GS in pairs 1-4, 2-3 and 5-6 
suddenly dropped down to 3.4 A at t = 3, 10 and 120 ns, respectively, which is the 
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Figure 5.10: (a) Evolution of GS self-assembly during the MD simulation with the aPMMA-GS 
nanocomposite in the absence of py-PMMA-py chains. (b) Same as with (a) but for the case where 
4096 of aPMMA chains have been replaced by functionalized py-PMMA-py chains. 


characteristic distance between consecutive graphene planes in a typical graphitic 
structure, and remained to 3.4 À throughout the simulation. In contrast, for system 2, 
no such phenomenon was observed. The reason for this behavior is the strong adsorp- 
tion of functionalized py-PMMA-py chains on the surface of GS by their end-pyrene 
groups. A characteristic example is depicted in Figure 5.12, showing a GS on the surface 
of which 10 pyrene groups from different py-PMMA-py chains have been adsorbed. 
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Figure 5.11: Time evolution of the distances between the centers of mass of graphene sheets in all 
GS pairs in the course of the MD simulation with (a) system 1 and (b) system 2. 
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Figure 5.12: Multiple pyrene adsorption on the two faces of a GS. We can discern the parallel 
arrangement of pyrene groups on the surface of the GS due to n—m stacking, with the rest of the 
polymer chain extending away from the GS. Hydrogen atoms have been omitted for clarity and 
adsorbed pyrene groups are depicted in yellow. 


A more detailed analysis revealed four different types of adsorbed structures: 

a) py-PMMA-py chains adsorbed on a GS by only one of their two ends (the other 
end remains free) forming a dangling end. A typical example is shown in 
Figure 5.13a. 

b) py-PMMA-py chains adsorbed on the same face of the same GS by both of their 
end-pyrene groups forming a loop. A typical example is shown in Figure 5.13b. 

c) py-PMMA-py chains adsorbed on both faces of the same GS forming an ex- 
tended loop. A typical example is shown in Figure 5.13c. 

d) py-PMMA-py chains adsorbed with their end-pyrene groups on two different GS 
forming a bridge. A typical example is shown in Figure 5.13d. 


The time evolution of the population of the four different types of adsorbed confor- 
mations in the course of the MD simulation is depicted in Figure 5.14. We see that, 
initially, all py-PMMA-py chains are free (i.e., nonadsorbed). This was done on pur- 
pose since we did not want to bias the py-PMMA-py chains to be close to GS. 

Then, during the first 15 ns of the simulation with system 2, the population of 
dangling ends increases rapidly, and this continues (albeit with a smoother rate) at 
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Figure 5.13: Examples of typical py-PMMA-py conformations adsorbed on GS: (a) a dangling end, 
(b) a loop, (c) an extended loop and (d) a bridge. 
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Figure 5.14: Time evolution of the relative population of adsorbed (dangling ends, loops, extended 
loops and bridges) and nonadsorbed (i.e., free) py-PMMA-py conformations in the course of the MD 
simulation. 
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Figure 5.15: Typical conformations of adsorbed and free py-PMMA-py chains on GS (in yellow) in 
the simulation cell after several nanoseconds of simulation time. Different conformations are 
shown in different colors: closed and extended closed loops in blue, bridges in red, dangling ends 
in green and free chains in white. 


all later times. Also increasing with time are the numbers of loops, extended loops 
and bridges. At the same time, the fraction of nonadsorbed chains decreases. 

Figure 5.15 presents the configuration of system 2 at the end of the (approxi- 
mately 500 ns long) MD simulation. Only graphene, pyrene and backbone atoms of 
the py-PMMA-py chains are shown for clarity. With blue, red and green colors, we 
depict the py-PMMA-py chains that participate in loops, bridges and dangling ends, 
with white color we depict free py-PMMA-py chains, and with yellow and black col- 
ors we depict carbon atoms in GS and pyrenes, respectively. At the end of the simu- 
lation, only few py-PMMA-py chains have remained free (i.e., nonadsorbed); the 
rest have formed bridges, extended bridges and dangling ends. 


5.4 Conclusions 


Within the space limitations of this chapter we hope that we have given a good ac- 
count of the unique, extraordinary properties of graphene that have justified why it 
is widely considered as the material of the future. We have highlighted its great 
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potential for new applications and have presented the recent progress in the design 
of graphene-based polymer nanocomposites both from an experimental and a 
computational point of view. 

We have also presented the basic concepts of molecular simulations focusing 
particularly on the MD technique, and simulation results for the structural, confor- 
mational and mechanical properties of a test nanocomposite system based on PMMA 
filled with two different types of GS: simple (i.e., unfunctionalized) GS and function- 
alized (i.e., GO). Our study was driven by recent experimental reports according to 
which addition of a small fraction of GO in a polar polymer matrix such as PMMA can 
lead to a remarkable enhancement of its elastic constants. 

We started with simulations at a high temperature (T = 500 K) well above the 
melting point of PMMA, where system equilibration was easier to achieve and studied 
several physicochemical properties of the corresponding nanocomposites (density, 
polymer dimensions, local structure, type and statistics of adsorbed conformations). 
We found that for both types of fillers, the polymer local density in the proximity of 
GS is significantly higher than in the pure bulk. In the case of the GO-sPMMA nano- 
composite, we observed the formation of a large population of hydrogen bonds be- 
tween the ester branches of sPMMA and the hydroxyl groups of GO NPs. 

In a second step, a good number of selected configurations from these high- 
temperature simulations were cooled down to room temperature, re-equilibrated and 
subjected to a rigorous statisticomechanical analysis which provided predictions of 
their elastic constants. For the pure glassy sPMMA matrix, these were found to be in 
very favorable agreement with values measured in the experiment. For the nanocom- 
posites, we found that GO sheets, even when present in very small concentrations in 
the sPMMA matrix, lead to a remarkable enhancement of the mechanical perfor- 
mance. This is due to the surface chemistry of GO which promotes strong interfacial 
interactions (mainly through hydrogen bonds) with the polymer matrix, thus signifi- 
cantly increasing the cohesive energy of the system and therefore its mechanical 
strength. Unfunctionalized GS also enhance mechanical behavior but to a much 
lesser degree than GO. 

We also highlighted the tendency of GS to come close to each other due to very 
favorable n—n interactions and form stacked (layered) structures in the polymer ma- 
trix. To overcome this, a new methodology was discussed based on the end- 
functionalization of a good fraction of PMMA chains in the nanocomposite with pyr- 
ene functional groups to form py-PMMA-py chains. Our atomistic simulations 
showed that with this new nondestructive and thus highly attractive technique (co- 
valent functionalization of graphene is completely avoided) GS self-assembly is 
avoided. This happens because pyrene groups adsorb quickly on GS and, owing to 
the intervening polymer mass, GS are kept separated, forming a stable and uniform 
dispersion. Our simulations revealed the creation of several interesting adsorbed 
structures in such a system: (a) closed loops (formed by py-PMMA-py chains ad- 
sorbed on the same face of a GS through both of their end pyrene groups), (b) 
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extended closed loops (formed by py-PMMA-py chains adsorbed on the two faces of 
the same GS), (c) bridges (formed by py-PMMA-py chains adsorbed on two different 
GS) and (d) dangling ends (corresponding to the free ends of py-PMMA-py chains 
adsorbed on a GS with only one of their pyrene groups). All these characteristic 
structures, especially loops and bridges, keep GS apart and well separated in the 
PMMA matrix, even at very high wt.96 loadings of the matrix in GS, ensuring a ho- 
mogeneous GS dispersion in the nanocomposite. 
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6 Carbon from waste source for Li-ion battery 


6.1 Introduction 


Energy generation and use depends heavily on fossil fuels, and it severely impacts 
the global economics and living sustainability [1]. Completely shifting to clean en- 
ergy production and consumption is vital for the present and future generations. 
Electrochemical energy storage and conversion can be implemented by various de- 
vices such as (rechargeable) batteries, fuel cells and electrochemical capacitors 
(ECs). They are considered fundamental elements for developing a sustainable en- 
ergy society. In particular, a battery [2] is an electrochemical power source compris- 
ing several cells, which are connected in series or parallel or in combination in 
order to deliver the required output voltage and energy. 

In ancient Greece, Thales observed that rubbing amber could generate an elec- 
tric charge [3]. In between 250 BC and AD 640, there are evidences for the use of a 
storage system known as the “Baghdad batteries" [4]. In 1748, Benjamin Franklin 
first coined the term “battery” to describe an array of charged glass plates. In the 
1780s, Luigi Galvani demonstrated the electrical basis of nerve impulses, which is 
widely called as the frog leg experiment. The Italian physicist Alessandro Volta in 
1800 invented the chemical battery, Voltaic pile. Volta discovered that an electrical 
current is generated when metals come in contact with chemicals. John Frederic 
Daniell in 1820 improved the voltaic cell using copper (Cu) and zinc (Zn) plates 
as electrodes and sulfates of respective metals as electrolytes, and termed as the 
Daniell cell, which was later used for operating telegraphs and doorbells. In the 
1830s, Michael Faraday conducted experiments with a ferrite ring and device called 
a battery. Carbon-zinc battery was developed by Georges Leclanché in the 1860s. 
The wet cell consisted of electrodes that were plunged into a body of electrolyte 
fluid, which was manufactured easily and exhibited acceptable shelf life [5]. An im- 
proved version known as a dry cell was later developed by sealing the cell and 
changing the fluid electrolyte to a wet paste. Raymond Gaston Planté in the 1860s 
invented the first practical storage lead acid battery. He immersed two thin solid 
lead plates separated by rubber sheets in a dilute sulfuric acid solution to make a 
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secondary battery. Later in 1900, Thomas Edison developed the nickel (Ni) storage 
battery and later in 1905 he developed the Ni-iron (Fe) battery [6]. 

A cell is made of three major components: a positive electrode (cathode), a neg- 
ative electrode (anode) and a separator. A separator is soaked in organic carbonate- 
based liquid electrolyte. The electrolyte contains dissociated salts, which facilitate 
ion transfer between the two electrodes. Once these electrodes are externally con- 
nected to a device, the chemical reactions occur at the electrodes. The process gen- 
erates electrons and enables the current to flow. The terms anode and cathode are 
applicable generally for primary batteries, but are not suitable for secondary (re- 
chargeable) batteries. During charging and discharging the direction of electron 
flow changes, however, the electrodes retain their polarity. Thus, the terms negative 
and positive electrodes are most suited to describe the terminals of a rechargeable 
battery system. 

The energy storage system can be primary or secondary. Primary cells are de- 
signed to discharge once and can be discarded after use. If one attempts to recharge 
a primary cell, then it can present serious safety hazards. The examples of primary 
cells include carbon-zinc and alkaline cells. Primary cells having lithium (Li) as an 
anode are called Li cells irrespective of the nature of the cathode. Such cells are 
separated into two types: cells with solid cathodes (e.g., Li-MnO;) and cells with 
liquid cathodes (e.g., Li-SO;) [7]. Reserve batteries are a special-purpose primary 
cell, which are used for applications that require a long inactive shelf period fol- 
lowed by intense discharge and unusual operational conditions. The oldest form of 
a rechargeable battery system in use is the wet cell lead acid battery. It represents 
one of the worst energy-to-weight ratios. However, it delivers good power-to-weight 
ratio along with good energy-to-volume ratio compared to other type of batteries. It 
is cheap and can supply high surge currents needed in starter motors along with 
low manufacturing cost. 

Li metal [8] has been selected as a good choice due to its low mass density, 
smaller size compared to other known battery materials as charge carrier, and a 
high mobility and fast transport speed. These features give Li-based secondary bat- 
teries a distinct advantage over other battery systems such as Ni-cadmium (Cd) and 
lead acid batteries. Cd is highly toxic and heavier, and the same problems are ap- 
plied to lead as well. The discovery (1970s) of several intercalating inorganic com- 
pounds was monumental in the development and widespread use of high-energy 
rechargeable Li-based batteries. Whittingham et al. introduced the concept of elec- 
trochemical intercalation and its potential use [9]. The discovery and the under- 
standing of the intercalating compounds led to the development of rechargeable Li 
batteries using Li insertion compounds as positive electrodes and Li metal as the 
negative electrode (Li metal battery). The first Li metal battery cell was introduced 
by Exxon et al. [10] in which they used titanium disulfide (TiS;) as positive elec- 
trode, Li metal as negative electrode and Li perchlorate in dioxolane as the electro- 
lyte. The safety concerns related to the use of Li metal as anode was addressed via 
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alternative approaches. It led to the modification of either the electrode or the elec- 
trolyte. Murphy et al. [11] and Scrosati et al. [12] demonstrated that Li metal can be 
replaced by an insertion host having low standard reduction potential. In the early 
1990s, the concept of using electrodes having Li*-ion insertion capabilities led to 
the development of today's Li'-ion batteries. The presence of Li in the anode in 
ionic form rather than in metallic form of a Li*-ion cell resolved the dendrite prob- 
lem, which was inherited by Li metal cells. 

In the 1980s, Goodenough et al. [13, 14] proposed the high-voltage metal oxides 
(Li,MO,, where M = Ni, cobalt (Co) or manganese (Mn)) as cathodes for Li*-ion cells. 
The discovery of highly reversible, low-voltage Li*-ion intercalation/deintercalation 
process in graphite led to the creation of the first Li*-ion cell based on C/LiCoO;, 
which was commercialized by Sony Corporation. Such cells exhibited an opera- 
tional voltage of 3.6 V (3 times vs. alkaline systems), gravimetric energy densities of 
150 Wh/kg (2-3 times vs. Ni-Cd (cadmium) batteries) and a cycle life of at least 500 
discharges. 

Presently, Li*-ion batteries are the fastest growing storage devices. The salient 
features of today's commercial Li*-ion batteries are 

- high operating voltage between 3.6 and 4.2 V, 
- compact, 

- lightweight, 

— high energy density, 

— fast charging, 

— low self-discharge, 

- not suffer from memory effect, 

— use between -20 and +60 °C and 

— recharge for more than 1,000 cycles. 


The Li-ion rechargeable cell working mechanism depends on an “intercalation/de- 
intercalation" process [15]. The storage processes, which do not involve intercala- 
tion are called “insertion/de-insertion” process. Indeed, the term “intercalation” is 
used for layered-type host substances (e.g., graphite, TiS,, molybdenum disulfide 
(MoS,) or lithium cobalt oxide (LiCoO;)), while “insertion” is used for 3-D materials 
that have narrow channel cavities. The charge/discharge process is reversible 
where the Li* ions move back and forth between the anode and cathode. The pro- 
cess presented the Li*-ion battery names such as rocking chair [16], swing [17, 18] or 
shuttlecock [19]. 
In general, for a Li*-ion cell: 
- The electrode materials must possess high specific capacity and durability. 
- The electrolyte must exhibit high ionic conductivity at ambient and subam- 
bient temperatures. 
— The electrolyte must have good interfacial properties. 
— The materials should be economical, eco-friendly, easily available and safe. 
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Thus, the selection of right materials is important to produce a sustainable and an 
efficient battery pack. 
A successful electrode should meet a number of criteria: 
- The intercalation/insertion compound should have large number of available 
Li sites. 
- TheLi'ions insertion/de-insertion process should be reversible. 
—- The host material should not undergo structural changes. 
— The host material should be accessing the multiple oxidation states of the tran- 
sition metals. 
- The compounds should possess good electronic and ionic conduction. 
- The electrodes must be stable over the working voltage range. 
— The electrodes should not undergo dissolution into the electrolyte. 
- It should be environmentally friendly and easy to produce. 


Li metal is an ideal anode for rechargeable batteries due to its high potential and 
specific capacity. It is widely used as an anode in primary Li cells. However, den- 
dritic growth, cost, safety hazards and poor recharge ability prevents its practical 
use in secondary batteries. The issues are addressed by the use of Li alloys (with 
Al or Si). Unfortunately, the alloy anodes undergo severe volume changes during 
the cycling process. Nevertheless, the urge for the development of solid-state (ce- 
ramic/polymer) electrolytes are a promising approach to realize the use of safe Li 
metal rechargeable batteries. The standard and commercialized anode for Li-ion 
batteries is carbon in its different layered structures. 


6.2 Role of carbon material in batteries 


Li metal anodes with an organic electrolyte result in nonuniform formation of a pas- 
sive film on the anode surface [20-22], which causes dendrite growth of Li metal. 
Longer charging-discharging cycles are one of the important requirements for mak- 
ing a commercially viable Li-ion battery. This problem is tackled by utilizing the 
electrochemical intercalation of Li in carbon. It produces a significant negative po- 
tential close to that of Li, but this is less reactive and easily reversible [23]. 

Various kinds of carbon materials have been employed as anodes in various 
electrolytes for Li rechargeable batteries. Carbon offers advantages over other types 
of materials used in rechargeable batteries such as high operating voltage, longer 
life cycle and lack of *memory" effects [24—26]. The structural characteristics of car- 
bon are believed to be the major factor that controls the performance of Li-ion bat- 
teries. The intercalation of Li* ion with carbon depends upon many factors such as 
the preparative history, carbon precursor and structure. The variations in the aver- 
age layer spacing and the voltage profile for the carbon with intercalating depend 
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on the degree of graphitization [27]. Literature survey on carbon in battery develop- 
ment revealed that graphite is the most commonly used anode for Li*-ion second- 
ary battery. Many studies employed graphite as an anode, which has theoretical 
capacity of 372 mAh/g based on the stoichiometry of Li carbon [28-34]. 


6.3 Graphene materials for highly performing 
electrodes 


The widespread interest in graphene includes its potential as a suitable material for 
energy storage and particularly for Li'-ion batteries. Intrinsically, graphene is 
aimed at providing excellent performance based on a combination of physicochem- 
ical properties, its extremely high surface area of 2,630 m^/g [35], high electrical 
conductivity and chemical inertness [36]. 

Crystalline carbon nanomaterial, carbon nanotube [37, 38] and graphene had 
their expected applications in electrochemical energy storage. This was precisely 
envisioned for both Li*-ion and Li-sulfur batteries and supercapacitor devices [39]. 
Advantageously with respect to carbon nanotube, graphene is the basic unit of sp? 
carbons (graphitic), which are the most conventionally/commonly used carbons for 
batteries (or other high surface area applications). 

Production wise, the cost of using a certain graphene material could be compared 
with conventional and cost-efficient traditional macroscopic methods, based on top- 
down approach or on raw materials, typical chemical or metallurgical approaches to 
obtain graphite (and graphene-like materials) synthesis and cell-making methods. 
Examples of current trends and recent achievements include pyrolytic graphite from 
graphene oxide [40] and kish graphite [41], which are obtained from waste graphite 
flakes from steel production. Scalable recycling and reusing approach applied to 
these abundant raw materials will be cited in the following sections. 

In the following section, a nonexhaustive collection of bottom-up methods, se- 
lected techniques for production of high-quality graphene materials as well as ex- 
amples of graphene-decorated electrodes are described. For a comprehensive topic 
review on graphene synthesis, there exist a number of extended resources such as 
how to achieve specialized morphologies, particularly, those specifically useful/ap- 
propriate for EC energy storage [42]. Interestingly, some general vision and hints on 
using waste and recycled or recyclable carbon sources are included for the de- 
scribed methods. Besides, the quest for valorizing waste materials for energy stor- 
age is not new or exclusive for novel nanomaterials [43]. The use of hierarchical 
structures such as combining graphene with other materials, for example, oxides 
[44], is beyond the scope of this text. 
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In summary, graphene is still considered to ultimately outperform existing solu- 
tions as it would provide improvements in terms of increased surface area (i.e., en- 
abling higher capacity and also lighter devices), conductivity (faster charge- discharge 
and efficiency) and increased chemical and mechanical resistance. Relevant evalua- 
tion and benchmarking of graphene electrodes as electrochemical energy storage ele- 
ments in affordable graphite materials, to critically assess their real potential and, 
eventually, quantify the advantages of graphene in operation that is needed in both 
energy and power density are carried out. The expected exponential increase in de- 
mand and performance of batteries, resulting from their use in electrified vehicles, re- 
quires advanced and sustainable solutions. Fundamental investigations on graphitic 
structure from waste material for electrochemical energy storage are focused on eluci- 
dating whether they provide unique features and advantages to affordable carbons, 
rather than being competitive to existing commercial solutions. Therefore, they are 
generally based on the "ideal" concept of graphene. High-quality graphene material, 
in contrast to common graphite/graphitic materials and low-quality graphene materi- 
als, can be defined as pure carbon, large crystal domain size and isolated films, that 
is, typically consisting of monolayer graphene or, eventually, bilayer or few layers' gra- 
phene. In fact, classification of graphene based on these three parameters has been 
established for proper reference in scientific papers and enabling some standards for 
industrial purposes, such as production [45]. Advantageously, laser power seems to 
provide a way to transform a variety of heterogeneous materials into two graphitic 
forms [46]. Similarly, one could envision that the resources potentially include waste 
materials. 

Recently, attention has been drawn to use disordered carbon materials. It may 
store Li via a mechanism that is completely different from that associated with 
graphite. Additional studies indicate that Li insertion into an unorganized carbon 
loses the layer structure [47]. Variation in the polarization resistance with intercala- 
tion revealed the processes in various disordered carbons. It demonstrated the sin- 
gle-phase reactions with different stoichiometry of Li intercalation. Electrochemical 
intercalation of Li into carbons has been studied using mesophase pitch-based car- 
bon fibers with different heat-treatment temperatures, coke [48] and graphite as an- 
odes [27, 47, 49]. It has been found that the activation energy for the diffusion 
process decreases by increasing the degree of graphitization [27]. The studies also 
demonstrated that intercalation of Li with camphor-based carbon nanobeads stabi- 
lize the cell for 10-20 days [50, 51]. The reversible Li*-ion intercalation capacity was 
found to be 45-6196 of that obtained with graphite [51]. Researchers also reported 
the result of Li intercalation by carbon nanomaterial obtained from the pyrolysis of 
biochar such as tealeaves [52]. 

Research on Li diffusion in carbon is of both theoretical and practical interest. 
The process was largely studied in connection with the development of Li*-ion batter- 
ies, whose negative electrode is based on graphite or some other carbon material 
[27, 53-61]. In some works, Li diffusion was investigated during the reverse process, 
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that is, the anodic deintercalation [47, 62-66]. A complex impedance model for spher- 
ical particles was used to determine the Li*-ion diffusion coefficient in graphite as a 
function of the state of charge and temperature [67]. Practically, in all cases, studies 
were concentrated on the diffusion of Li, which was electrochemically incorporated 
during cathodic polarization of a carbon material in an aprotic electrolyte [68]. 

However, carbon pyrolysis at low temperature shows attractive behavior as 
their theoretical capacity exceeds that of graphite as anode material [52, 69-75]. 
This observation indicates that there may be some other different mechanisms oper- 
ating in these carbons in terms of Li intercalation. However, the density of low- 
temperature carbon is lower than that of graphite [76]. Therefore, the capacity of 
low-temperature carbon should be large enough to compensate for its smaller den- 
sity. The capacity of an anode for Li-ion battery depends on the amount of carbon 
that can be placed within the limited volume of a typical battery container. Carbon 
obtained from the same precursor at different temperatures has different structural 
orientation. Therefore, trial and error is often needed to establish a suitable prepa- 
ration method for optimizing the properties of carbon [77, 78]. Many researchers 
studied the relation between the properties of carbon and its electrochemical char- 
acteristics relevant to Li*-ion rechargeable [79-85]. 

As for the source of atomic carbon, methane is the commonly used gas as its de- 
composition products are the simplest. The reuse of greenhouse gases such as the 
very methane for graphene production would follow automatically. Similarly, green 
approach is also aimed at reusing carbon dioxide [86], yet it is more complex in 
terms of chemistry and efficiency of the reaction. Likewise, a number of additional 
works have actually been reported for the (re)use of greener carbon sources, which 
include molecules or polymers [87], as well as various recycled waste materials, such 
as plastics but also organic materials or industrial coal (Figure 6.1) [88-90]. 

In this study, we examined the effects of carbonization atmosphere on the 
charge-discharge characteristics of carbon from waste source as an anode for a Li - 
ion rechargeable battery. This work targeted the *green" approach for energy stor- 
age. The carbonaceous materials such as carbon nanowalls (CNWs) and carbon 
fiber sheet used in this study were made from different precursors such as organic 
waste (methane) and polymer waste [nanoscale microfibrillated cellulose (NMFC)]. 


6.4 Vertically oriented micro-mesoporous CNWs 


Graphene materials obtained by chemical vapor deposition (CVD) present several 
limitations for EC energy storage. First, morphologically the effective graphene sur- 
face could be reduced up to half of its theoretical value as CVD-obtained graphene 
is bonded/stuck onto the metal substrate. Besides, the high temperature used in 
CVD is not desirable to obtain wrinkled or crumbled carbon, especially if power 
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Figure 6.1: Graphene films are converted from coal tar pitch on the top and bottom surface of Ni 
layer [90]. (Image courtesy: https://www.nature.com/articles/srep16710/figures/1) 
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consumption reduction is aimed, but also for process compatibility, with any sub- 
strate such as plastics. Apart from providing better volume-area ratio, an increase in 
charge storage capacity is expected from the large amount of graphene flake edges to 
the electrolyte. Consequently, alternatives of producing porous or microstructured 
carbon materials, ideally based on scalable techniques, have been studied [91]. 

As an analogy with CNTs, this is the case of using physically assisted CVD sys- 
tems, such as fitted with radiofrequency and microwave (MW) sources. In particular, 
for certain solid sources, the reactivity of physically assisted CVD such as based on 
the use of a powerful MW field may restrict the growth process to be operated in con- 
ditions of relatively low pressure. For instrument and operator's safety in the pres- 
ence of chemical elements (carbon source), many systems and synthesis principles 
do make use of hydrogen as a carries gas to avoid oxygen contamination. With the 
plasma systems, the reuse of methane and other recycled carbon sources may be pos- 
sible to be applied. Although the compatibility of recycled carbon sources could be 
more difficult in general, it has actually been achieved [92]. Yet, more examples can 
be found in the literature [93]. 

Figure 6.2 provides the image of CNWs deposited on Cu substrate, which has 
the possibility of being used as anode for Li*-ion batteries. Owing to the orienta- 
tion and morphologies, obtained materials are typically referred in the literature 
with a variety of terms, including carbon nanosheets [94] or CNWs [95]. These 
graphene flakes or CNWs are randomly distributed, perpendicularly oriented to 
the substrate and/or curvy. Their dimensions are maximum of approximately 
1-1.5 pm in height and 10-20 nm in thickness, which accounts for multilayered 
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Figure 6.2: Microscopic images of CNWs deposited on Cu: (left and center) FESEM and (right) TEM 
image. 


graphene, and are characterized by very small crystal size, in the order of nanometers. 
This is further confirmed with micro-Raman analysis. 

The observed Raman spectra (Figure 6.3) is typical of a good quality graphitic 
material. Both first- and second-order Raman effects can be seen. First-order Raman 
peak is presented as a so-called G-peak around 1,580 cm“. This peak is commonly 
observed in all graphitic materials being often the sharpest and intense one. It 
arises from the first Brillouin zone center of the graphene unit cell and it is due to 
the double degeneration of in-plane transverse optic [96, 97] and longitudinal optic 
modes [98, 99]. All other peaks are due to second-order Raman effects. The main 
peak belonging to this category is D-peak around 1,350 cm‘. This peak is normally 
associated with the disorder of the graphitic structure. The principal reasons for the 
presence of this peak are the combination of grain boundaries, vacancies and other 
structural defects. The D-peak is particularly intense and very sharp: this can be ex- 
plained not mainly by the density defects but related to the particular orientation of 
CNWs, which are vertically aligned or oriented (mainly edges exposed to laser 
beam). This statement is further confirmed by looking at the D + G peak. It does not 
show any abnormally great number of defects. In conclusion, D-peak intensity is 
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Figure 6.3: Raman analysis of CNWs. (Image courtesy by Electrochimica Acta (Elsevier) to reuse in a 
book under the copyright license number 4454750182186) 


correlated with the number of defects in the structure, while the full-width half- 
maxima of D-peak is correlated with the grain size. The presence of a right- 
shouldered G-peak around 1,620 cm™ is called second-order D'-peak which appears 
as a shoulder peak of G-peak. This peak is observed in some disordered graphitic 
carbon structures, including graphene, but cannot be seen in defect-free structure. 
Its intensity is usually lower as compared to the G-peak but becomes more impor- 
tant especially in CNW structure. Other second-order features are indicated by the 
remaining peaks named 2D, D + G, 2G. The peak around 2,700 cm ! is often called 
2D. Generally, this peak fall around a zone of the Raman graph where the 2D 
Raman shift doubles the Raman shift of the D-peak [100, 101]. The origin of peaks is 
not related to the presence of defects or grain boundary. It arises from the interac- 
tion of two phonons. The peak is very intense, especially in planar structure like 
graphene, CNWs. The peak shape of Raman can give us information about the num- 
ber of layers present in the structure [102]. The peak around 2,700 cm ! observed 
here is D + G peak, which is generally associated with the damaged graphene. It is 
less evident or almost absent in high-quality and defect-free carbon structure [103]. 
The peak around 3,240 cm is sometimes called 2G and is present in both planar 
(graphene, graphite, CNWs) and nonplanar (CNT) structures [104]. 


6.5 Growth of CNWs 


Vertically arranged graphene materials, arrays of CNWs, can be grown onto as- 
purchased flexible metal polycrystalline Cu foil. It is directly tested as anode in Li-ion 
batteries without any prior treatment [70, 105]. A commercial MW plasma-enhanced 
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CVD technique (CVD-CN-100 from Ulvac, Japan Ltd.) [106-108] is used for the 
growth of CNWs. The whole process and growth conditions typically consist of, 
first, exposure of Cu foil to a hydrogen (H3) flux; then, methane (CH,) as a reaction 
gas is introduced in combination with H- (H2:CH,) with the ratio of 20:20 sccm. At 
growth stage, bias voltage of 200 V is applied. The whole process is done at a con- 
stant pressure of 1.5 mTorr pressure. In the process, hydrocarbon gas decomposes 
into graphitic carbon. The process needs to be catalyzed and is therefore applied on 
affordable Cu-metal surface, which provides a better control on the number of layers 
with electrochemically active surface. Figure 6.4 shows the charge-discharge cycling 
of CNWs on Cu at different current rates. The behavior of current is slightly different 
from the standard behavior of carbon. This may be due to the presence of disordered 
carbon in the CNW deposits (as learnt by Raman spectroscopy). However, one of the 
highlights of the material is that even if the current rate is increased to higher orders, 
the capacity is not dropping down. 
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Figure 6.4: (a) Graph showing the galvanostatic charge/discharge curves of CNWs deposited on Cu: 
inset - morphological characterization of the nanostructured carbons deposited on Cu by FESEM 
(ZEISS Supra 40); (b) capacity versus number of cycles in a cell Li/electrolyte/CNW-Cu at room 
temperature. (Image courtesy by Electrochimica Acta(Elsevier) to reuse in a book under the 
copyright license number 4454750182186) 


6.6 Carbonized cellulose paper 


Use of paper in our everyday life has turned out to be an original answer to the de- 
mand of eco-friendly, flexible, lightweight and cost-effective materials for electron- 
ics [108, 109] and energy-storage devices [110, 111]. There are promising reports on 
paper/cellulose-based composite electrodes for supercapacitors and Li batteries 
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[112, 113]; the rough and porous surface of paper has turned out to be a great ad- 
vantage, where both high surface area and electrolyte absorption are in demand. 
Recyclability and biodegradability are significant benefits, considering that cost 
of waste management and impact on the environment could be considerably low. 
It has been proven by Liu et al. [109] that paper printed circuit boards have about 
two orders of magnitude less impact on the environment than the ordinary printed 
ones. 

The work on this material has demonstrated the feasibility of manufacturing 
graphite-based negative electrodes on a pilot-line exploiting microfibrillated cellu- 
lose as binder, through the spray-coating method, and the use of water-based pa- 
permaking technologies for large-scale production of Li batteries [114]. Beneventi 
et al. have successfully assembled a flexible, truly solid Li*-ion cell based on self- 
standing electrodes manufactured using conducting paper as substrate. This ap- 
proach makes use of a carbonization step of the nanocellulose-binded electrodes at 
high temperature to obtain conducting carbon nanofiber binder concurrently re- 
moving water and impurities, so that no successive drying step was necessary, thus 
making it more cost-effective in terms of energy and time. An NMFC-reinforced 
polymeric membrane which was activated by incorporating a freshly prepared 
room temperature ionic liquid-based green and safe electrolyte solution was used, 
and the resulting truly solid system assembled in a “pouch cell” envelop demon- 
strated a remarkably stable cycling behavior upon prolonged cycling at ambient 
temperature even at elevated current rates as well as in bent configuration. It was 
also noted that all cell components could be easily recovered using common water- 
based recycling procedures used in paper mill. The waste was minimal; the process 
was of low cost and low energy demanding, which led to the worldwide sustainable 
Li-based battery use and recovery. 


6.7 Preparation of self-standing electrodes 


Figure 6.5 demonstrates the different processes involved in obtaining pyrolyzed 
nanocellulose-binded electrodes, which are then used to assemble the solid poly- 
meric Li-ion cell. The anode slurry was made of graphite powder (GP), carbon 
black (CB), NMFC, carboxy methyl cellulose (CMC) and alum. The cathode slurry 
was composed of lithium iron phosphate (LiFePO,)/carbon, CB, NMFC, CMC and 
alum. Both the slurries, having a dry solid content ranging between 25 and 15 wt.%, 
were prepared by dissolving CMC in the aqueous NMFC slurry, followed by the addi- 
tion of CB. 

A commercial airless spray system (Wagner, Project Pro 119) was used to spray 
the anode slurry on both sheets and on other substrates (i.e., smooth polyester 
films and copy paper). An aqueous suspension of carbonized NMFC (2 wt.%) was 
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Figure 6.5: Different phases of the process involved in obtaining pyrolyzed nanocellulose-binded 
electrodes, which are then used to assemble the solid polymeric Li*-ion cell. (Image courtesy by 
Carbon (Elsevier) to reuse in a book under the copyright license number 4454160899455) 


prepared for reinforcing the polymer electrolyte and also the same was used as a 
binder for retaining the active materials within the electrode films (Figure 6.6). 


6.8 Pyrolysis of cellulose paper 


Pyrolysis method [115] is adopted for carbonization of waste NMFC. This involves the 
pyrolysis of self-standing anode and cathode samples at 800 and 600 °C, respectively. 
To avoid thermal degradation of the LiFePO, active material, lower pyrolysis tempera- 
ture was applied to the cathode material. The pyrolysis chamber is normally flushed 
with argon inert gas to avoid oxygen traces. The furnace is heated with 10 ?C/min. 
Pyrolysis is carried out for 30 min after achieving the required temperature. These 
steps ensured the completion of the carbonization procedure. Figure 6.7 shows the 
different forms of an NMFC sheet that has undergone pyrolysis. The pyrolysis method 
is normally followed by the cooling of the furnace to room temperature to recover the 
carbonized materials [69]. 
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Figure 6.6: Different steps illustrating the facile recovery of pyrolyzed nanocellulose-binded 
electrode materials by exploiting the standard water-based recycling procedure used in paper 
industries. (Image courtesy by Carbon (Elsevier) to reuse in a book under the copyright license 
number 4454160899455) 


E 


>» 


Figure 6.7: Self-standing ability of (a) pristine NMFC paper, (b) pyrolyzed NMFC paper, (c) pyrolyzed 
NMFC paper in its bend configuration, (d) representative nanocellulose-binded LiFePO,-based 
cathode in its top and (e) bottom side views, (f) pyrolyzed LiFePO,-based cathode in its bottom and 
(g) top side views, (h) pyrolyzed GP-based anode. (Image courtesy by Carbon (Elsevier) to reuse in a 
book under the copyright license number 4454160899455) 
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6.9 Analysis of pyrolyzed cellulose paper 


Figures 6.8 and 6.9 shows the FESEM images of GP-based anode pyrolyzed at 800 ?C 
and LiFePO,-based cathode pyrolyzed at 600 °C. The structure of the carbonized elec- 
trode (carbonized anode) is more accessible with respect to the pristine spray- 
deposited paper electrodes, and the presence of both LiFePO, and CB particles is well 
evidenced. The carbonized cellulose fiber network successfully binds the electrode 
components effectively, maintaining the preburning structure, despite slight shrink- 
age. It was observed that LiFePO, particles perfectly maintained their characteristics 
after heat treatment without suffering any kind of modification (Figure 6.9d). The 
fundamental triphasic system composed of the active material grain (LiFePO, parti- 
cle), the electronic conductivity enhancer (CB spherical particle, locally enhancing 
the electronic conductivity of LiFePO,) and the binder (electronic conductive carbon- 
ized cellulose fiber, which in turn also ameliorates the electronic conductivity of the 
system) was still present as highlighted by the yellow circle. 
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Figure 6.8: FESEM at different magnifications of GP anode pyrolyzed at 800 °C. (Image courtesy by 
Carbon (Elsevier) to reuse in a book under the copyright license number 4454160899455) 
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Figure 6.9: FESEM at different magnifications of LiFePO, cathode pyrolyzed at 600 °C. (Image 
courtesy by Carbon (Elsevier) to reuse in a book under the copyright license number 
4454160899455) 


Figure 6.10 shows the Raman spectrum of GP-based anode, LiFePO,-based cath- 
ode and NMFC paper. A-800-TS shows a good separation and sharp D, G and Go 
peaks when compared with A-800-BS. It evidences that the carbon material in sample 
A-800-TS is better graphitized than that of A-800-BS. Raman analysis of top rough 
(C-600-TS) and bottom smooth (C-600-BS) surfaces of the LiFePO,-based cathode 
shows different behavior (Figure 6.10c and d). Intensity ratios shown in Figure 6.10f 
indicate that C-600-TS has lower graphitized matrix as compared to C-600-BS. The 
Raman spectrum of pyrolyzed NMFC paper (Figure 6.10e) shows no peak intensity 
and positioning variation on either side of the paper due to the uniform characteris- 
tics of both the surfaces. In general, the Raman analysis of pyrolyzed NMFC paper 
shows similar graphitized structure as that of C-600-BS. It can be stated that A-800- 
TS and C-800-TS have lower degrees of graphitization after analyzing the Raman 
spectra. Raman analysis also confirmed that the temperature variation did not reflect 
on the variation of graphitization of the sample. Raman spectroscopy by itself does 
not provide a complete picture of carbon quality. 
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Figure 6.10: Raman spectra of (a) GP-based anode pyrolyzed at 800 ?C analyzed on top and (b) 
bottom surfaces, (c) LifePO,-based cathode pyrolyzed at 600 °C analyzed on top and (d) bottom 
surfaces, (e) NMFC paper carbonized at 600 °C. (f) Samples versus /p//g and /e:/Ip ratios. (Image 
courtesy by Carbon (Elsevier) to reuse in a book under the copyright license number 


4454160899455) 
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Figure 6.11: Different assembly stages of an all-paper-based solid polymer cell in pouch 
configuration: (a) assembly of the electrodes on the current collector grids, (b) packing followed by 
sealing, (c) final Li*-ion cell connected to the instrument for testing. (Image courtesy by Carbon 
(Elsevier) to reuse in a book under the copyright license number 4454160899455) 


Figure 6.11 provides the different assembly stages of an all-paper-based solid 
polymer cell in pouch configuration. The composite polymer electrolyte is sand- 
wiched between the carbonized nanocellulose-binded LiFePO,-based cathode and 
graphite-based anode paper electrodes and housed in a “pouch cell” envelop with 
Al and Cu grid current collectors (for cathode and anode, respectively). Different 
phases of the assembly are shown in Figure 6.11a-c. The cell is kept for ambient 
temperature galvanostatic charge/discharge cycling at various current rates. 
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Figure 6.12: Current charge/discharge cycling behavior in terms of specific energy versus cycle 
number at different current regimes and charge and discharge potential versus time profiles 
extracted from the 400th cycle at C/5 rate. (Image courtesy by Carbon (Elsevier) to reuse in a book 
under the copyright license number 4454160899455) 
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Figure 6.12 provides the electrochemical response of the Li*-ion polymer cell 
in terms of specific energy (Wh/kg calculated on the weight of the active cell com- 
ponents excluding current collectors and packaging) versus cycle number at dif- 
ferent rates. The potential versus time profiles (400th constant current charge and 
discharge profiles at C/5 current rate) are very well consistent with the potential 
trend of the corresponding cell with natural cellulose bonded electrodes and liq- 
uid electrolyte [116]. The potential drop between charge and discharge plateaus 
was found to be limited, which accounts for low internal resistance, thus good 
electrode/electrolyte interface. The cell showed good rate capability and high cou- 
lombic efficiency approaching 10096. Specific energy values were found to be defi- 
nitely interesting for a solid polymeric device, approaching 300 Wh/kg at 1 C and 
even exceeding 350 Wh/kg when reducing the current rate to C/5. Even more im- 
portant, the cycling response is maintained steady for at least 500 cycles. This ob- 
servation is definitely remarkable for a solid polymeric device. The exceptional 
capacity retention indicates that no contact loss occurred between the active ma- 
terial particles. The interface between active material particles and pyrolyzed 
nanocellulose binder remains intact even after over 500 charge/discharge cycles. 

Figure 6.13 provides the electrochemical stability window of the composite poly- 
mer electrolyte evaluated by linear sweep voltammetry. It was found to exceed 4.6 V 
versus Li*/Li in the anodic region. The Li plating/stripping indicates an overall good 
cathodic electrochemical stability, which is evident at around O V versus Li*/Li. 
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Figure 6.13: Arrhenius plot showing the ionic conductivity versus temperature; inset: 
electrochemical stability window of the polymer electrolyte. (Image courtesy by Carbon (Elsevier) to 
reuse in a book under the copyright license number 4454160899455) 
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Figure 6.14: Ambient temperature electrochemical characterization of the carbonized electrodes in 
Li metal cell with standard liquid electrolyte: (a) potential versus time plot of a GP-based anode 
pyrolyzed at 800 °C showing typical charge and discharge profiles at C/5; (b) specific capacity 
versus cycle number plot of the same anode at different discharge/charge current regimes; 

(c) potential versus time plot of a LiFePO,-based cathode pyrolyzed at 600 °C showing charge and 
discharge profiles at 1C; (d) specific capacity versus cycle number plot of the same cathode at 
different charge/discharge current regimes. (Image courtesy by Carbon (Elsevier) to reuse in a 
book under the copyright license number 4454160899455) 


Figure 6.14 provides the potential versus time plot of the carbonized GP-based 
anode. Figure 6.14a clearly reflects the characteristics of the electrode active material. 
It shows the typical discharge (intercalation of Li* ions) and charge (deintercalation of 
Li* ions) profiles expected for a highly crystalline graphite material, where the different 
insertion steps toward the so-called stage-1 Li-graphite intercalation compound, 
that is LiCg, are clearly visible [117]. From the analysis shown in Figure 6.14b, the 
irreversible capacity loss during the first cycle can be attributed to the expected 
side reactions with the components of the electrolyte (solid electrolyte interface 
layer formation). Typically, organic solvent molecules and salt anions are reduced 
on the active material's surface. Thus, forming insoluble Li salt that precipitates to 
form a passivating film. This is helpful in preventing further reaction between com- 
ponents of the electrolyte and graphite active material [118]. 
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The reversible specific capacity obtained at low C/10 rate is close to the theoret- 
ical value for graphite, reflecting the practical specific capacity of the commercial 
material used. This can be explained by the fact that the carbon obtained upon cel- 
lulose pyrolysis also takes part in the electrochemical process to a certain extent. 
After the initial cycles, the Coulombic efficiency rapidly increases to above 99%. 
Subsequently, it remains highly stable throughout the cycles. This indicates that 
the formed surface film lasted intact and showed excellent reversible cycling after 
surface reactions were completed. At a relatively high current of 1C, the anode can 
deliver a capacity of 150 mAh/g, which is definitely higher if compared to the re- 
sults obtained for the same pristine electrode reported by Beneventi et al. [114]. 
Thus, it can be summarized that the carbonization procedure at high temperature 
successfully enhances the electrochemical behavior of paper electrodes. This can 
be due to the improved electronic conductivity resulting from the presence of the 
carbonized network of fibers, which contributes to the overall electrode specific ca- 
pacity. The performance of carbonized LiFePO,-based cathode in lab-scale Li cell 
was outstanding in terms of both overall specific capacity and rate capability 
(Figure 6.14d). The capacity obtained at low C/10 rate was around 160 mAh/g, 
which is the specific capacity of the commercial material used. It is also interesting 
to note that after the specific carbonization procedure adopted, this kind of paper- 
like cathode is able to operate at very high 10C current regime, still providing spe- 
cific capacity values exceeding 100 mAh/g, thus allowing a higher power output 
of the final device. These results are difficult to be achieved by the corresponding 
*as-prepared" spray deposited LiFePO,-based electrodes as described by Zolin et al. 
[116]. This signifies the role of pyrolyzed nanocellulose fibers, which combine the 
functions of strong binder and conducting additive without negatively affecting cy- 
cling stability and rate performance [119]. It is worthy to note that this system does 
not show any performance decay; indeed, reducing the C-rate after more than 120 
charge/discharge cycles completely restores the specific capacity, thus achieving ex- 
cellent stability. 

The overall observation concludes that the accomplishment of the “12 principles 
of Green Chemistry" is a fundamental goal to be pursued by means of designing/opti- 
mizing environmentally “conscious” processes, materials and devices to reduce their 
environmental impact throughout all the phases of their life cycle, thus benefiting 
the economy, protecting people and achieving the real goal of a truly sustainable 
world. This goal can be achieved by exploiting carbonized nanocellulosic binders 
and water-based papermaking coupled with high-performing green composite poly- 
mer electrolyte. 
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